MST125 
Essential mathematics 2 


Computer algebra guide 


The Open 
University 


'This publication forms part of an Open University module. Details of this and other 

Open University modules can be obtained from the Student Registration and Enquiry Service, The 
Open University, PO Box 197, Milton Keynes MK7 6BJ, United Kingdom (tel. +44 (0)845 300 6090; 
email general-enquiriesQopen.ac.uk). 


Alternatively, you may visit the Open University website at www.open.ac.uk where you can learn 
more about the wide range of modules and packs offered at all levels by The Open University. 


To purchase a selection of Open University materials visit www.ouw.co.uk, or contact Open 
University Worldwide, Walton Hall, Milton Keynes MK7 6AA, United Kingdom for a brochure 
(tel. +44 (0)1908 858779; fax +44 (0)1908 858787; email ouw-customer-servicesQopen.ac.uk). 


The Open University, Walton Hall, Milton Keynes, MK7 6AA. 
First published 2015. 
Copyright © 2015 The Open University 


All rights reserved. No part of this publication may be reproduced, stored in a retrieval system, transmitted 
or utilised in any form or by any means, electronic, mechanical, photocopying, recording or otherwise, without 
written permission from the publisher or a licence from the Copyright Licensing Agency Ltd. Details of such 
licences (for reprographic reproduction) may be obtained from the Copyright Licensing Agency Ltd, Saffron 
House, 6-10 Kirby Street, London ECIN 8TS (website www.cla.co.uk). 


Open University materials may also be made available in electronic formats for use by students of the 
University. All rights, including copyright and related rights and database rights, in electronic materials and 
their contents are owned by or licensed to The Open University, or otherwise used by The Open University as 
permitted by applicable law. 


In using electronic materials and their contents you agree that your use will be solely for the purposes of 
following an Open University course of study or otherwise as licensed by The Open University or its assigns. 


Except as permitted above you undertake not to copy, store in any medium (including electronic storage or 
use in a website), distribute, transmit or retransmit, broadcast, modify or show in public such electronic 
materials in whole or in part without the prior written consent of The Open University or in accordance with 
the Copyright, Designs and Patents Act 1988. 


Edited, designed and typeset by The Open University, using the Open University TEX System. 
Printed in the United Kingdom by Latimer Trend and Company Ltd, Plymouth. 


SUP 03479 5 
1.1 


Contents 


Introduction 


1 


9 


Installing, testing and configuring Maxima 
1.1 Introducing and installing Maxima 
1.2 Starting, testing and configuring wxMaxima 


Getting started with Maxima and wxMaxima 
2.1 Calculating with Maxima 

2.2 Reusing and editing commands 

2.3 Assigning values to variables 

2.4 System variables 

2.5 Annotating your wxMaxima worksheet 

2.6 Saving and printing your work 

2.7 Getting help 


Using Maxima 

3.1 Algebra and functions 
3.2 Plotting graphs 

3.3 Solving equations 

3.4 Calculus 


Number theory 


Parametric equations 
5.1 Plotting curves from parametric equations 
5.2 Plotting conics parametrically 


Partial fractions and polynomial division 
6.1 Partial fractions 
6.2 Polynomial division 


Differential equations 

7.1 Solving differential equations analytically 
7.2 Plotting direction fields 

7.3 Solving differential equations numerically 


Eigenvalues and eigenvectors 
8.1 Matrices 
8.2 Characteristic equations, eigenvalues and eigenvectors 


Solving linear recurrence systems 


Maxima accessibility guide 


Computer methods for CAS Activities in Books A-D 


Contents 


O XI A a 


14 
14 
20 
23 
26 
2r 
28 
30 


32 
32 
3T 
45 
51 


58 


60 
60 
65 


67 
67 
69 


69 
70 
74 
79 


82 
82 
86 


90 
95 
118 


Contents 


Solutions to Computer activities 


Maxima reference guide 


Acknowledgements 


Index 


122 
161 


170 


171 


Introduction 


During this module you will use the Maxima computer algebra system 
(CAS) introduced in MST124 to investigate mathematical concepts and 
perform calculations that may be too laborious or difficult to do by hand. 


'This Guide is designed to be read in conjunction with the main study 
texts and is not meant to be read all at once because its later sections use 
mathematical concepts that you may not have yet met. As you work 
through the main study texts in Books A-D, you will come across sections 
or activities marked with the icon shown here. Some of these ask you to 
study part of this Guide. Other activities ask you to solve a problem using 
Maxima commands that you have previously learned. While the 
mathematical answers to such activities are included within the 
appropriate book, the computer methods needed to solve them are given in 
the section 'Computer methods for CAS activities in Books A-D' towards 
the end of this Guide. 


Maxima can be used on many different types of computer through a 
variety of interfaces. This Guide concentrates on using Maxima on a 
Microsoft Windows personal computer through the wxMaxima interface. 
This is the interface that you are recommended to use. All the Maxima 
commands mentioned in this Guide can be used with any interface to 
Maxima, unless indicated otherwise. 


Section 1 gives an overview of Maxima and includes activities in which you 
can install, test and configure Maxima on your computer. If you already 
have a working copy of Maxima ready to use, you can skip this section. 


Section 2 describes how to use the Maxima user interface for basic 
calculations and Section 3 gives an overview of using Maxima for 
manipulating algebraic expressions, plotting graphs and calculus. These 
sections relate to Unit 1 of MST125 and revise skills introduced in 
MST124. If it is some time since you studied MST124, or you have never 
studied it, you may want to work through these sections in detail, 
otherwise you should use them to refresh your memory as necessary. If you 
require further details on any of the topics covered, you might like to read 
the appropriate section of the MST124 Computer algebra guide available 
on the OU Maxima website at learni.open.ac.uk/site/maxima. If you 
are studying MST124 and MST125 starting at the same time, you won't 
need to study these sections, as MST124 will introduce the concepts as you 
need them. 


The final section of this Guide contains a summary of all the Maxima 
commands introduced. This should be useful for reference, both as you 
work through the module and afterwards. 


Note that this Guide does not attempt to describe all the features of 
Maxima. If you need more information, please consult the extensive 
documentation provided on the Maxima website at 
maxima.sourceforge.net or use Maxima's help system, which is 
described in Subsection 2.7 of this Guide. 


Introduction 


Introduction 


The *Maxima accessibility guide' (see page 95) is primarily aimed at those 
who may have difficulty using Maxima because of a disability, and provides 
some advice and suggestions. 


Maxima's origins lie in a system called Macsyma, a well-respected 
computer algebra system first developed in the late 1960s at the 
Massachusetts Institute of Technology and funded by the United 
States Department of Energy. From 1982 until his sudden death in 
2001, William (Bill) Schelter, a professor of mathematics at The 
University of Texas at Austin, maintained a version of Macsyma 
which we now know as Maxima. In 1998 he obtained permission from 
the Department of Energy for the system to be made freely available 
to all. Maxima is now dedicated to his memory. 


A word of warning 


Using a computer algebra system is not a substitute for learning the 
mathematical techniques taught in this module! You need to become 
proficient at using the methods taught so that you can apply them in the 
examination and in your future studies, and also to help you develop your 
mathematical understanding and intuition. 


Maxima is not intelligent — it is just a tool, similar to your calculator. You 
need to use your own mathematical skills to guide it and interpret its 
results critically. 


Conventions 


Here are the various computing terms used in this Guide, and how they 
should be interpreted on Microsoft Windows and Apple Mac computers. 


Computer terminology 


Term Microsoft Windows Apple Mac 
Click Click the left mouse/touchpad button Click the mouse/touchpad button 
Right click | Click the right mouse/touchpad button Hold the key while clicking 
the mouse/touchpad button 
Select Either click to select an option, Either click to select an option, 
or select text by dragging the cursor over or select text by dragging the cursor over 


the text while pressing the left mouse button | the text while pressing the mouse button 


Ctrl-Q Hold down while pressing (9) Hold down while pressing (9) 


Etc. Etc. Etc. — but there are exceptions to this 
pattern that we'll mention as they arise. 


William Schelter (1947-2001) 
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On some computer keyboards the key might be labelled ; 
or simply , and the key might be labelled : 


In this Guide, these keys are referred to as [ Enter ) and Shift), 


respectively. 


Sometimes (Ere and are two separate keys. For such keyboards 
references in this Guide to the key should be interpreted as 
references to the key. 


1 Installing, testing and configuring 
Maxima 


If you already have a working version of Maxima installed on your 
computer, for example, from your study of MST124, then you can continue 
to use that version and can skip this section. 


1.1 Introducing and installing Maxima 


Maxima can be used on many different types of computer, and there are 
several different intermediary systems, known as interfaces or 
front-ends, that provide a means of communication between you, the 
user, and the Maxima system itself. This is illustrated in Figure 1. The 
interfaces enable you to enter commands that are then processed by 
Maxima, and they display the results obtained. 


Interfaces 


wxMaxima 


XMaxima 


Command line 


Figure 1 Different interfaces for Maxima 


1 Installing, testing and configuring Maxima 


'The recommended interface is wxMaxima, which is shown in Figure 2. 
(The letters ‘wx’ in the name of this interface refer to the software used to 
create it.) 


S wxMaxima 13.04.2 [ unsaved* ] cei x 


File Edit Cell Maxima Equations Calculus 


Algebra 


Simplify Plot 


Numeric Help 


¡BS 


¿IXDOlQIQ/ÍD Oo 


e 


($11) 
($01) 


expand ( (2*x+1) * (5-3*x"2)); 
-6 x?-3 x*+10 x*5 


Ready for user input 


Figure 2 The wxMaxima interface 


If you use a screenreader, then the basic command-line interface shown 
in Figure 3 may be the easiest interface for you to use. 


E Command line Maxima 


Maxima 5.30.0 http://maxima.sourceforge.net 
using Lisp GNU Common Lisp (GCL) GCL 2.6.8 (a.k.a. GCL) 
Distributed under the GNU Public License. See the file COPYING. 


[a] > 


Dedicated to the memory of William Schelter. 
The function bug report() provides bug reporting information. 
(%i1) 


Figure 3 The command-line interface to Maxima 


There are more details on choosing an interface, configuring the 
wxMaxima interface to meet your needs and using the command-line 
interface in the ‘Maxima accessibility guide’ (see page 95). If you think 
you need this advice, you should read the start of that section now. 


Before you can use Maxima, you need to install it on your computer. You 
will do this in the next activity. 


1 Installing, testing and configuring Maxima 


Computer activity 1  /nstalling Maxima 


(a) Go to the OU Maxima website at learni.open.ac.uk/site/maxima 
and follow the instructions there to download and install Maxima on 
your computer. (You can also find a link to the OU Maxima website 
on the MST125 website.) 


(b) Create a folder on your computer in which to store your Maxima work 
for this module. You might like to put this on your computer desktop, 
or in your usual documents area. 


If you are studying within a prison or other closed institution, you can 
install Maxima from the Offline resources disc. 


1.2 Starting, testing and configuring wxMaxima 


If you are using the command-line interface you should not study this 
material, but study the alternative material in the ‘Maxima accessibility 
guide’ (see page 95) instead. If you are using a different interface, you 
should make sure that you can achieve similar outcomes with that interface. 


Computer activity 2 Starting wxMaxima 


Start wxMaxima on your computer, and keep it open to work with as you 
study this section. 


If you cannot remember how to do this on your computer, see the 
OU Maxima website: learni.open.ac.uk/site/maxima. 


The wxMaxima interface to Maxima, illustrated in Figure 4, consists 
mainly of a large white working area (called the worksheet) where you 
can type Maxima commands. These commands either instruct Maxima to 
perform calculations or change system settings. 
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G) wxMaxima 13.04.2 [ unsaved" ] 


File Edit Cell Maxima Equations Algebra Calculus Simplify Plot Numeric 


Oe®eoxkKiXDolgaieiboa 


o start using wxMaxima right away, start 
ing your command. An input cell should 
appear. Then press Shift-Enter to 
evaluate your command. 


— 


Status line 


Figure 4 The initial wxMaxima window 


When the software first starts, a secondary window showing the ‘Tip of the 
Day’ is also displayed. You can close it by clicking the ‘Close’ button. If 
you want to prevent this window from appearing the next time you start 
wxMaxima, click the box in its bottom left-hand corner, so that it is no 
longer ticked, before clicking ‘Close’. (You can always reinstate these 
messages at any time by selecting Show tips from the Help menu.) 


Along the top of the main window (or along the top of the screen in the 
case of an Apple Mac) is a list of menu names, including the File menu 
which allows you to save or print your work and re-open previously saved 
sessions (there is more on this later). You can also close the software from 
the File menu by selecting Exit. (On an Apple Mac you can do this by 
selecting Quit wxMaxima from the wxMaxima menu.) 


Below the row of menu names is a bar called the toolbar, containing a set 
of icons that provide shortcuts to commonly-used functions. 
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1 Installing, testing and configuring Maxima 


The following activity shows you how to use wxMaxima to perform a 
simple calculation and checks that your system is working properly. Follow 
the steps using wxMaxima yourself as you read the activity. 


If you encounter unexpected problems when working through this activity, 
check the Frequently asked questions (FAQ) section on the OU Maxima 
website: learni.open.ac.uk/site/maxima. 


Computer activity 3 Your first Maxima calculation 


Follow these steps to use Maxima to calculate 2 4- 3 x $. 


1. Click anywhere on the large white area of the wxMaxima window, to 
ensure that the text you type next is directed to the correct place. 


2. Type the following exactly as it appears here: 
2+3*4/5; 
including the semicolon (;) at the end. 


AAAA, 
Notice that as you type your command, it is shown in 
wxMaxima preceded by a red arrow: -->. 


This is the input prompt and shows where commands can be typed. 


'To the left of the input prompt is the cell marker 


[ 


which will be described in more detail later. 


3. Press Shift-Enter, that is, hold down the key while pressing 


CE). 


(Remember: if, as on some Macs, your keyboard has separate 


and keys, then you should hold down the key while 
pressing .) 
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If you make a mistake, try clicking elsewhere on the workspace, below 
your error, and try again. 


Pressing Shift-Enter instructs wxMaxima to send the expression to 
Maxima to calculate. This is known as evaluating the expression. 


The input prompt (-->) is now replaced by (Aii). This is how 
Maxima labels lines. The i stands for ‘input’: this is input line 
number 1. The percentage symbol (4) is used by Maxima to indicate 
objects built into the system. You will learn more about this later. 


The result is displayed under the input line, preceded by the label 
(401). The o stands for ‘output’, so this label identifies the line as 
output line 1. It gives the output corresponding to input line 1. 


G wxMaxima 13.042 [ unsaved* ] 


File Edit Cell Maxima Equations Algebra Calculus Simplify Plot Numeric Help 
DEAS XDR |e 
(%i1) 2+3*4/5; 
($01) = 
me 5 


If, in Computer activity 3, Maxima behaved as described, then you have a 
working system to use throughout MST125. Before finishing though, you 
should make a couple of small changes to how the wxMaxima interface 
behaves. 


In Computer activity 3, you used Shift-Enter rather than to 
instruct Maxima to evaluate expressions, and the input prompt (-->) 


appeared only once you started typing. Both these behaviours can be 
disconcerting, but can be changed, as described in Computer activity 4. It 
is recommended that you configure wxMaxima to work in this modified 
way, and in the remainder of this Guide it is assumed that you have done 
so. 
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Computer activity 4 Configuring wxMaxima 


1. From the Edit menu, select Configure (or on an Apple Mac 
computer, from the wxMaxima menu, select Preferences). 


Alternatively, click the following icon on the toolbar. 


This opens the ‘wxMaxima configuration’ window. 


2. Tick the ‘Enter evaluates cells’ box in the Configuration window by 
clicking it. 


This changes the behaviour so that pressing alone evaluates 


commands. 


3. Tick the ‘Open a cell when Maxima expects input’ box in the 
Configuration window by clicking it. 


This ensures that the input prompt is shown whenever 


Maxima is waiting for input. 


4. Click ‘OK’ (or close the window if using an Apple Mac). 


Closing wxMaxima 


When you have finished using Maxima, you can close Maxima and the 
wxMaxima interface by using one of the following methods. 


e Select Exit from the File menu (or Quit wxMaxima from the 
wxMaxima menu on an Apple Mac computer). 


e Press Ctrl-Q. (Here, Q stands for quit.) 


e Ona Windows computer, click the usual small cross button at the top 
right-hand corner of the wxMaxima window. 


You may be asked if you want to save your changes to the worksheet 
before closing. At this stage, click ‘No’. You will learn how to save your 
work in Subsection 2.6. 
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2 Getting started with Maxima and 
wxMaxima 


In this section you will revise how to use the wxMaxima interface to 
perform calculations with Maxima. If you are familiar with Maxima from 
MST124, you should be able to quickly skim through this section, and 
Section 3, and then try some of the activities in Section 3 to check your 
skills, revising them more thoroughly as necessary. 


If you require further details on any of the topics covered, you might like 
to read the appropriate section of the MST124 Computer algebra guide 
available on the OU Maxima website at 
learni.open.ac.uk/site/maxima. 


If you encounter unexpected problems when working through this section, 
remember to check the Frequently asked questions (FAQ) section on the 
OU Maxima website. 


If you are using the command-line interface there is an alternative version 
of the material in this section entitled ‘Getting started with Maxima using 
the command-line interface’ in the ‘Maxima accessibility guide’ (see 

page 95) which you should read instead of this section. If you are using a 
different interface, you should make sure that you can achieve similar 
outcomes with that interface. 


2.1 Calculating with Maxima 


The way in which a calculation or a command has to be entered in 
Maxima is known as its syntax. The Maxima syntax for the basic 
mathematical operations and other functions is summarised as follows. 


You can practice using some of these commands in Computer activity 5. 
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2 Getting started with Maxima and wxMaxima 


Mathematical operations and functions 


Operation or function Syntax Example 
Addition t 443; 
Subtraction - 4-3; 
Multiplication * 4*3; 
Division y 4/3; 
Brackets C and ) 2* (3+4); 
Powers, for example 2? ^ 209: 
OT ** 2**3; 


Square root, for example V5 sqrt( ) sqrt(5); 


Exponential, for example e? | %e” he’3; 

or expC ) | exp(3); 
Natural logarithm, ln log (Y log (8) ; 
Magnitude (absolute value) | abs( ) abs(-3); 
sin sin(| ) sin(1); 
cos cos( ) cos (3*%pi/2) ; 
tan tan( ) tan(Api/4); 
cosec csc (gm csc(2); 
sec sec( ) sec (%pi/3) ; 
cot cot (|) cot (4pi/4) ; 
sin”! asin( ) asin(sqrt(3)/2); 
cos”! acos( ) acos(1/sqrt(2)); 
tan”? atan( )  |atan(1/2); 


In this Guide, we will use ^ for powers. Many Maxima commands take the 
form of a command name, such as sqrt or abs, followed by a pair of round 
brackets containing one or more objects, such as numbers, that the 
command operates on. Such an object is called an argument of the 
command. In the summary tables in this Guide, arguments of commands 
are often represented by F, as illustrated in the previous table. 


Like many computer systems, Maxima assumes that all angles are 
measured in radians. The names of the inverse trigonometric functions 
sin”, cos”! and tan”! in Maxima are asin, acos and atan respectively. 
This notation is used by many computer systems; it is short for arcsin, 
arccos and arctan, the alternative names often used for sin^!, cos”! and 


tan” 


The symbols 4e and 4pi in the table above are the Maxima representations 
of the mathematical constants e and 7. The % symbol indicates that this is 
the name of a special quantity built into Maxima. When displaying 
mathematics in an output line, wxMaxima displays %pi as 7. 
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2 Getting started with Maxima and wxMaxima 


The constants e and 7 


€ ^e 4e^2; 
T Api 2*/pi; 


You should end all Maxima commands with either a semicolon (;) or a 
dollar sign ($). Ending a command with a dollar sign tells Maxima to 
perform the instruction, but not to display the answer. When using 
wxMaxima, if you forget to finish your command with a semicolon or 
dollar sign, then the system will add a semicolon for you. 


Maxima generally ignores all spaces that you type while entering a 
calculation or command, so you can use spaces to make a command easier 
to read. You can also enter line breaks, using Shift-Enter, which are 
similarly ignored by Maxima. 


Warning: Do not omit multiplication signs! 


A common error when first using Maxima is to omit multiplication 
signs. For example, you might type 


e  2sqrt(2) rather than 2*sqrt (2), or 
e 2(3+4) rather than 2* (344). 


Doing this will result in an error. The displayed error message might 
include one of the following statements: 


e  parser: incomplete number; missing exponent? 


e incorrect syntax: Syntax error 
incorrect syntax: Too many )'s 
incorrect syntax: Premature termination of input at ;. 


e  ... is not an infix operator 


If you receive one of these error messages, first check for missing 
multiplication signs! 


Warning: Take care with question marks (7?) 


Do not type question marks when asking Maxima to evaluate an 
expression. The symbol ? is a special symbol in Maxima which 
provides access to underlying systems. 


The one exception to this is when you are using the commands for 
obtaining help from the system, which you will see in Subsection 2.7. 


You can practice using Maxima for mathematical calculations in the 
following activity. 
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2 Getting started with Maxima and wxMaxima 


Computer activity 5 Calculating with Maxima 


Calculate the following using Maxima. 


(a) V32 


Since 7/32 = 325 this can be found by entering 
39^(175)5 


(The brackets in this expression are required since Maxima calculates 
expressions using the BIDMAS rules for the correct order of 
evaluation.) 


Remember to press to evaluate the command. 


After you press , the status message in the bottom right-hand 
corner of the wxMaxima window first changes to ‘Maxima is 
calculating’, then ‘Reading Maxima output’, before the result of the 
calculation is displayed. The message then displays ‘Ready for user 
input! when the system is ready to accept another command. 


Note that the first calculation of any session may take a little time. 


Remember that to enter the constant 7 in Maxima, you type Api. 


Maxima returns an exact answer, in the form of a surd. 


(c) tan 302 


Remember to convert the angle to radians, by multiplying 30? by iy 


Remember, the Maxima symbol for e is ^e. 


Since there is no simpler exact way of expressing e^, Maxima leaves 
this expression unevaluated. 
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(e) 1n3.4 


Remember that the Maxima command for a natural logarithm is 
log( ). Maxima has no command 1n( ). 


Here, since the argument is a decimal number, Maxima returns a 
decimal approximation to the answer. 


During the above activity you may have noticed that when you enter an 
opening bracket, the wxMaxima interface automatically adds a closing one 
for you. This behaviour can be turned off by unticking Match 
parenthesis in text controls in the Configuration window. (See 
Computer activity 4 for how to open the Configuration window.) 


Corresponding matching input and output lines in wxMaxima form a cell 
identified with the marker shown in Figure 5. 


Figure 5 T'he wxMaxima cell marker 


The top part of the cell marker (between the upper two horizontal marks) 
corresponds to the input line, and the bottom part corresponds to the 
output line. The cell marker turns red while an input line is being edited, 
and is shown surrounded by a box while the cell is being evaluated, as 


illustrated in Figure 6. 


Figure 6 The cell marker during evaluation 


Cells can be compressed (and uncompressed) to help make a complicated 
worksheet easier to read by clicking on the small triangle at the top of the 
cell marker. A cell can be deleted from a worksheet by clicking part of the 
cell marker other than the top triangle to highlight it, then press the 


keyboard key. 


In Computer activity 5 you saw that Maxima usually tries to calculate 
quantities exactly rather than approximating them by decimal numbers. It 
does this by manipulating the symbols in the expression according to the 
rules of mathematics. This is sometimes known as symbolic 
computation. If an exact answer cannot be obtained, the expression will 
be returned unevaluated, as in Activity 5(d). The exception to this is if 
the expression to calculate contains a decimal number, such as 

Activity 5(e), when Maxima will work in terms of decimal numbers. 


18 
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If you wish to have a result displayed as a decimal number, you can 
instruct Maxima to do so with the float command. 


Converting to decimal numbers 


Convert to a decimal number | float (F) | float(sqrt(2)); 


The name of the command float arises from the fact that the method 
used by computers to store decimal numbers internally using the binary 
digits 0 and 1 is called a floating-point representation. The float 
command instructs Maxima to convert a number stored symbolically as a 
mathematical expression to one stored using a floating-point 
representation. 


Maxima usually displays decimal numbers to 16 significant figures (though 
it suppresses trailing zeros at the end of the decimal part of a number; for 
example 0.125 is displayed simply as 0.125 rather than as 
0.1250000000000000). This is also the precision to which Maxima 
performs decimal calculations. You will see in Subsection 2.4 of this Guide 
how to change the number of significant figures displayed. 


Computer activity 6 Simplifying and approximating numerical 
expressions 


For each of the following expressions, use Maxima to simplify it and then 
find a decimal approximation for it. 


(a) V325 (b) 5200 


The expression in part (b) doesn’t involve decimal numbers and so 
the result is evaluated exactly. Notice, however, that the middle 80 
digits of the answer are not displayed because the system does not 
expect you to read them all! 


The decimal form of the result, 6.2230152778611415 101% is given 
in wxMaxima’s way of displaying scientific notation. It means 
6.2230152778611415 x 101%, 
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Troubleshooting Maxima 


If you are using Maxima and it does not seem to be responding, try the 
following suggestions. 


Troubleshooting Maxima 


Cell marker | Status message 


Boxed: 


Boxed: 


[ Maxima is calculating 


[ Ready for user input 


Maxima process terminated 


Comments 


Maxima is evaluating a command that is 
taking a long time to complete. 


Maxima is waiting for you to type 
something before continuing. 


In either case you might wish to abort the 
command in one of the following ways. 


X 
e Click Q on the toolbar. 


e Select Interrupt from the Maxima menu. 
e Type Ctrl-G (on an Apple Mac hold 

down the key while pressing CJ Vs 
Maxima might take some time to respond to 
this request. 


Maxima (but not the wxMaxima interface) 
has stopped working. 

Select Restart Maxima from the Maxima 
menu. 


If you encounter other problems, check the Frequently asked questions on 
the OU Maxima website: learni.open.ac.uk/site/maxima. 


2.2 Reusing and editing commands 


In Computer activity 6, you may have entered each expression twice, once 
to evaluate it exactly and once for a decimal answer. This is, however, not 
necessary. When you have entered and evaluated an expression, you can 
edit it and then evaluate it again, as shown in the following activity. 


Computer activity 7 Editing a command 


(a) Calculate 23 + 5/277 by entering the following line into Maxima. 


2*sqrt (3)*5*sqrt (27) ; 


(b) Edit your entry in part (a) to calculate 2/3 — 5/27 as follows. 
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D 


Click on the expression you entered in part (a), or move to it 
using the up and down keyboard arrow keys. 


2 Getting started with Maxima and wxMaxima 


As you do this the cell marker will turn red and you will see a small 
flashing vertical line (the editing cursor) appear in the expression. 


($11) 2*sqrt(3) +5*sqrt (27)5 
$01) 1743 


2. Edit the expression by typing new text and using the 
(backspace) or key to make deletions. 


3. Press ; 


Pressing re-evaluates the cell. When a cell is re-evaluated in 


this way, the input and output line numbers are updated as if a new 
cell had been added. 


Should you click in the space between two cells, a horizontal cursor will be 
shown, as in Figure 7. This indicates the position at which a new cell will 
be created, if you begin typing here. This is also shown as you move 
between cells using the up and down keyboard arrow keys. 


(311) Spi + %pi/3 + $e^2; 


4 
($01) +0? 


($12) float (%pi + $pi/3 + $e^2); 
($02) 11.57784630371704 


Figure 7 The wxMaxima horizontal cursor 


In wxMaxima, you can also copy and paste expressions, or parts of 
expressions, as you might copy and paste when using other software; that 
is, using Copy and Paste from the Edit menu, or the Ctr1-C and Ctrl-V 
keyboard shortcuts. 


Maxima also permits the result of one calculation to be used in another 
calculation. The symbol 4 represents the result of the last cell evaluated. 
Due to the updating of line numbers when a cell is re-evaluated, this might 
not be the last cell in the worksheet, but it will be the cell with the highest 
line numbers. So, for example, if you enter sqrt (4) ; Maxima calculates 
the square root of the last evaluated output. You can also use the output 
(or input) of other cells, by typing the cell line label in a calculation. For 
example, if you enter 4o05^3; Maxima calculates the cube of the expression 
in output line 5. 
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Computer activity 8 Using and editing previous calculations with 
wxMaxima 


(a) Enter cos(5*%pi/6) ; and press to evaluate it. 
(b) On the following line type 2*%; and press ; 


Maxima returns — v3, which is twice the previous answer. 


(c) Edit the line where you typed cos (5*/4pi/6) ;, change it to 
sin(5*%pi/6); and re-evaluate the cell by pressing ; 


Notice that the cell number and output of the line you edited is 
updated, but the result of 2*4; does not change. This is because 


when the 2*4; cell was evaluated, the symbol % represented the result 
of the old version of the edited command. 


(d) From the Cell menu, select the Evaluate All Visible Cells 
option, or press Ctr1-R. 


This re-evaluates all visible cells in the worksheet, in the order in 
which they appear, updating all the line numbers as it does. 


Notice that the result of 2*4; has now changed. At the time it was 
re-evaluated the symbol % represented the value of the last evaluation, 
which was the revised version of the previous line. 


(e) In a new cell at the bottom of your worksheet enter %o ^2; 
with | replaced by the output line number corresponding to the result 
of sin(5*%pi/6) ; entered earlier. 


This command finds the square of the result of the line referenced. 


Re-evaluating a worksheet 


Selecting the Evaluate All Visible Cells option from the Cell 
menu re-evaluates all visible cells in the worksheet, in the order in 
which they appear. 


This can also be achieved using Ctr1-R. 
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2.3 Assigning values to variables 


Maxima allows you to assign a value to a variable, and then use the 
variable in further calculations. 


For example, you can assign the value 23 to the variable a, and then 
evaluate an expression in a, such as a? — 3a 4- 2. You can also assign 
expressions, both numerical and (as you will see later) algebraic, and 
indeed any other Maxima object, to a variable, and use that variable in 
subsequent commands. 


To assign a value to a variable you use a colon (:). For example, the 
command 


a:23; 


assigns the value 23 to the variable a. The colon can be read as ‘is 
assigned the value’. 


Warning 


You cannot use = to assign a value to a variable. This symbol is used 
for equations, as you will see later. 


Variable names can be any combination of letters and numbers that begins 
with a letter. For example, Maxima will accept any of the following 
variable names. 


a A solution solution2 x2b 


Variable names are case-sensitive; for example, the variable x is different 
from the variable X. If a variable name is the name of a Greek letter, then 
wxMaxima will display it in output lines as the proper Greek character; for 
example, the variable alpha will be displayed as a. 


Once a variable has been assigned a value, Maxima will remember this 
value for the rest of your session. This may sometimes surprise you. It is 
easy to forget that earlier in the day you assigned a value to, say, x; you 
may then be confused later on when a calculation involving x gives an 
unexpected result. You can, however, tell Maxima to ‘forget’ about a 
variable to which you have previously assigned a value. You can also ask it 
to list all the variables that you have assigned values to. These commands 
are demonstrated in the following activity. 
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Computer activity 9 Working with variables 


(a) Assign the value of v8 to the variable a by entering 
a: sqrt(8); 


The value of a is displayed as output, in an equivalent form. 


(b) Enter a; 


This displays the current value of a. 


(c) Assign the value of ay2 to b. 


Don't forget to include a multiplication sign when inputting a2. 


The value of b is simplified when displayed. 


(d) Edit the line where you assigned v8 to a, so that the value /7 
becomes assigned to a. 


(e) Enter b; to display the value of b. 


Notice that the value of b has not changed. The variable b was 
defined when a had its original value. 


(f) Enter values; to list the names of all the variables that are currently 
assigned values. 


'The variable names are given within square brackets, separated by 


commas. 'This is how Maxima displays lists. 


(g) Enter kill(a); to remove the variable a from the system. 


You should obtain the output done. 


(h) Enter a; to display the value of a. 


'The variable name is output, as a no longer has a value. 
4 


2 Getting started with Maxima and wxMaxima 


(i) Enter b; to display the value of b. 


b is still defined. 


Notice that in Computer activity 9(d) and (e), when the value of a was 
changed, the value of b was not affected, since it was defined when a had 
its original value. If you want to update the value of b using a different 
value of a, then you should enter a new value of a then re-evaluate the line 
defining b (by selecting the line and pressing ), or re-evaluate the 
whole worksheet by selecting the option Evaluate All Visible Cells 
from the Cell menu. 


Working with variables 


Operation Command Example 
Assign a value to a variable : a:23; 
Display the value of a variable | variable a; 

List all user assigned variables | values values; 
Remove an assigned variable kill(variable) | kill(a); 
Remove all assigned variables | kill(all) kill(all); 


Note: here the placeholder variable represents any variable name. 


You can practise using variables in the following activity. 


Computer activity 10 Using variables 


(a) Define the variable a to have the value 42. 


(b) Define the variable b to be equal to sin? a + cos? a. 


2 


'To enter, for example, sin^ a in Maxima, you have to type sin(a)^2 


(c) Find the decimal approximation of the value of b. 
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2.4 System variables 


There are some variables built into Maxima whose values affect the 
behaviour of the system. These are called system variables. You met 
one system variable in Computer activity 9; the variable values holds a 
list of the names of all the variables you have defined. 


Another system variable, fpprintprec, specifies the number of significant 
figures of decimal numbers that are displayed. The name fpprintprec is 
an abbreviation of ‘floating point print precision’. 


To change the system behaviour so that, for example, only 4 significant 
figures are displayed, use the command fpprintprec:4; to assign the 
value 4 to the variable fpprintprec. 


You can set the variable fpprintprec to any value between 2 and 16. 
Also, setting it to 0 restores the default behaviour of displaying 16 
significant figures. The value of fpprintprec sets the requested number of 
significant figures, but Maxima does not always exactly meet the request. 


Computer activity 11 Displaying a specified number of significant 
figures 


Display the decimal approximation of 7 to 8 significant figures by doing 
the following. 


(a) Set the system variable fpprintprec to be 8. 
(b) Display the decimal approximation of 7. 


You can reset the values of fpprintprec and all other system variables to 
their original values by using the reset O command. Try entering 
reset(); followed by float (Api); and check that the value of 7 is 
displayed to 16 significant figures. Note that the reset command also 
resets your Maxima line numbers to start at 1 again. 


Resetting system variables 


Reset all system variables 


Note that the list of variables displayed when you enter values; does not 
include any system variables. If one of your variables does not appear in 
the list displayed by values;, then you have probably used a variable 
name that is also the name of a system variable and hence changed the 
value of that system variable. In extreme circumstances this may change 
the behaviour of Maxima, which can be restored by resetting all system 
variables as described above. 
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2.5 Annotating your wxMaxima worksheet 


It is important to be able to clearly explain and present your mathematical 
work, both so that someone else can understand it, and so you can 
understand what you have done if you return to your work after some 
time. To help you do this, wxMaxima allows you to enter text comments 
within your worksheet, as illustrated in Figure 8. 


@ wxMaxima 13,04.2 [ unsaved" ] 
File Edit Cell Maxima Equations Algebra Calculus Simplify Plot Numeric Help 


OU wm esxix0adiiale!»o: |e 


Circles 


(441) r:3; 
($01) 3 


($12) c: float(2*$pi*r); 
($02) 18.84955592153876 


($i4) a: float($pi*r^2); 
($04) 28.27433388230814 


E --> 


Figure 8 Adding text to a wxMaxima worksheet 


To enter text into your worksheet, first you need to create a suitable cell in 
which to type. To do this first you position the cursor, then you select the 
appropriate one of the following commands from the Cell menu. 


e Insert Text cell 

e Insert Title cell 

e Insert Section cell 

e Insert Subsection cell 


'The new cell will be created under the cell in which the editing cursor 
appears, or between two existing cells at the position indicated by the 
horizontal cursor, if shown. Remember, the horizontal cursor appears 
whenever you click between two cells, or move within the worksheet using 
the up and down keyboard arrow keys. 


As you type your text, to start a new line within the cell press a 
When you have finished typing the text, click elsewhere on the worksheet 
or use the up or down keyboard arrow keys to move out of the text cell. If 
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your text cell is at the bottom of the worksheet, and you click or move 
below it, an input prompt may not be shown until you start typing again. 


Section and Subsection cells are automatically numbered. 


Computer activity 12  /ncluding text in a worksheet 


(a) Create the worksheet shown in Figure 8. Use a Title cell for the title 
‘Circles’ and Text cells for the remaining text. Remember to click 
elsewhere on the worksheet or move position with the up or down 
keyboard arrow keys when you have finished entering each piece of 
text. 


After evaluating a Maxima command, you will need to correctly 
position the horizontal cursor before inserting the next text cell. 


(b) Use your worksheet to find the circumference and area of a circle of 
radius 5. (Remember to re-evaluate the worksheet, for example, by 
using Ctrl-R, after changing the value of r.) 


Each Title cell, Section cell and Subsection cell has a small square in its 
top left corner. Clicking this square hides (or reveals) the contents of the 
worksheet from that point onward, until the next cell of the same type. 


2.6 Saving and printing your work 


After working with Maxima you will probably want to save your 
calculations so that you can reuse or read them at a later date. You can do 
this in wxMaxima using the Save or Save As options from the File menu. 


There are three different formats in which you can save your work. These 
are as follows. 


e wxMaxima document, with the .wxm file type. This saves all your 
input but not any calculated output. You have to recalculate the 
output when you re-open the file. 


e wxMaxima xml document, with the .wxmx file type. This saves all 
the input and the output. 


Note that when you load a worksheet previously saved as a wrMazima 
zml document (*.wxmz), although both the saved input and output is 
displayed, the input has not been evaluated by Maxima, so you cannot 
immediately use any variables defined in the worksheet. To re-evaluate 
the worksheet, and re-define all the variables used, select Evaluate 
All Visible Cells from the Cell menu, or press Ctr1-R. 


e Maxima batch file, with the .mac file type. This is for more 
experienced users of Maxima. 
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It is recommended that you save your wxMaxima work as a 
wxMarima xml document, with the .wxmx file type. 


You can open a saved wxMaxima worksheet by using the Open option from 
the File menu, or by double-clicking on the saved file on your computer. 


Computer activity 13 Saving your work 


Save your current worksheet then reload it, by following these steps. 
(a) Select Save As from the File menu. 
(b) In the window that appears: 


(i) Choose the folder in which to save your work. You might like to 
save it in the folder that you created in Computer activity 1. 


(ii) Choose a name for the file in which the worksheet will be saved. 


(iii) Make sure the file type is set to wxMaxima xml document 
(iv) Click ‘Save’. 
(c) Close wxMaxima, by selecting Exit from the File menu, or typing 


Ctrl-Q, or clicking the small ‘x’ at the top right of the window. (On a 
Mac select Quit wxMaxima from the wxMaxima menu or hold down the 


key while pressing (Q) .) 


(d) Open your worksheet in wxMaxima again, either by 
(i) double clicking on the file you just saved, or 
(ii) following the steps below: 
e Start wxMaxima again. 


e Select Open from the File menu. Maxima will ask if you 
want to save your just opened, blank worksheet. Click ‘No’. 


e Find the file containing your work and click ‘Open’. 


(e) Re-evaluate the worksheet by selecting Evaluate All Visible Cells 
from the Cell menu or pressing Ctrl-R. Then type values; to list the 
known variables. 


All the variables used in the worksheet should now be known to 


Maxima and can be used in new calculations. 


You can print a wxMaxima worksheet by selecting Print from the File 
menu. 
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Alternatively, you can obtain an image of the entire wxMaxima window by 
taking a ‘screenshot’. To do this on a Microsoft Windows computer, first 
click on the wxMaxima window, then press Alt-PrintScreen (that is, 
hold down the keyboard key while pressing the PrintScreen key, 
which may be labelled { PrtScr ) , {Prt Scrn Jor something similar). This stores 
the image, which you can then paste into a suitable application. If you are 
using Windows Vista, Windows 7 or Windows 8, you might like to use the 
‘Snipping Tool’, available by typing ‘snip’ into the Start menu. To take a 
screenshot on an Apple Mac computer, press Cnd-Ctrl-Shift-4 (that is, 
hold down the ; and keyboard keys while pressing the (4) 
key). Next press the space bar, move the camera pointer over the 
wxMaxima window and click. This stores the image which you can then 
paste into a suitable application. You may also like to use the Grab 
application from the Utilities folder. 


If you want to copy a single output line from your worksheet into, for 
example, a word-processed document, then you can select that part of the 
worksheet, and choose Copy As Image from the menu obtained by right 
clicking on what you have selected. You can then paste the resulting image 
into a suitable application. 


2.7 Getting help 


There are several ways in which you can obtain help with using Maxima, 
beyond this Guide. In wxMaxima, you can access the complete Maxima 
manual by selecting Maxima Help from the Help menu. This opens the 
window shown in Figure 9. You can use the Contents, Index and Search 
tabs to look for the information that you want. 


[Top] [Contents] [Index] [?] 
Maxima 5.30.0 Manual 


Maxima is a computer algebra system, implemented in Lisp. 


Maxima is derived from the Macsyma system, developed at 
MIT in the years 1968 through 1982 as part of Project MAC. 
MIT tumed over a copy of the Macsyma source code to the 
Department of Energy in 1982; that version is now known as 
DOE Macsyma. A copy of DOE Macsyma was maintained by 
Professor William F. Schelter of the University of Texas from 
1982 until his death in 2001. In 1998, Schelter obtained 
permission from the Department of Energy to release the DOE 
Macsyma source code under the GNU Public License, and in 
2000 he initiated the Maxima project at SourceForge to 
maintain and develop DOE Macsyma, now called Maxima. 


Maxima infrastructure 


1. Introduction to Maxima Sample Maxima sessions. 
2. Bug Detection and Finding and reporting bugs in 


Figure 9 The Maxima help system 
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Maxima also includes commands that you can use to obtain help. These 
are demonstrated in the following activity and summarised in the table 
that follows. 


Computer activity 14 Getting help! 


(a) Enter ? float; to show the help for the float command. 


Note the space between the ? and float, which is needed. This is one 
of the few times in Maxima when spaces matter! 


The help information for the command may include more detail than 
you need! 


After giving the help information, Maxima displays the output true. 
This means that the command you entered was successful. 


(b) Suppose that you cannot remember the Maxima command for square 
root, but you do remember that it was something like sqr. 


Type ?? sqr; to list all the help information titles that contain the 
letter sequence ‘sqr’. 


A list of possible titles is displayed. 


Title number 1 is what you were looking for, so enter 1; to display the 
relevant information. 


If more than one title is displayed in response to a ?? command, then 
Maxima will not continue until you enter one of the following: 


e the appropriate number, to display your chosen information 


e several numbers, separated by spaces, to see help information 
on several topics 


e all to display the information on all topics listed, or 
e none to see no information. 


A box will be shown around the cell marker while Maxima is waiting 
for your response. 
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Help commands 


Operation Command Example 


Get help on a command ? command ? float; 
or describe( command ) describe(float); 


Find help information whose ?? ?? flo; 
title contains the given text | or describe( , inexact) | describe(flo, inexact); 


Note 1: the space after the ? is needed. 


Note 2: the second pair of commands displays a list of the titles of all 
help information pages whose title contains the given text. To obtain 
help on a particular topic listed, type the number corresponding to the 
required title in the list, followed by a semi-colon. 


There is further help with Maxima in the Frequently asked questions 
section of the OU Maxima website, at learni.open.ac.uk/site/maxima. 


3 Using Maxima 


In Section 2 of this Guide you revised how to use wxMaxima and how to 
perform numerical calculations. Here, you will revise how to use Maxima 
for algebraic manipulations, equation solving, calculus and graph plotting. 


If you are familiar with Maxima from MST124, you should be able to 
quickly skim through this section and try some of the activities to check 
your skills, revising them more thoroughly as necessary. If you require 
further details on any of the topics covered, you might like to read the 
appropriate section of the MST124 Computer algebra guide available on 
the OU Maxima website at learn1.open.ac.uk/site/maxima. 


3.1 Algebra and functions 


You can enter algebraic expressions into Maxima in a similar way to 
numerical expressions, using the operations and functions listed in 
Subsection 2.1. You can assign algebraic expressions to variables, using the 
assign (:) command and you can also define your own functions using the 
:= operator. These are demonstrated in the following activity. 
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Computer activity 15 Defining algebraic variables and functions 


(a) Define the variable p in Maxima to be the algebraic expression 
2x + V 1-4 3x 


(b) Find the value of p when x — 3 by entering 
subst(3,x,p); 


The command name subst is short for ‘substitute’. 


The command subst (3,x,p); substitutes the value 3 for the variable 
x in the expression p. 


z*-c-r41l 


r3 — 13x +12’ 
by entering the following command. (Use the characters : followed 


by = to type the := operator.) 


(c) Define the function f in Maxima, with the rule f(x) = 


f (x) :=(x^4+x+1)/ (x^3-13*x*12) ; 


(d) Calculate the value of f (2), by entering £(2) ; 
(e) Try to calculate the value of f (3). 


'This gives an error. 


This is because when Maxima tries to evaluate f(3), it finds that this 
involves dividing by zero, since the denominator of f is zero when 

x = 3, so it returns an error message. The function f is undefined at 
r=3. 


(f) Enter g(1); 


Maxima returns the expression g(1), unevaluated. 


Since the function g has not been defined, Maxima cannot calculate a 
value. It knows only that the value is g(1). 


(g) Enter functions; 


'This displays the value of the system variable functions, which is a 
list of all your currently defined functions in this Maxima session. 


3 Using Maxima 
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Notice the difference between defining a function using := and assigning 
an expression to a variable using :. A function has an argument and can 
be evaluated at different values of the argument, while a variable does not. 


Once you have defined a function in Maxima it is remembered by the 
system for the rest of your session, unless you remove it using the kill 
command that you revised in Subsection 2.3 of this Guide. 


Substituting into an expression 


Substitute | subst (value , variable, | subst (4, x, x^241); 
expression) which substitutes 4 for x in x^2*1. 


Note: expression indicates that the argument of the command must 
be an expression, or a variable whose value is an expression. 


Working with functions 


Operation Command Example 
Define a function = f (x) :=2*x+3; 
Evaluate a function at a value | function ( ) f£); 

List all user defined functions | functions functions; 


Remove a defined function 


kill(function) | kill(f); 


Maxima can also re-express algebraic expressions in different, but 
equivalent forms. It is hard for a computer to automatically simplify 
algebraic expressions, because often there is no single ‘simplest form’ for a 
particular expression, and the computer has no way of knowing which form 
you prefer. 


Instead, a range of commands are provided (which are listed in the 
following table) which you can use to instruct Maxima to attempt to 
re-express an expression in different ways. You can use these commands to 
simplify an expression in the way you want. 
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Algebraic and trigonometric manipulation commands 


Operation Example 


Expand brackets expand( ) expand( (x+1)%2 ); 

Factorise factor( ) factor( 2*x+4*x%2 ); 

Simplify fullratsimp( ) | fullratsimp( (2*x+4*x%2)/x ); 

Simplify something involving radcan( ) radcan( log(x^2) ); 
exponentials and logarithms 

Combine logarithms logcontract(| ) | logcontract (log(a)+log(b)) ; 

Expand trigonometric functions trigexpand( ) |trigexpand(sin(A+B)); 


of sums, differences and 
multiples of angles 
Express powers of sines and trigreduce( ) | trigreduce(sin(x)^2); 
cosines in terms of sines and 
cosines of multiple angles 


Simplify trigonometric expressions | trigsimp( ) trigsimp(sin(x)^24cos(x)^2); 
Simplify algebraic fractions trigrat( ) trigrat((sin(2x))/cos(x)); 
containing trigonometric 


functions 


Computer activity 16 — Rearranging algebraic expressions 


(a) Enter 
expand( (1+sin(x))*(x+4)%2 ); 


to expand the expression (1 + sin x)(x + 4)?. 


'The command expand multiplies out the expression that is its 
argument. 


Spaces have been included simply to help with the reading of the 
command. 


(b) Enter 
factor( 3*x%3+11*x%2+8*x-4 ); 


to factorise the expression 3x? + 11z? + 8x — 4. 


The command factor factorises, where possible, the expression that 


is its argument. 


(c) Enter 
fullratsimp( (2x*x^3-3*x^241)/(x-1) ); 
2x3 — 33? +1 


to simplify the expression 1 
A be — 
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'The command fullratsimp simplifies, where possible, the expression 
that is its argument. It does this by using techniques for multiplying 
out brackets, collecting terms over a common denominator and 
cancelling out common factors. Simplifications requiring other 
techniques may not be handled by fullratsimp in the way you would 
expect. 


The name fullratsimp is short for ‘full rational simplify’. 


(d) Enter 
radcan( log(x^3-3*x-2)-2*1log(x*1) ); 


to simplify the expression In(z? — 3x — 2) — 21n(z + 1). 


The command radcan uses the rules of indices and logarithms to 
simplify, if possible, an expression containing roots (for example 
square roots or cube roots), powers and logarithms. 


Note that roots are often known as radicals. The radcan converts 
expression involving radicals (and related functions such as powers 
and logarithms) into a canonical (or standard) form. 


(e) Enter 
logcontract( log(2*x)-2*1og(x)-*1log(y) ); 
to simplify the expression In(2x) — 21n x + In y. 


The command logcontract combines, where possible, several terms 


involving logarithms. 


(f) Enter 


trigexpand( tan(3*x) ); 


'The trigexpand command attempts to express trigonometric 
functions of sums, differences and multiples of angles in terms of 
powers of trigonometric functions of individual angles. 


(g) Enter 


trigreduce( cos(x)^4 ); 


'The trigreduce command attempts to express powers of sines and 


cosines of angles in terms of sines and cosines of multiples of the 
angles. It reduces the powers of the functions. 
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(h) Enter 
trigsimp( cos(x)^2*(1*tan(x)^2) ); 


to simplify cos? x (1 + tan?(z)). 


'The trigsimp command attempts to simplify expressions involving 


trigonometric functions using the fact that sin? x -- cos? x — 1. 


(i Enter 
trigrat( tan(x)/sin(2*x) ); 
ng 


to simplify <2 
o simplify ———.. 
RM sin(2x) 


'The trigrat command attempts to simplify algebraic fractions 


involving trigonometric functions. 


You can also access many of these commands for manipulating expressions 
from the wxMaxima menus. For example, to factorise an expression you 
can first enter the expression in the worksheet, select it by highlighting the 
text using the mouse or clicking on the corresponding output line number, 
then select Factor Expression from the Simplify menu (or from the 
menu obtained by right clicking on a highlighted output line number). 


Other options on these menus, such as Substitute... open a window for 
you to enter the command options. Options that work in this way are 
identified by ... at the end of the option name. 


3.2 Plotting graphs 


You can plot graphs of expressions or functions in wxMaxima using the 
wxplot2d command. This command inserts the graph into your worksheet, 
as demonstrated in the following activity. Maxima plots the graph by 
calculating a large number of points on the curve and joining them up. 
The ‘2d’ in the name of this command comes from the fact that its graphs 
are two-dimensional plots. The ‘wx’ is because it is a wxMaxima command. 


If you are using another interface to Maxima, you should use the command 
plot2d instead. You use it in the same way as wxplot2d, but it will plot 
the graph in a different computer window, which you then need to close 
before you can create another plot. 


3 Using Maxima 
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Computer activity 17 Plotting graphs 


(a) Plot the graph of the equation y = xsin x for values of x such that 
—5 <a € 5, using the command 


wxplot2d(x*sin(x), [x,-5,5]); 


The first argument of the command is the expression to be plotted, 
and the second argument specifies the variable in the expression 
together with the least and greatest values that you want it to take in 
the graph. These are grouped in square brackets and separated by 
commas, which is the Maxima syntax for a list. You will revise more 
details of lists later. 


Note that wxMaxima displays the graph next to a line label 
beginning %t, such as (/t1). This is how wxMaxima shows 
intermediate output, that is, output other than the final result of the 
command. The wxplot2d command has no final result, just 
intermediate output, which is the graph. 


2 
+2 
(b) Define the function f with the rule f(t) = WIS 


graph of f for values of t such that 0 € t < 3, using the command 
vwxplot2dC0f(t), [t,0,3]); 


and plot the 


Notice that the first argument of the wxplot2d command must be 
f (t) rather than just f and that since the independent variable is t, a 
range of values of t is given as the second argument. 


The graph produced shows a large spike near t = 1. This is because 
f (t) is not defined at t = 1, since its denominator is zero there. 


'To plot the graph, Maxima has evaluated f at values of t in the range 
0 € t € 3, obtained large values near t = 1 and joined up the 
corresponding points. 


'To obtain a better graph, you should restrict the range of the vertical 
axis. You are asked to do this next. 


(c) Plot the graph of f over the range 0 € t € 3, with the vertical range of 
the graph restricted to —10 € y € 10 using the following command. 


wxplot2d(f(t), [t,0,3], [y,-10,10]); 
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The range of the vertical axis is restricted by including an additional 
argument which is a list starting with the keyword y and followed by 
the minimum and maximum values to be shown on the axis. 

A warning: 

plot2d: some values were clipped. 

is given. This is because some of the points on the graph with 


x-coordinates between 0 and 3 could not be displayed due to the 
range of values of y specified. 


You can also plot multiple graphs in the same diagram by including a list 
of expressions as the first argument of the command, as shown in the 
following activity. 


Computer activity 18 Plotting multiple graphs 


(a) Define the functions p and q to be 


p(x) = x cos x q(x) =e 


(b) Plot graphs of the two functions over the range —3 € x < 3 on the 


same diagram by entering: 


vxplot2d([p(x) ,q(x)], Dx, -3,3]) ; 


Notice that the first argument of the command wxplot2d is a list 
containing the two functions to be plotted. 


The curves are automatically plotted in different colours, and a legend 
showing which function corresponds to which curve is included at the 
top right-hand corner of the graph. 


Maxima lets you control many properties of the graph, for example the 
colour of the curves plotted and the text in the legend. This is done by 
including additional arguments in the wxplot2d command, each being a 
list beginning with a keyword describing the property to be changed and 
followed by the required value of the property. Where properties apply to 
individual curves, a value should be listed for each curve plotted. This is 
demonstrated in the following activity. 
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Computer activity 19 Changing the properties of a graph 


Plot the graph of the functions p and q from Computer activity 18 a 
second time, with 


e p coloured black and q green, 

e the graph legend changed to refer to the names p and q rather than 
the expressions they represent 

e the vertical axis labelled y 

by entering 
wxplot2d([pGO ,qGO ], [x,-3,3], [color, black, green], 
[legend, "p","q"], [ylabel,"y"]); 


'The additional settings are given as lists, each beginning with a 
keyword. 


The colours of the curves are given in a list beginning color 
(American spelling), and one colour is listed for each curve. Maxima 
uses the first colour listed for the graph of the first expression listed 
and so on. If you do not include enough colours, Maxima starts using 
the listed colours again. 


The list beginning legend specifies the legend text for each curve, 
each enclosed in double quote marks. 


The list beginning ylabel sets to text to be displayed on the vertical 
axis. If required, text for the horizontal axis label can be specified 
using a list beginning with the keyword xlabel. 


Warning 


When plotting graphs, do not assign the expression to be plotted (or 
any other expression) to the variable y as this will cause an error when 
the keyword y is used in the specification of the range of the y-axis. 
Similarly, the variable in the expression to be plotted should not have 
been previously assigned to a specific value. (Remember, a variable 
with an assigned value can be deleted using the kill command.) 


It is, however, possible to assign the expression to be plotted to the 
function y. 


The different ways in which you can use the wxplot2d command are 
summarised as follows, along with some of the optional arguments. If you 
are using a different interface to Maxima, you can use the plot2d 
command in the same ways. 
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Plotting graphs 
Operation 


Plot a graph 


Plot several graphs 


Command 


wxplot2d (expression , 
horizontal range,.. 


wxplot2d(list of expressions , 
horizontal range,.. 


Option 


Vertical range 
Curve colour 
Legend text 

Turn legend off 
Horizontal axis label 
Vertical axis label 


Note: for a full list of the possible options, view the Maxima help 


Graph plotting options 


Argument 


[y, RB. NI 
[color,  ,...] 
[legend," ",...] 
[legend, false 
[xlabel, " "] 
[ylabel, "| "] 


Example 


Example 


[y, 0, 5] 

[color, red] 
[legend, "Car A"] 
[legend, false] 
[xlabel, "t (h)"] 
[ylabel, "s (km)"] 


obtained using the ?? plotting options; command. 


Plotting curves represented by implicit equations 


'The wxplot2d command can only be used to plot curves represented by 
equations in ezplicit form; that is, written in the form y — f(x), where f is 


a function. 


To plot a curve represented by an equation in implicit form, such as 
z? + y? = 1, you need to use the wximplicit plot command. This is 


demonstrated in Computer activity 20. 


The wximplicit plot command is not part of the basic Maxima collection 
of commands, so you have to add it by loading an additional Maxima 


package, the implicit.plot package, which was installed on your 


computer when you installed Maxima. The Maxima command to load an 
additional package is load, which is also demonstrated in Computer 
activity 20. Each package needs to be loaded only once per Maxima 


session. 


If you are using a Maxima interface other than wxMaxima, then you 
should use the implicit. plot command instead of the wximplicit plot 


command. You load and use it in exactly the same way as 


wximplicit.plot. 


wxplot2d(x^2, [x,0,1]); 
z) which plots the graph of x? 
for x between 0 and 1. 


wxplot2d(x^2,[x,0,11, [y,0,2]) ; 
which plots the graph of x? 
for x between 0 and 1, and with 
vertical axis values between 0 and 2. 


wxplot2d([x^2, 2*x], [x,0,1]); 
>) which plots graphs of z? and 2x 
for z between 0 and 1. 


3 Using Maxima 
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Computer activity 20 Plotting a circle 
Plot the circle represented by the equation x? + y? = 1, as follows. 
(a) Load the implicit.plot package by entering 
load(implicit.plot); 


Remember to type the underscore symbol (_) in the middle of the 
package name. 


'This package provides the wximplicit.plot command. 


(b) Plot the circle by entering 
wximplicit_plot (x*2+y*2=1, [x,-2,2], [y,-2,2]); 


Again, remember to type the underscore symbol (_) in the middle of 
the command name. 


The first argument of the command is the equation to be plotted. 
Notice that in Maxima, equations are given using the usual equals 
sign. 


The second argument specifies the variable to be plotted on the 
horizontal axis, and the range of values of that variable to be plotted. 


The third argument specifies the variable (and its range) to be 
plotted on the vertical axis. You must include this third argument. 


The circle with equation £? + y? = 1 has centre (0,0) and radius 1. 
So the circle lies within —2 € x € 2 and —2 < y < 2 and hence these 
seem suitable ranges to use. 


'The curve that Maxima plots, however, does not look very circular! 
This is discussed below. 


Note that the first argument of the wximplicit plot command is an 
equation, whereas the first argument of the wxplot2d command is an 
expression (or a function or variable representing an expression). If an 
expression is given as the first argument, wximplicit.plot will plot the 
graph of the equation obtained by setting the expression equal to zero. 


The curve that you plotted in Computer activity 20 does not look circular 
because it has different scales on the axes. Whatever ranges of x- and 

y-values you specify for a plot, by default Maxima makes the y-axis appear 
3 


5 as long as the z-axis. 


To specify the use of equal scales, you can add the following additional 
argument to the wximplicit_plot command: 


[gnuplot_preamble, "set size ratio -1"] 
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Gnuplot is the underlying software that Maxima uses to plot graphs. In 
general, if you set the size ratio to be —r, where r is a positive number, 
then the length that represents 1 unit on the y-axis is r times as long as 
the length that represents 1 unit on the z-axis. On the other hand, if you 
set the size ratio to be r, where r is a positive number, then the y-axis 
appears r times as long as the z-axis. 


Computer activity 21 Plotting a circular circle! 
Plot the circle z? + y? = 1, using equal scales on the z- and y-axes by 
entering 


wximplicit.plot(x^2*y^2-1,[x,-2,2], [y,-2,2], [gnuplot.preamble,"set size ratio -1"]); 


Including the gnuplot_preamble argument in the wximplicit plot 
command can be tiresome if you are plotting lots of curves and need 
equally scaled x- and y-axes for each. An alternative is to change the 
default behaviour of Maxima so that it uses axes with equal scales for the 
rest of your session. You will see how to do this in the following activity. 


If you have configured Maxima to change the font and line styles used 
when plotting graphs, as described in the Accessibility section of the OU 
Maxima website, then the command given in the following activity will 
revert these back to their original styles. To retain your graph style 
settings, use 


set plot option([gnuplot preamble, OUEqualScales]); 


instead. 


Computer activity 22 Plotting another circular circle! 
Plot the circle (x — 1)? + y? = 3, using equal scales on the z- and y-axes as 
follows. 


(a) Tell Maxima to always use equal scales on the z- and y-axes by 
entering 


set plot.option([gnuplot.preamble, "set size ratio -1"])$ 
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This sets [gnuplot preamble, "set size ratio -1"] to bea 
default option for all plots during the rest of your Maxima session. 
The line ends with $ to prevent the display of the output of the 
command. When present this output consists of a long list of all the 
current default settings, which can safely be ignored. 


(b) Plot the circle by entering 
wximplicit.plot((x-1)^2*y^2-3, [x,-1,3], ly,-2,2]); 


If you wish to stop plotting graphs with equal scales and restore the 
default Maxima plotting behaviour, then set the gnuplot_preamble 
plotting option to be blank using 


set_plot_option([gnuplot_preamble, ""])$ 


(or use set_plot_option([gnuplot_preamble, OUAxisSettings])$ if 
you have changed your preferred font and line styles). 


To plot more than one curve at once you can list their equations within 
square brackets as the first argument of wximplicit plot, as you did with 
the wxplot2d command. The equations do not necessarily have to be in 
implicit form. This is demonstrated in the following activity. 


Computer activity 23 Plotting a parabola and a circle 


Plot the parabola y = x? — 6x + 9 and the circle (x — 3)? + (y — 2)? = 3 
over the range 0 < x < 6,0 € y € 6 by entering 


wximplicit plot([y2x^2-6*x-9, (x-3)^24(y-2)^2-3], [x,0,6], [y,0,6]); 


The commands introduced in this section are summarised as follows. 


Loading additional packages 


Load a package | Load (package name) | load(implicit plot); 
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Plotting curves represented by equations in implicit form 


Operation Command Example 

Plot a curve represented | wximplicit plot(equation, | wximplicit plot(x^2*y^2-1, 
by an equation in horizontal range, [x,-1,1], [y,-1,11); 
implicit form vertical range ,...) 

Plot several curves wximplicit plot( wximplicit.plot( 
represented by list of equations, [x^242*y^2-4, y-x^3], 
equations in implicit horizontal range, [x,-1,1], [y,-1,1])D; 
form vertical range ,...) 


Note 1: the implicit. plot package needs to be loaded before this 
command can be used. 


Note 2: the first argument is an equation or a list of equations. 
Alternatively it can be an expression, or a list of expressions, in which 
case the equation formed by setting the expression equal to zero is plotted. 


Note 3: equations in explicit form can also be plotted using this command. 


Graph plotting options 


Option Argument Example 

Change the relative lengths/scales | [gnuplot preamble, [gnuplot preamble, 
of graph axes. (A positive "set size ratio number "] "set size ratio -1"] 
number specifies the ratio of which sets the axes to have 
axis lengths, a negative number equal scales. 
specifies the ratio of axis 
scales.) 


Setting default plotting options 


Set the default plotting options | set plot option( option) | set plot option([gnuplot.preamble, 
"Set size ratio -1"]); 


3.3 Solving equations 


You can use the Maxima solve command to find the exact solution(s) of 
equations, and to change the subject of an equation, as demonstrated in 
the following activity. 


Remember, to type an equation in Maxima, you use an equals sign in the 
usual way. For example, you can type 3*x*2-4. You can also assign 
equations to variable in a similar way to expressions. For example, to 
assign the equation 3x + 2 = 4 to the variable q, enter 


q: 3*x+2=4; 


45 


3 Using Maxima 


Computer activity 24 Solving and rearranging equations 


(a) Solve the equation z? — 4x? + 2x + 4 = 0 by entering 


solve( x*%3-4*x%2+2*x+4=0 ); 
This is a cubic equation, since the highest power of x is 3. Each 
cubic equation has at most three solutions. 


Notice that Maxima displays the solutions of the equation as a list. 


(b) Assign the equation y = In(ux + v) to the variable eqn. 


Remember, the Maxima function for In is log. 


Rearrange the equation to make x the subject by entering 


solve(eqn, x); 


'The second argument of the solve command gives the variable we 


wish to ‘solve’ for. 


(c) Solve the equation tan(3x — 1) = V3. 


When solving equations involving trigonometric functions Maxima 
will only try to find one solution, even though many others may exist. 


'The warning 


solve: using arc-trig functions to get a solution. 
Some solutions will be lost. 


is given to alert you to this. 


(d) Solve the equation x? — 4x + 20 = 0. 


The solutions are given in terms of %i. 


This is how Maxima represents the imaginary number 7 whose square 
is equal to —1. There are no solutions to this equation that are real 
numbers, but Maxima has found solutions that are complex numbers. 


(e) Try to use Maxima to solve the equation z? + 2z — 1 = 0. 
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Here, Maxima returns (within a list) the equation, slightly 
rearranged, but unsolved. 


'The solve command always tries to find an exact solution of the 
equation using algebraic manipulations. If it is unable to do this, then 
it returns the original equation, as here. 


You will revise later in this section how to find approximate numerical 
solutions of such equations. 


The use of the solve command illustrated in the previous activity is 
summarised below. 


Solving and changing the subject of an equation 


Solve an equation, exactly solve( equation) solve (2*x*2-1=0) ; 
Solve an equation for a variable | solve (equation, variable) | solve(2*a*b-3*b=0, a); 


You saw in Computer activity 24 that the output of the solve command is 
a list of equations. To extract the individual solutions you need to be able 
to manipulate lists. This is demonstrated in the following activity. 


Computer activity 25 Working with lists 


(a) Assign the solutions of the equation z? — z — 1 — 0 to the variable 
solns by entering 


solns:solve(x*2-x-1=0); 


Remember, the solutions of an equation are given as a list, so solns 


is a list. 


(b) Display the first element of the list solns by entering 
solns[1]; 


The individual elements of a list can be obtained by adding the 
element number (which is known as the index of the element) after 
the name of the list and within square brackets. 


V 5-1 


In this case, the first element of the list is the equation x=-———. 


The value of the first solution of the original equation is the 
right-hand side of this. 
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(c) Display the exact value of the first solution of x? — x — 1 = 0 by 
entering the command 


rhs(solns[1]); 


The command rhs gives the right-hand side of an equation. 


(d) Display the exact value of the second solution of the equation similarly. 


Lists 


Extract an element of a list A[2] ; 


Extracting the left- and right-hand sides of an equation 


Left-hand side of an equation | lhs(equation) | Ihs (4*x41-2*x-2) ; 
Right-hand side of an equation | rhs (equation) | rhs (4*x+1=2*x-2) ; 


Another use of the solve command is to solve two or more simultaneous 
equations. To use the command in this way, you must enter the equations 
as a list, and enter the variables to be solved for as another list, as 
indicated below. 


Solving simultaneous equations 


Solve simultaneous equations | solve (Tist of equations, | solve([2*x+y=4, x-2*y-1], 
list of variables) [x,y]); 


Computer activity 26 Solving simultaneous equations 


(a) Use Maxima to solve the simultaneous equations 
3x + 2y=1 
52+ 3y = 3 


as described below. 


(i) For simplicity, assign the first equation to the variable eq1 and 
the second equation to the variable eq2. 


Maxima might reorder the terms of the equations when displaying the 


output. 
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(ii) Solve the pair of equations (eq1 and eq2) for the two unknowns 
(x and y) by entering 


solve( [eqi,eq2], [x,y] ); 


Maxima’s output is [[x=3, y=-4]] 


The solution list contains the element [x=3, y=-4] which is itself a 
list of the pair of values of x and y that solve the equations. 


(b) Solve the equations 
y=2*-—6r+9 
(a — 3)? + (y- 2}? =3 


which gives the points of intersection of the parabola and circle plotted 
in Computer activity 23. 


Notice there are four solutions, Maxima gives a list of each pair of 


values of x and y which satisfy the equations. 


In Computer activity 24(e) you saw an example where the solve command 
failed to find the exact solution of an equation. In such circumstances you 
can use another Maxima command, find.root, to try to find an 
approximate solution. 


Whereas the solve command manipulates the equation mathematically to 
find a solution, the find root command performs a series of numerical 
calculations to find a decimal approximation to a solution. This is known 
as solving the equation numerically, or finding a numerical 
approximation to a solution. 


The solutions of an equation of the form f(x) = 0 are the x-coordinates of 
the points at which the graph of y = f(x) crosses the x-axis. If f is a 
continuous function and |a, b] is an interval included in its domain such 
that one of f(a) and f(b) is negative and the other is positive, then the 
graph of f(x) must cross the x-axis at one or more points in the interval 
(a, b), as illustrated in Figure 10. So there must be at least one solution of 
the equation f(x) = 0 in the interval (a, b). 


3 Using Maxima 
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y^ y^ 


S 
yY 


Figure 10 There is at least one solution of f(x) = 0 in the interval (a,b) 


if f(a) and f(b) have different signs 


The find-root command looks for a single solution of an equation of the 
form f(x) = 0 that lies within an interval that you specify. The command 
requires that the values of f at each endpoint of the interval have different 
signs. If this is not the case, then Maxima gives an error message. 


Solving an equation numerically 
Operation Command 


Solve an equation numerically | find_root (expression , 
variable S, 

interval start value, 
interval end value) 


Example 


find root(2xx^2-1,x,0,1); 


Note 1: the command finds a value of the variable at which the 


expression is zero. 


Note 2: the values that the expression takes at the endpoints of the 


interval must have different signs. 


Computer activity 27 Solving equations numerically 


(a) Find a numerical approximation to a solution of the equation 
x? + 2x — 1 = 0 seen in Computer activity 24(e) as follows. 


(i) First, define the function f in Maxima, with the rule 
f(z) 5 z? + 2x — 1. 


(ii) Next, plot the graph of f(x) over the interval —1 € x < 1, say. 


This will help you to find an interval that contains a 


solution in the interval (0, 1). 


so on, until you find a suitable interval, if one exists. 
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From the graph, you can see that f(0) < 0 and f(1) > 0, so there is a 


If there were no solution in the interval over which you plotted the 
function, then you could plot the function over a larger interval, and 


solution. 


(iii) Find a numerical approximation to the solution of the equation in 
the interval (0, 1) by entering 


find_root (f(x), x, 0, 1); 


Remember to type the underscore symbol (_) in the middle of the 
command name. 


Notice that the first argument of the find_root command is an 
expression (or a function defined by an expression) not an equation. 
Maxima finds an input value at which the expression takes the 
value 0. 


(b) Using similar methods, find a solution of the equation cos x = yz. 


Hint: Consider the function g(x) = cos z — yz. 


3.4 Calculus 


The Maxima commands for differentiating and integrating expressions are 
diff and integrate. These are demonstrated in the following activity. 


Note that sometimes Maxima may give the result of differentiating or 
integrating an expression in a different, but equivalent, algebraic form to 
the result that you might obtain when differentiating or integrating by 
hand. If this happens, then you might like to check that the two answers 
are equivalent. You can do this by finding the difference between the two 
answers and checking that it simplifies to zero (in the case of 
differentiation) or a constant (in the case of integration). 


A rougher check is to use Maxima to plot graphs of the two answers and 
check they appear to be the same (for differentiation) or vertical 
translations of each other (for integration). 


3 Using Maxima 
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Computer activity 28 Differentiating and integrating 


COST conten 
a enterin 
Ps” R 


diff(cos(x)/(x^243), x); 


(a) Find the derivative of f(x) = 


The first argument of the diff command is the expression to be 
differentiated, and the second argument is the variable to differentiate 
with respect to. 


You must include the second argument. If you do not, then the 
answer returned by Maxima will be the derivative multiplied by 
del(x). This notation is outside the scope of MST125. 


As you would expect, the first argument can also be a variable that 
you have previously defined to be an expression, or a function you 
have previously defined. 


(b) Find the second derivative of g(u) = u? 


diff(u^2*sin(u), u, 2); 


sinu, by entering 


The third argument of diff specifies the number of times that the 


expression is to be differentiated. 


l 
(c) Find the third derivative of h(t) = E 
t 


Remember, the Maxima function for In is log. 


(d) Evaluate / z sin?(z) dx by entering 


integrate(x*sin(x)^2, x); 


The first argument of integrate is the expression to integrate, and 
the second is the variable to integrate with respect to. 


'The output returned by Maxima does not include an arbitrary 
constant. That is, Maxima gives an antiderivative of the expression, 
not its indefinite integral. When you write down an indefinite integral 
that you have found using Maxima, you need to remember to add an 
arbitrary constant yourself. 


1 
(e) Evaluate the definite integral 7 tet dt by entering 
0 


integrate(t^2*4e^t, t, 0, 1); 
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'The third and fourth arguments of integrate are the lower and 
upper limits of integration respectively. 


Notice that, as usual, Maxima returns an exact result. Use the float 
command to find a decimal approximation if required. 


2 
1 
(f) Evaluate the definite integral I — ds 
1 T 


1 
(g) Try to evaluate f tan(e *) de 
0 


Maxima cannot calculate this integral, so it is displayed unevaluated. 


As you saw in Computer activity 28(g) it is not always possible to express 
the result of an integral exactly. Sometimes, Maxima may give a result in 
terms of the so-called special functions. These cannot be simply 
expressed in terms of functions you have already met, but many can be 
defined in terms of antiderivatives of simpler functions. Such functions 
include erf, which is known as the the Error function, and 
gamma_incomplete, the incomplete gamma function. 


3 Using Maxima 


You will revise how to calculate numerical approximations to definite 
integrals that cannot be calculated symbolically later in this section. 


The diff and integrate commands are summarised in the following box. 


Calculus 
Operation 


Differentiate 


Differentiate multiple times 


Integrate 
(find an antiderivative) 


Evaluate a definite integral 


Command 

diff( , variable) 

diff I, ESSI, 
positive integer) 


integrate( , variable) 


integrate( , variable, 
lower limit ,upper limit) 


Example 


diff(sin(x^2),x); 


diff (log(x) ,x,3); 


integrate (x%2,x); 


integrate(sin(x) ,x,0,1); 


If you wish to assign the result of a diff or integrate command to a 
named function, then you need to take care that the function is defined to 
be equal to the result of the operation, not the operation of differentiating 


or integrating itself. 


This is demonstrated in the following activity. 
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Computer activity 29 Defining a function to be the derivative or 
integral of another function 


(a) Define the function f(x) = zV1-4z?. 
(b) Try to define d£ (x) to be the derivative of f(x), by entering 
df GO :=diff (f(x) ,x); 


(c) Try to find the value of the derivative of f at x = 3 by evaluating 
df (x) at x = 3; that is, by entering 


df(3); 


This gives an error. The problem is that when you entered the 
command 


df (x) :=diff (f(x) ,x); 


Maxima did not perform the differentiation. Instead, it just noted 
that df (x) means diff (f(x) ,x). Then, when you entered df (3), 
Maxima understood this to mean diff (f£ (3) ,3), which is 
meaningless, as you cannot differentiate with respect to a particular 
number. 


To correctly define df (x) to be the result of differentiating f(x), you 
need to force Maxima to perform the differentiation at the point 
df (x) is defined. 


You can do this using the syntax ''()) (which contains two single 
quotation marks, not a double quotation mark, and a pair of round 
brackets). This is demonstrated below. 


(d) Define df (x) to be the result of differentiating f with respect to x, by 
entering 


df (x) :='' (diff (f(x),x)); 


Remember that '' is two separate quotation marks. Use the upright 
quotation mark on your keyboard ('), not the sloping one (^). Don't 
forget to include the outer round brackets too. 


Notice that this time the output line shows that df (x) is indeed 
defined to be the result of differentiating f (x). 


(e) Find the value of the derivative of f at x = 3 by entering 
df(3); 


(f) Use a similar approach to define F(x) to be an antiderivative of f(x). 


(g) Calculate F(1)-F(0) and confirm this is the same value as Ja f(x)dz. 
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The command used to force the evaluation of another command is 
summarised as follows. 


Forcing the evaluation of a command 


Force the evaluation of a | ' '( ) g(x) :-2''(diff(x^4,x)); 
command 


Assumptions 


Sometimes you may want to integrate a function that contains a variable 
other than the one you are integrating with respect to. For example, the 


integral | x” dx contains the variable n as well as the variable x. When 


you use Maxima to find such an integral, you sometimes need to provide 
information about the other variable as shown in the next activity. 


Computer activity 30 Giving Maxima more information 


(a) Try to use Maxima to find the integral f x” dz, as follows. 


(i) First, ensure that n has no predefined value, by entering 


kill(n); 


(ii) Then, enter 


integrate(x^n, x); 


Maxima asks the question: Is n+1 zero or nonzero? 


This is because the answer to / x" dx depends on the value of n. 


1 
If n = —1 (so n + 1 = 0), the integral is i — dz, which is equal to 
z 


In |x| + c, whereas if n is any other value, then the integral is equal to 
qnl 
+c. So Maxima cannot calculate the integral without knowing 


n+1 
more about n. 


Enter the answer n; 
(which is short for nonzero) 


You could also answer the question with z; to mean zero. 


Maxima displays the expected result. 


3 Using Maxima 
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3 Using Maxima 


(b) An alternative approach is to give Maxima information about n before 
you use the integrate command. You can do this as follows. 


(i) Tell Maxima that n 4 —1 by entering 


assume (notequal(n,-1)); 


The assume command tells Maxima information about variables. For 
example, assume (a>0) ; tells Maxima that a is positive. 


Maxima has no simple symbol for *not equal to', so you have to use 
the notequal command for this, as above. 


(ii) Now integrate x” using 


integrate(x^n, x); 


No question is asked this time, since Maxima knows all it needs to 


know about n. 


(iii) Enter 
facts(); 


to see all the additional facts about variables known to Maxima. 


This lists all such facts. The only one shown should be the one you 
gave earlier. 


(iv) Tell Maxima to forget about the assumption on n by entering 


forget(notequal(n,-1)); 


(v) Type facts O ; again to see all the facts now known. 


No facts are displayed, since Maxima has now forgotten the 


assumption about n. 


The commands used in the previous activity are summarised below, 
together with the equal command, which is used in a similar fashion to 
notequal. 
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Facts about variables 


Operation Command Example 
Make an assumption about a | assume (F) assume (a>0) ; 
variable 
State two things are equal equal( , ) assume (equal (n,3)); 
State two things are not equal | notequal( , ) assume (notequal(n,-1)); 
Forget a property forget (property) | forget(a>0); 
List all known facts facts facts O ; 
List all known facts about facts(variable) | facts(a); 
a particular variable 


Quadrature 


You saw in Computer activity 28(g) that Maxima cannot always find an 

antiderivative for a function. However if Maxima cannot find an 

antiderivative of a function f, it can often find an approximate value for a 
b 


definite integral of the form 1 f(x) dx by calculating a sophisticated 


numerical approximation to the area between the graph of the function 
and the x-axis. The command for doing this is quad. qags, which is 
demonstrated in the following activity. 


The ‘quad’ part of the name of this command arises from the fact the 
process of finding the value of a definite integral (or any area) is sometimes 
called quadrature. The letters q, a, g and s in the second half of the 
name specify the particular method Maxima uses to find the approximate 
value of the definite integral. This method gives good results for a wide 
variety of integrals. 


Computer activity 31 Finding approximate values of definite 
integrals 


1 
(a) Try to calculate f In(sin(z?)) dx by entering 
0 


integrate(log(sin(x^3)), x, 0, 1); 


'The integral is returned unevaluated, since Maxima cannot 


calculate it. 


(b) Find an approximate value of the integral by entering 


quad. qags (log(sin(x^3)), x, 0, 1); 
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'The quad. qags command calculates an approximate value for a 
definite integral. 


Notice that the output of the quad. qags command is a list. The first 
element of the list is the approximate value of the definite integral, the 
second element is an estimate of the accuracy of the approximation, 
the third element is the number of values at which the function to be 
integrated was evaluated during the calculation, and the final element 
is an error code. An error code of 0 means that no errors occurred. 


Finding approximate values for definite integrals 


Find an approximate value | quad. qags( , variable, quad.qags (exp(-x^2) ,x,0,1); 
of a definite integral lower limit , upper limit) 


If you began studying this section from Activity 46 of Unit 1, then this 
completes your study of the unit. 


4 Number theory 


Many of the calculations that you have seen in Unit 3 can be performed 
quickly and accurately by a computer, which is particularly useful when 
the numbers are large. Maxima is ideally suited to this, since it can do 
exact calculations with extremely large numbers. However, for really large 
numbers some calculations may take a long time! 


'The Maxima commands for finding the quotient and remainder when one 
integer is divided by another, and the highest common factor (HCF) of two 
integers, are demonstrated in the following activity. 


Computer activity 32 Finding quotients, remainders and HCFs 


(a) Find the quotient on dividing 59 by 8, by entering 
quotient(59,8); 


(b) Find the remainder on dividing 59 by 8, by entering 


remainder(59,8); 


(c) Find the quotient and remainder when 3!9? is divided by 21°, 
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(d) Find the highest common factor of 108 and 93, by entering 
gcd(108,93) ; 


'The name of the Maxima command for the HCF is gcd, because an 


alternative name for the highest common factor is ‘greatest common 
divisor’. 


(e) Find the highest common factor of 31% and 21% + 11. 


You know from Unit 3 that if a and b are integers, not both 0, and d is 
their highest common factor, then there exist integers v and w such that 
av + bw = d. This is known as Bézout’s identity. Maxima has a command 
for finding d, v and w for a given a and b, which is demonstrated in the 
following activity. 


Computer activity 33 Using Bézout's identity 


(a) Find the highest common factor d of 108 and 93, and two integers v 
and w such that 108v + 93w = d, by entering 


gcdex(108,93); 


The result of this calculation is a list: 


[-6,7,3] 


The last number in the list is the highest common factor of the 
numbers, which agrees with the result of Computer activity 32(d). 
The first two numbers are the required integers v and w, so 


108 x (-6) +93 x 7=3. 
= ee 


(b) Find the highest common factor d of 309 and 210? + 11, and two 
integers v and w such that 3199y + (2100 + 11)w = d. 


The commands introduced in this section are summarised as follows. 


4 Number theory 


59 


5  Parametric equations 


Quotients, remainders and highest common factors 


Operation Syntax 

Quotient on dividing one number by another quotient( , ) 
Remainder on dividing one number by another | remainder(| , ) 
Highest common factor gcd (g. lm 


Example 


quotient(6,4); 
remainder(6,4); 
gcd(6,4); 


(greatest common divisor) of two numbers 
Bézout's identity: given integers a and b gcdex( , ) 
find the highest common factor d and 

integers v and w such that av + bw — d 


If you began studying this section from Activity 6 of Unit 8, now return to 
the unit to continue your study. 


5 Parametric equations 


In this section you will learn how to use Maxima to plot curves given by 
parametric equations. 


5.1 Plotting curves from parametric equations 


Subsection 3.2 of this Guide reviewed how to use the wxplot2d command 
to plot graphs of equations given in the explicit form y — f(x), and how to 
use the wximplicit.plot command to plot graphs of equations given in 
implicit form. In this subsection, you will learn how to use the wxplot2d 
command to plot curves given by parametric equations. 


Note that if you are using a Maxima interface other than wxMaxima, you 
should use the plot2d command instead of wxplot2d. The two commands 
are used in similar ways. 


Computer activity 34 Plotting a curve from parametric equations 


(a) Plot the curve given by the parametric equations 
=t+4, y=2t+3 
for values of the parameter t in the range —1 < t < 2 by entering 


wxplot2d( [parametric, t+4, 2*t+3, [t,-1,2]] ); 
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gcdex (6,4); 


Here, the first argument of the wxplot2d command is a list whose 
first element is the keyword parametric. The second and third 
elements of the list are the parametric expressions for x and y 
respectively. The fourth element of the list is itself a list containing 
the name of the parameter followed by the minimum and maximum 
values of the parameter for which the curve should be plotted. It is 
important that this list is within the list beginning parametric. 


The curve plotted is part of the straight line that you plotted in 
Activity 23 of Unit 4. 


It is a line segment corresponding to the specified range of values of t, 
namely —1 < t € 2. Maxima automatically chooses ranges for the x- 
and y-axes to ensure that the whole of the line segment is visible. To 
see more clearly that the curve is a line segment rather than a whole 
line, you can change the ranges of the x- and y-axis, which you can do 
using the usual Maxima graph plotting options as shown in part (b) 
below. 


Notice also that the expressions for the x- and y-coordinates of the 
points on the line have been automatically added as axis labels. You 
can change these using the standard graph plotting options, as also 
shown in part (b). 


(b) Plot the line from part (a) again, this time with the x- and y-axes 
having the ranges 0 € x € 8 and 0 € y € 8, and with the axes labelled 
‘x’ and ‘y’ respectively, by entering 
wxplot2d( [parametric, t+4, 2*t+3, [t,-1,2]], 
[x,0,8], [y,0,8], [xlabel,"x"], [ylabel,"y"] ); 


'The graph now clearly shows the line segment described by the given 


equations and range of t. 


(c) Plot the curve given by the parametric equations 


r=3°, y=6t (24:22) 
This is the curve that you plotted by hand in Activity 24 of Unit 4. 


When you are plotting curves from parametric equations, it is often 
convenient to assign the expressions for x and y to variables or functions, 
and use these in the wxplot2d command. Note, however, that you should 
not assign them to the variables x and y, as this will give errors when you 
specify axis ranges. You can, however, assign them to the functions x and 
y, as demonstrated next. 


5 Parametric equations 
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5 Parametric equations 


Computer activity 35 Using functions to plot curves from 
parametric equations 
(a) Plot the curve given by 
r—12—12sim?t, y=12cost (0<t<m), 
as follows. 


(i) Define the function x(t) to be 12 — 12sin? t and the function 
y(t) to be 12cost. 


Remember, to enter sin? t, type sin(t)^2. 


(ii) Plot the curve, using the z-axis range —5 € x < 15 and the y-axis 
range —15 < y < 15, by entering 
wxplot2d( [parametric, x(t), y(t), [t,0,%pil], 
[x,-5,15], [y,-15,15], [xlabel,"x"], [ylabel,"y"] ); 


This is the same curve as plotted in Computer activity 34(c). It is 


possible to represent the same curve with different parametric 
equations. 


(b) By first defining the functions x(t) and y(t), plot the curve given by 
x=tcost, y=tsint (051 < 7). 


Let Maxima choose appropriate axis ranges. 


In some cases, the default settings for plotting curves from parametric 
equations result in plots that are not very representative of the actual 
curve. This is demonstrated in the following activity, where the technique 
needed to resolve this issue is explained. 


62 


Computer activity 36 — /nvestigating the number of points plotted 


(a) Plot the curve given by 
r—10, y=sin(10t) (0<t<5) 
as follows. 


(i) Assign the expressions for x and y to the functions x (t) 
and y(t). 


(ii) Plot the curve by entering 


5 Parametric equations 


wxplot2d( [parametric, x(t), y(t), [t,0,5]], [xlabel,"x"], [ylabel,"y"] ); 


The curve plotted is not very smooth. 


To plot a curve from parametric equations, Maxima calculates points 
corresponding to values of the parameter spread throughout the 
specified range, and joins them with straight lines. By default, it uses 
29 points. In this example, this is not sufficient to adequately 
represent the curve. The individual straight lines joining the points 
are easily seen. 


You can change the number of points plotted by including another 
argument of the wxplot2d command, namely a list whose first 
element is the keyword nticks and whose second element is the 
number of points required. 


(iii) Plot the curve again, this time using 200 points, by including the 
additional argument [nticks,200] in the wxplot2d command. 


(b) Plot the curve given by 
x = tsin(2t), y= tcos(2t) (0<t<10), 
ensuring that a smooth curve is obtained. 
(c) Plot the curve given by 
x = sin(2t), y=cos(3t) (0< t< 20) 


ensuring that a smooth curve is obtained. 


This curve is an example of a Lissajous curve. In general, a Lissajous 
curve is a curve given by parametric equations of the form 
x = Asin(at + d), y= Bsin(bt) 


where A, B, a, b and d are numbers. Different values of these numbers 
give different curves. If a/b is a rational number, then a closed curve 
is obtained, provided it is plotted for a sufficiently large range of t. 
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You can plot a curve from parametric equations on the same diagram as 
another curve from parametric equations, or on the same diagram as any 
other graph. You do this by listing the curves to be plotted as the first 
argument of the wxplot2d command. In this case, to prevent the 
commands getting too long, it is often helpful to assign the list defining a 
parametric curve to a named variable, as demonstrated below. 


Computer activity 37 Plotting more than one curve 


(a) Plot the curve given by 
r—30, y=8t? (-1<t<1) 
and the curve that is the graph of the function 
fix) = 22? 
on the same diagram, as follows. 


(i) First, assign the expressions for x and y for the first curve to the 
functions x(t) and y(t), and define f (x) to be 227. 


(ii) Assign the list required to specify the first curve to the variable p 
by entering 


p:[parametric, x(t), y(t), [t,-1,1]]; 


(iii) Plot the curves by entering 
wxplot2d( [p,f(x)], [x,-2,2], [nticks,50] ); 


Here, the Maxima specifications of the two curves are given in a list 
as the first argument of wxplot2d. This is a similar syntax to that 
used when plotting the graphs of more than one function together. 


Note that since one of the plots is specified by a function, a range of 
values of x over which to plot is needed, as this is a requirement of 
wxplot2d when plotting functions. This must be given immediately 
after the list of curves to plot, as the second argument of the 
wxplot2d command. 


(b) Plot on the same diagram the two curves given by 
r-—sin?t, y=cos*t (0<t< 2r) 
and 


z—cos*tsint, y-sin?tcost (0<t< 2m) 


The use of the wxplot2d command to plot curves from parametric 
equations is summarised as follows. 
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Plotting curves from parametric equations 


Plot a curve from wxplot2d([parametric, wxplot2d([parametric, 


parametric equations expression for x, expression for y, 2*t, 3*t*1, [t,0,1]]); 
[parameter , min, max]], ...) 


Note: the parameter can be any variable, but x and y should not be 
used. 


Graph plotting options 


Change the number of points plotted ona | [nticks, |] [nticks, 100] 
curve from parametric equations 


If you began studying this subsection from Activity 25 of Unit 4, now 
return to the unit to continue your study. 


5.2 Plotting conics parametrically 


You can plot conics parametrically using the wxplot2d command as 
described in Subsection 5.1 of this Guide. 


However, you need to take care when plotting a hyperbola, since the 
expressions for z and y in the standard parametrisation of a hyperbola in 
standard position which is of the form 


T 
x =asect, y=btant GIL EE) 


are not defined at the endpoints t = —7/2, t = 1/2 and t = 37/2 of the 
ranges of values of t. If you try to use Maxima to plot a hyperbola using a 
parametrisation of this form, then the system will attempt to evaluate the 
expressions at values of t near —7/2, 7/2 and 37/2, and will obtain values 
with very large magnitudes, resulting in a misleading graph. You can 
correct this by imposing suitable ranges for the x- and y-axes, as 
demonstrated in the next activity. 


Another advantage of imposing suitable ranges for the z- and y-axes is 
that, as demonstrated in the next activity, you do not have to plot two 
curves, one for —7/2 < t < 7/2 and one for 7/2 < t < 37/2. Instead you 
can plot a single curve for —7/2 < t < 37/2. This is because the imposed 
ranges for the x- and y-axes prevent Maxima displaying points whose 
coordinates have large magnitudes arising from evaluating the expressions 
at values of t near 7/2. 
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Parametric equations 


Computer activity 38 Plotting conics parametrically 


(a) Plot the hyperbola given by 


3+2 t 5+ tant Aude T guo Bt 
=> = — an =— = — — 
x Sect, y 2 2" 5) 2 


by entering 
wxplot2d( [parametric, 3+2*sec(t), -5ttan(t), [t,-pi/2,3*pi/211, 
[xlabel , "x"] 7 [ylabel , "spit ) ; 


This is the hyperbola considered in Activity 40(a) of Unit 4. 


The graph is distorted by the large values calculated near t = - 


T 3T 
t= audiet. 
a 2 


(b) Plot the hyperbola in part (a) again, this time restricting the x- and 


y-axes to —2 < x < 8 and —10 € y < 0. Make sure that sufficient 
points are plotted to produce a good representation of the curve. 


(c) Plot the hyperbola in parts (a) and (b) a third time, now including the 


E: 
asymptotes, which are y — and y — 237 5 9n the diagram. 


x 
2 2 


(d) Plot the parabola given by 


r=-—1+4, y=4+8t (—2<t<2) 


This is the parabola considered in Activity 40(b) of Unit 4. 


(e) Plot the ellipse given by 


x=1+2cost, y=2+3sint (0<t< 27) 


If you began studying this subsection from Activity 41 of Unit 4, this now 
completes your study of the unit. 
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6 Partial fractions and polynomial division 


6 Partial fractions and polynomial 
division 


In this section you will learn how to use Maxima to find partial fraction 
expansions of rational expressions, and to find the quotient and remainder 
when you divide one polynomial expression by another polynomial 
expression. 


6.1 Partial fractions 


'The Maxima command for finding the partial fraction expansion of a 
rational expression is partfrac. This command is demonstrated in the 
following activity. 


Computer activity 39 Finding a partial fraction expansion 


(a) Find the partial fraction expansion of the rational expression 
Ta 


(x + 5)(z — 2) 
partfrac( 7*x/((x+5)*(x-2)), x ); 


by entering 


'The first argument of the partfrac command is the rational 
expression that you want to write in terms of its partial fractions. The 
second argument is the variable used in the rational expression. All 
the rational expressions we consider here contain only one variable. 


Notice that in this example the denominator of the rational 
expression is enclosed in brackets, since without the brackets 
Tx/ (x*5) *(x-2) would be interpreted as 
7x(a — 2) 
r+5 ` 
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(b) Confirm that the partial fraction expansion is equal to the original 
expression by entering 


factor (%); 


This command adds together the partial fractions obtained in 
part (a) and expresses the result with a factorised denominator. 
Remember that % refers to the result of the last evaluated command. 


The output is the rational expression that you started with. 


The use of the partfrac command is summarised below. 


Partial fraction expansion of a rational expression 


Find the partial fraction | partfrac(expression, variable) | partfrac(1/(x*2-1) ,x); 


expansion of a rational 
expression 


In the following activity you can use the partfrac command to see the 
form of the partial fraction expansions of various rational expressions. 


Computer activity 40  /nvestigating partial fraction expansions 


Find the partial fraction expansions of the following rational expressions. 


How are the denominators of the partial fractions related to the 
denominators of the original expressions? 
9 3x 4r +8 

A b) 2% ASEO 
(a) (2x + 1)(z — 1) (b) ge+an—2 (c) (4x — 3)? 

8a +4 26 
do CU MERE CENE ee 
(d) (5x — 3)(a — 1)? (e) (a 4- 5)(z? +1) 


Notice that all the partial fractions that you obtained in parts (a)-(e) 
are of the form 


A Ax+B 
— o — 
(ax +b)” (ax? + bx +c)” 
where n is a positive integer and a, b, c, A and B are constants, as 
described in Unit 7. 


If you began studying this subsection from Activity 8 of Unit 7, now return 
to the unit to continue your study. 
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6.2 Polynomial division 


You can find the quotient and remainder on dividing one polynomial 
expressions by another using the same commands that you use to find the 
quotient and remainder when dividing integers, which you met in Section 4 
of this Guide. The use of these commands for polynomial expressions is 
summarised below and is demonstrated in the activity that follows. 


Polynomial division 
Operation Command Example 


Quotient on dividing one polynomial quotient( , ) quotient(x^2,x-1); 
expression by another 


Remainder on dividing one polynomial | remainder( , ) remainder(x^2,x-1); 
expression by another 


Computer activity 41 Dividing polynomial expressions 
(a) Find the quotient and remainder on dividing z? + x +1 by z? 4- 2, by 
entering 
quotient (x*5+x+1, x^342); 
then 
remainder (x*5+x+1, x%3+2); 
(b) Find the quotient and remainder on dividing a(x) by b(x) for each of 
the following pairs of polynomial expressions. 
(i) a(x) =3zrt— 72, b(a)=x=-—1 
(ii) a(z)=2°+2°, d(x) =22741 


If you began studying this subsection from Activity 16 of Unit 7, now 
return to the unit to continue your study. 


7 Differential equations 


In this section, you will learn how to use Maxima to find analytic solutions 
of first-order differential equations (where this is possible), plot direction 
fields of first-order differential equations, and calculate approximate 
numerical solutions of first-order differential equations. 
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7 Differential equations 


7.1 Solving differential equations analytically 


The Maxima command for solving a differential equation analytically is 
ode2. Here, ‘ode’ stands for ordinary differential equation. All the 
differential equations that you have met in Unit 8 are so-called ordinary 
differential equations, which simply means that the unknown function is a 
function of only one variable. Equations in which the unknown function is 
a function of more than one variable are called partial differential equations 
and you may meet these in higher-level modules. The ‘2’ in the command 
name means that this command is capable of analytically solving (where 
possible) differential equations of order up to 2; that is, it can solve first- 
and second-order differential equations. However, as in Unit 8, we will 
consider only first-order differential equations here. 


When you enter a differential equation in Maxima, you need to use a 
special syntax to enter the derivative in the equation. Consider, for 
example, the differential equation 

e y = e” sing 

dz ' 
which you met in Example 11 of Unit 8. You might think that you can 
enter the term dy/dx by typing diff (y,x), since the command for 
differentiating an expression in Maxima is diff. (This is revised in 
Subsection 3.4 of this Guide.) However, entering diff(y,x) in Maxima 
gives the result 0. You might like to try this. The reason is that when you 
enter diff(y,x), Maxima will try to perform the differentiation, and since 
there is nothing to tell it that y is a function of x, it will obtain the 
result 0. 


To enter the operation of differentiation in Maxima, while preventing 
Maxima from trying to actually perform the operation, you need to prefix 
the diff command with a single quote mark ('). You may recall (again, as 
revised in Subsection 3.4 of this Guide) that two quote marks are used to 
force the evaluation of a command. Here, one quote mark is used to 
prevent the evaluation of a command. 


The use of this method of entering differential equations and the ode2 
command are demonstrated in the following activity. 


Computer activity 42 Solving differential equations analytically 


(a) Solve the differential equation 
dy 

dx 

as follows. 


—y=e" sina 
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(i) First assign the equation to the variable eqn by entering 


eqn:'diff(y,x) - y = exp(x)*sin(x); 


Notice the use of 'diff(y,x) for the derivative dy/dx, to prevent 
Maxima from performing the differentiation immediately. 


The quote mark to use is the upright one (') on your keyboard, not 
the sloping one (^). 

d 
wxMaxima displays 'diff(y,x) as m when showing the result of 


this command. Also, the expression exp(x) is displayed as %e*, where 
^e is Maxima’s notation for the constant e. 


(ii) Now solve the equation by entering 


ode2(eqn, y, x); 


The first argument of the ode2 command is the equation to be solved, 
or, as in this case, the name of a variable representing that equation. 
The second argument is the name of the dependent variable, and the 
third argument is the name of the independent variable. 


'The result of this command is the general solution of the differential 
equation. The symbol 4c is used to represent the arbitrary constant. 
Notice that the solution is an equation relating the dependent 
variable, y, to the independent variable, x. 


The solution should be the same as that obtained in Example 11 of 
Unit 8, namely y = e*(— cos x + c). 


(b) Use Maxima to solve the differential equation 


d 
gt +y = ry’ 
dz 


You might like to use the fullratsimp command to simplify the 


result. 


(c) Use Maxima to try to solve the differential equation 


dy COS T 
I$ o = | 
dx E 


This differential equation is directly integrable, but Maxima cannot 


calculate one of the integrals needed for the solution. The answer is 
displayed with this integral unevaluated. 
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(d) Use Maxima to try to solve the differential equation 


d 3 
E = —1 
dt i 


Notice that the independent variable is t, not x. 


The result of the ode2 command is false, which means that Maxima 
is unable to solve the equation using this command. 


The ode2 command uses a range of methods to try to find an analytic 
solution of a differential equation. These include the methods that you 
have studied in Unit 8 and some others. However, the command cannot 
always find an analytic solution to an equation, as you saw in Computer 
activities 42(c) and (d). One reason for this is that not all differential 
equations have analytic solutions. In Subsections 7.2 and 7.3 of this Guide 
you will see some ways in which you can investigate solutions if analytic 
solutions cannot be found. 


As you have seen in Unit 8, having found the general solution to a 
differential equation, you often need to find a particular solution that 
satisfies a given condition. The command ic1 can do this, as demonstrated 
in the following activity. The ‘ic’ in the name of this command stands for 
‘initial condition’, which is how we often describe the additional 
information needed to find a particular solution. The number ‘1’ means 
that this command is appropriate for first-order differential equations, 
where one initial condition is needed to allow you to find a value for the 
one arbitrary constant appearing in the general solution. 


Computer activity 43 Solving an initial value problem 


Find the solution of the initial value problem 


dy 
dt 
as follows. 


t— +2y =t, where y(1) — 1, 


(a) Assign the equation to the variable eqn. 


(b) Solve the equation, and assign the solution to the variable sol, by 
entering 


sol:ode2(eqn, y, t); 


(c) Find the required particular solution of the equation by entering 


ici(sol, y-1, t-1); 
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The given initial condition, y(1) — 1, is equivalent to saying that 
y —1 when t= 1. 


The first argument of the ic1 command is the general solution 
previously calculated. The second and third arguments give the initial 
condition by specifying corresponding values of the dependent and 
independent variables, in either order. These values are given by 
equations of the form 


variable name = value 


The solution obtained should be the same as the one you found by 
hand in Unit 8, Activity 27(b), namely y = t?/4 + 3t ?/4. 


The commands needed to solve differential equations analytically are 
summarised below. 


Solving differential equations analytically 
Operation Command Example 


Prevent the evaluation : 'diff(y,x) 
of a command 


Solve a differential equation | ode2( equation , dependent variable, ode2('diff(y,x)-x, y, x); 


analytically independent variable) 
Apply initial conditions to ici(solution ,variable- value, ici(sol, y=1, x=0); 
the analytic solution of a variable - value) 


differential equation 


Here is another problem for you to solve using Maxima. 


Computer activity 44 Solving another initial value problem 


(a) Solve the initial value problem 


Ln y—(x—2), where y(0) = 1.9. 
dx 


(b) Plot a graph of the solution over the range 0 € x < 5. 


(c) Find the exact value of y given by the solution when x = 4. 
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7.2 Plotting direction fields 


In Unit 8 you saw that one way of investigating the solutions of a 
differential equation when it is not possible to find analytic solutions is to 
plot a direction field, as illustrated in Figure 11. 


A direction field of a differential equation of the form 


where f(x,y) is a known function, is a set of short line segments equally 
spaced in the x, y-plane. The gradient of each segment is equal to the value 
of f(x,y) at the midpoint of that segment. The direction field gives you an 
idea of the shape of the graph of any particular solution of the differential 
equation. This is because at each point (x,y) on the graph of a particular 
solution the gradient is f(x,y). So where such a graph passes through the 
midpoint of a line segment, the graph will have the same gradient as the 
line segment, and where the graph passes close to a line segment, you can 
expect the graph to have nearly the same gradient as the line segment. 


The Maxima command for plotting a direction field is wxdrawdf, where 
‘df’ stands for ‘direction field’. This command is not part of the basic 
Maxima collection of commands, so before using it you have to load an 
additional Maxima package, the drawdf package. This package needs to be 
loaded only once in each Maxima session. Details of how to load the 
package are given in Computer activity 45. 


If you are using a Maxima interface other than wxMaxima, then you 
should use the drawdf command instead of wxdrawdf. If you have 
configured Maxima to automatically use different colours and styles for 
graphs, as described in the Accessibility section of the OU Maxima website, 
then you should use the commands ouwxdrawdf (if you use wxMaxima) or 
oudrawdf (otherwise) instead. All these commands are used in a similar 
way to wxdrawdf, and all of them require the drawdf package to be loaded. 


Note that the wxdrawdf command uses a different drawing system to the 
wxplot2d and wximplicit. plot commands that you have used previously, 
so many of the options that you have used with wxplot2d and 
wximplicit plot do not work with wxdrawdf. 


The use of the wxdrawdf command is demonstrated in the following 
activity, where you will plot the direction field of the differential equation 
that you solved analytically in Computer activity 44. 
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Figure 11 A direction field 
for dy/dx = sin(xy) 


7 Differential equations 


Computer activity 45  /nvestigating a direction field 


(a) First load the drawdf package by typing 
load(drawdf); 


Maxima will produce some output as a result of this command, which 
may vary depending on whether you have ever used this package on 
your computer before. 


The output will end with a line similar to 


(%01) C:/PROGRA~1/MAXIMA~1.0/share/maxima/5.30.0/share/ 
diffequations/drawdf.mac 


which indicates the location of the package on your computer. 


(b) Plot the direction field of the differential equation 


dy _ 2 


by entering 
wxdrawdf(y-(x-2)^2, [x,y], field .arrows-false); 


The first argument of the wxdrawdf command is the expression for 
d 
the gradient i at each point (x,y). 
T 


The second argument is a list containing the names of the two 
variables in the differential equation. In this case these are x and y. 
'The first variable listed is plotted along the horizontal axis and the 
second is plotted along the vertical axis. Note that if, as in this case, 
the variables are x and y and you want to plot x on the horizontal 
axis and y on the vertical axis, then this argument can be omitted. 


'The third argument, field arrows-false, indicates that you want to 
plot a set of line segments not arrows. As you saw in Unit 8, a 
direction field is drawn using line segments. In other applications, 
arrows are needed. 


Notice that the default ranges of the values of the variables plotted by 
the wxdrawdf command are from —10 to 10 for both variables. 


From the direction field displayed you can see that for much of the 
region plotted the line segments are nearly vertical. Something more 
interesting seems to be happening in one particular region of the 
direction field however, so in the next part of this activity you can 
change the region plotted to concentrate on this area. 
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Differential equations 


(c) Plot the same direction field, this time limiting the range of values of x 


c 


and y to —3 < x < 6 and —5 < y < 10 by entering 
wxdrawdf(y-(x-2)^2, [x,-3,6], [y,-5,10], field.arrows-false); 


(You may like to edit and re-evaluate the command that you entered 
in part (b) rather than entering the command again.) 


The ranges of values of x and y are specified by including additional 
arguments each of which is a list containing the name of the variable 
and the minimum and maximum values of this variable over which 
the direction field should be plotted. This way of specifying ranges is 
identical to the way that you have previously specified ranges of 
values when using the wxplot2d and wximplicit.plot commands. 


Note that when you specify ranges of values for both variables, you 
can omit the argument listing the variables, even if they are not x 
and y. 


In the plot, you can more easily see details of the behaviour of 
solutions of the differential equation in the region shown. 


(d) Use the direction field to try to work out what the graphs of particular 


solutions of the differential equation might look like. 


'To do this, start at a particular point on the direction field and move 
in the direction given by the line segment at that point. When you get 
to the next column of line segments, continue moving but in the 
direction suggested by the nearest line segments. Continue in this way, 
ideally ‘smoothing’ the changes of direction, until you reach an edge of 
the plot. 


You may like to print the direction field and draw such graphs of 
solutions on it. 


(i) What do you think the graph of the particular solution that 
passes through the point (0,0) might look like? 


(ii) What do you think the graph of the particular solution that 
passes through the point (0,4) might look like? 


(e) Add the following additional argument to the command entered in 
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part (c): 


[trajectory_at, 0,0] 
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'This argument adds the approximate graph of the particular solution 
that passes through the point specified. 


'The argument is a list beginning with the keyword trajectory.at. 
(Note the underscore character in this keyword.) The second and 
third elements of the list are the coordinates of the point that the 
required solution should pass through. 


You can add more than one approximate solution graph to the 
direction field by including more than one argument of this form. 


Each approximate particular solution is made up of a sequence of 
small straight line segments, which are found as follows. Maxima 
starts with the specified point. It then increases the value of the 
independent variable (in this case x) by a small amount, known as the 
step size, and calculates the corresponding change in the dependent 
variable (in this case y), by using the value of dy/dx at that point. 
This gives a new point. The process is then repeated starting with 
this new point, and repeated again a number of times, to give a 
sequence of points. Maxima then returns to the specified point and 
carries out the same process, this time decreasing the independent 
variable by the step size. The points obtained are joined up with 
short line segments. 


Maxima automatically changes the step size used when calculating 
approximate solutions depending on how rapidly the direction field is 
changing. You can set the maximum step size to be used by including 
an additional argument of the form 


[tstep, 0.001] 
This example sets the maximum step size to be used to 0.001. 


You can also change the number of steps taken by including an 
additional argument of the form 


[nsteps, 200] 
which sets the number of steps to be taken to 200. 


Note that these additional arguments only affect trajectory.at 
arguments that appear after them in the list of arguments of the 
wxdrawdf command. 


You can see that the particular solution that passes through the point 
(0,0) has y decreasing as x increases. 


(f) Add to the plot the particular solution that passes through the point 


TT 


7 Differential equations 


(g) Finally, add to the plot the particular solution that passes through the 


point (0, 1.9). 


conditions. 


This graph should match the graph of the particular solution that you 
obtained analytically and plotted in Computer activity 44. 


The solutions plotted here demonstrate that the behaviour of the 
solution to a differential equation can vary depending on the initial 


The fact that the behaviour of the solutions of a differential equation 
can be obtained from a direction field without the need to solve the 
differential equation makes a direction field a powerful tool. 


The commands needed to plot direction fields are summarised in the tables 


below. 


Plotting direction fields 


Plot the direction field of a | wxdrawdf (expression, 


differential equation 


RAS EN 


wxdrawdf (xx*y, [x,y], 
field_arrows=false) ; 


field_arrows=false, 


Note 1: the drawdf package needs to be loaded before this command 


can be used. 


Note 2: the first variable listed is plotted on the horizontal axis and 


the second variable listed is plotted on the vertical axis. 


Direction field plotting options 
Operation 


Set the range of a variable 
Plot the graph of a particular solution 


Set the maximum step size to use when 
approximating a solution 


Set the number of steps used when 
approximating a solution 


Argument 


[variable , I 


[trajectory.at, 


[tstep, 


[nsteps, 


] 


] 


] 


* 


Example 


[x,0,5] 
[trajectory.at, O, 1] 
[tstep, 0.001] 


[nsteps, 200] 


In the following activity you can use direction fields to help understand the 
behaviour of the solutions of two differential equations that Maxima is 


unable to solve analytically. 
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Computer activity 46 Using direction fields 


Plot a direction field for each of the following differential equations and 
then add the graphs of some approximate particular solutions and describe 
their behaviour. For example, you might like to start by considering the 
solution that passes through (0, 0). 


How does the behaviour of the particular solutions vary as the initial 
condition is changed? 

dy 

a — = 
(a) dz 


cos 4 


(Hint: use the default range of the direction field, —10 € x < 10, 
—10 < y € 10.) 
(This is the differential equation that Maxima was unable to solve 
analytically in Computer activity 42(c).) 
dy 

b ==ty-1 

(b) qe Y 
(Hint: restrict the range of the direction field to —2 « t < 2, 
—2<y<2.) 
(This is the differential equation that Maxima was unable to solve 
analytically in Computer activity 42(d).) 


7.3 Solving differential equations numerically 


As you saw in the previous section, you can find an approximate particular 
solution of a differential equation by using the direction field of the 


canon. Carl Runge (1856-1927) 
Such a solution is usually represented by a sequence of points that lie on 
the approximate solution. Such approximate solutions are known as 
numerical solutions. 


The Maxima command for numerically solving a differential equation of 
the form 
dy 


de fy) 


is rk. 
Here, ‘rk’ stands for Runge-Kutta. This is the name given to an 
often-used class of methods for solving differential equations numerically. 


They are named after the German mathematicians Carl Runge and 
Martin Wilhelm Kutta, who first developed them. 


The rk command is demonstrated in the following activity. 


To solve a differential equation numerically, we need both an equation and Martin Wilhelm Kutta 
an initial condition, since the approximate solution is an approximation of (1867-1944) 
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a particular solution of the equation. That is, we numerically solve an 
initial value problem. 


Computer activity 47 Solving an initial value problem numerically 


(a) Find a numerical approximation of the solution of the initial value 


problem 
— =y—(x—2)?, where y(0) = 1.9, 
for 0 € a € 5 by entering 


pts:rk(y-(x-2)^2, y, 1.9, [x, 0, 5, 1]); 


(This is the differential equation that you solved analytically in 
Computer activity 44.) 


by the equation to be solved. 


The second argument is the dependent variable (in this case y) and 


initial condition. 


The final argument is a list beginning with the name of the 


the step size to be used. The step size is the change in the 
independent variable between adjacent points in the numerical 
solution. 


the z- and y-coordinates of a point on the numerical solution. 


Notice that the x-coordinates start at 0 and increase in steps of 1 
until 5, as requested in the command entered. 


LII L———————————————————————————— —————————À 


(b) The analytic solution of the initial value problem considered in 
part (a) that was obtained in Computer activity 44 is 


ev 
p= +0 2042. 


Assign the analytic solution to the variable asol, and then plot this 


solution together with the numerical solution calculated above by 
entering 


asol:-exp(x)/10*x^2-2*x*2; 
then 


wxplot2d([asol, [discrete, pts]], [x,0,5], [legend, "Analytic", 
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d 
The first argument of the rk command is the expression for E given 
x 


the third argument is the initial value of this variable, as given by the 


independent variable (in this case, x) followed by the initial value of 
this variable, as given by the initial condition, then the last value of 
the variable for which you want to calculate the solution, and finally 


The result is a list of pairs of numbers, in which each pair consists of 


"Numerical"]); 
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Here, the first argument of the wxplot2d command is a list of the 
curves to plot. (This use of the command is revised in Computer 
activity 18 in Subsection 3.2 of this Guide.) The first element of this 
list is the analytic solution and the second element is a list that 
represents the numerical solution. This list consists of the keyword 
discrete followed by the list of points to be plotted. In the plot, 
these points will be joined by straight line segments. A legend has 
also been included. 


Notice that the numerical solution does indeed approximate the 
analytic solution, but the agreement gets worse as x increases. This is 
not surprising since only a crude numerical solution was calculated, 
using the relatively large step size of 1. 


(c) Calculate a revised numerical solution, using a step size of 0.1, and 
plot this solution together with the analytic solution. 


'The agreement between the two solutions is much improved. 


Usually, the smaller the step size in a numerical solution, the better 
the result. 


'The rk command for finding a numerical solution of a differential equation 
and the method used to plot the graph represented by a set of points are 
summarised below. 


Solving a differential equation numerically 


Solve a differential rk(expression, dependent variable, rk(x^2,y,1, 
equation numerically initial value, [independent variable, [x,0,10,0.1]); 


initial value, final value, step size]) 


Plotting graphs represented by a set of points 


wxplot2d([discrete, | wxplot2d([discrete, 
list of points]) [[0,0], [1,113 1); 


Here is an activity for you to try. 


Operation 


Plot a graph represented 
by a set of points 
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Computer activity 48 Finding more numerical solutions 


Find a numerical solution of each of the following initial value problems, 
over the range given. Plot your solutions and compare them to the 
direction fields plotted in Computer activity 46. 


d 
(a) i — 679? —] where y(0) 20, for0<x<10 
ls 


(b) E —ty!—1 where y(0)=1.5, for0<t<1 


If you began studying this subsection from Activity 37 of Unit 8, this now 
completes your study of the unit. 


8 Eigenvalues and eigenvectors 


In this section, you will learn how to use Maxima to find the characteristic 
equation, eigenvalues and eigenvectors of a matrix. Before doing that, in 
Subsection 8.1 you can revise how to enter and use matrices in Maxima. 
You learned how to do this in MST124 Unit 9, and if you are still familiar 
with it you may wish to skip this subsection. 


8.1 Matrices 


There are several ways to input matrices in Maxima. Two of these are 
demonstrated in the following activity. 


Computer activity 49 Entering matrices 


2 —3 3 
(a) Assign the matrix | -2 1 2] to the variable P by entering 
1-1 4 


P:matrix( [2,-3,3], [-2,1,2], [1,-1,4] ); 
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'The matrix command allows a matrix to be defined. You give each 
row of the matrix as a list (that is, within square brackets with the 
elements separated by commas). 


Notice that wxMaxima displays matrices using square brackets rather 
than the round brackets used in the unit: 


2 -3 3 
=2 4 2 
1 -1 4 


If you are using a different interface to Maxima, or have chosen 
different display settings, matrices may be displayed as follows. 


2 =3 3 


[ ] 
[ ] 
[ -2 1 2 ] 
[ ] 
[ ] 


(b) If you are using wxMaxima, you can also input a matrix using an 
on-screen form. Use this form to assign the matrix 


5 —3 
=I =2 
0 4 


to the variable Q as follows. 


(If you are not using wxMaxima, then enter this matrix using the 
matrix command, as in part (a).) 


(i) From the Algebra menu select Enter Matrix... 


Make sure you select Enter Matrix..., not Generate Matrix.... 


A window opens requesting properties of the matrix, as shown below. 


Matrix 


8 Eigenvalues and eigenvectors 
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(ii) The matrix given above has size 3 x 2, so enter 3 as the number 
of rows and 2 as the number of columns, replacing the default 
values shown. 


It is also possible to give information on the type of matrix, but for 


now leave this set to general. 


Finally, enter Q as the variable name to which the matrix will be 
assigned, and click ‘OK’. 


This opens a second window requesting the elements of the matrix. 


Enter matrix 


(iii) Complete entry of the matrix by filling in the elements in this 
second window. When you have finished, click ‘OK’. 


This has the effect of including a matrix input line in your worksheet, 


followed by the appropriate output. 


1 -2 -1 
(c) Assign the matrix 4 T7  0| to the variable R using an 
—5 1 3 


appropriate method. 


'The typed command for entering a matrix is summarised below. 


Entering matrices 


Specify a matrix A:matrix([1,2],[3,4]) ; 
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Once a matrix has been entered, you can use it in calculations, subject to 
the usual laws of matrix arithmetic. The Maxima syntax to use for matrix 
operations is summarised in the following box. You can practise using 
some of these in Computer activity 50. 


Matrix operations 


Operation Syntax Example 
Matrix addition + A+B; 
Matrix subtraction - A-B; 
Scalar multiplication * 2xA; 
Matrix multiplication : A.B; 
Matrix powers A^ A^^3; 
Determinant (of a square matrix) determinant (matrix) | determinant (A); 
Matrix inverse (of a square matrix) | invert (matrix) invert (A) ; 
or matrix’ (-1) AA 
Warnings 


If a matrix operation cannot be performed, for example, if you 
attempt to add two matrices of different sizes, an error message is 
given. 


Remember to use . (not *) for matrix multiplication and ^^ (not ^) 
to find the power of a matrix. 


If you use the symbol * when multiplying two matrices, or the 
symbol ^ to find a power of a matrix, then Maxima may return an 
answer, but in general it will be the wrong answer. (The symbol * 
used between two equal-sized matrices produces a matrix of the same 
size whose elements are the products of corresponding elements of the 
given matrices. The symbol ^ produces a matrix whose elements are 
the elements of the original matrix raised to the given power. These 
are sometimes called element-wise operations.) 


Computer activity 50 Using matrix operations 


In this activity, 


2 -3 3 Al D 9 =] 
P=|-2 1 2], Q=|-1 -2], adR=| 4 7 Of], 
1 —1 4 0 4 -5 1 d 


which are the matrices you input in Computer activity 49. 
(a) Find P +R. 
(b) Find 3P. 


This is a scalar multiple of a matrix, so use the * operator. 
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(c) Find RQ. 


This is matrix multiplication, so use the . operator. 


(d) Find P?. 


Since P is a matrix use the matrix power operator ^^. 


(e) Find the determinant of P. 
(f) Find the inverse of R. 


0 0 
(g) By multiplying P by | 1 |, show that | 1 | is an eigenvector of P and 
Ji 


find the corresponding eigenvalue. 


8.2 Characteristic equations, eigenvalues and 
eigenvectors 
You have seen in Unit 11 that the characteristic equation of an n x n 


matrix A is a polynomial equation of degree n whose solutions are the 
eigenvalues of A. For example, you saw in Example 2 of Unit 11 that the 
characteristic equation of the matrix 6 J is A? — 4A +3 = 0. The 
solutions of this equation are 1 and 3 which are the eigenvalues of the 
matrix. 


You have also seen that the characteristic equation of an n x n matrix A is 
given by 


det(A — AI) = 0, 
where I is the n x n identity matrix. 


The left-hand side of this equation, det(A — AI), is known as the 
characteristic polynomial of the matrix A. For example, the 
characteristic polynomial of the 2 x 2 matrix above is 


2 1 1 0 
«(t 9-9) 
This can be re-written as 


2—A 1 
aet ( 1 p 


ACIDS. 
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The Maxima command for finding the characteristic polynomial of a 
matrix is charpoly. It is demonstrated in the following activity. Once you 
know the characteristic polynomial of a matrix, it is straightforward to 
write down the characteristic equation. 


Computer activity 51 Finding the characteristic polynomial of a 


matrix 
0 1 1 
(a) Find the characteristic polynomial of the matrix A= |-1 2 1 
—1 -1 4 
as follows. 


(i) Enter the matrix into Maxima, and assign it to the variable A. 
(ii) Calculate the characteristic polynomial of A by entering 


charpoly(A,m) ; 


The first argument of the charpoly command is the matrix whose 
characteristic polynomial is wanted, and the second argument is the 
variable in terms of which the polynomial will be expressed. 


You have seen in Unit 11 that characteristic equations are usually 
expressed in terms of A, the lower-case Greek letter lambda. However, 
in Maxima lambda has a special meaning, so you need to use an 
alternative variable, such as m. 


The polynomial obtained is not given in its simplest form. You can 
simplify it by expanding the brackets using the expand command. 


(iii) Simplify the characteristic polynomial by entering 
expand (7) ; 


and hence write down the characteristic equation of A. 


Remember that % refers to the output of the last evaluated command. 


4 0 0 
(b) Find the characteristic polynomial of the matrix B= |3 1 -3 
3 3 1 


and hence write down its characteristic equation. 
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Once you have found the characteristic equation of a matrix, you could 
find the eigenvalues of the matrix by solving this equation, perhaps using 
the solve command. However, Maxima also has a command 
eigenvalues, which you can use to calculate the eigenvalues of a matrix 
directly. This is demonstrated in the next activity. 


First, note that the multiplicity of an eigenvalue of a matrix is the 
number of times that the eigenvalue occurs as a solution of the 
characteristic equation of the matrix. For example, a 2 x 2 matrix with 
characteristic equation (A — 3)? = 0 has the eigenvalue 3, with 
multiplicity 2. Similarly, a 2 x 2 matrix with characteristic equation 

(A — 5)(A — 7) = 0 has the eigenvalue 5, of multiplicity 1, and the 
eigenvalue 7, of multiplicity 1. Larger matrices can have eigenvalues of 
greater multiplicities. For example, a 5 x 5 matrix with characteristic 
equation (A — 5)?(A — 7)? = 0 has the eigenvalue 5, of multiplicity 2, and 
the eigenvalue 7, of multiplicity 3. 


Computer activity 52 Finding the eigenvalues of a matrix 


(a) Find the eigenvalues of the matrix A entered in Computer activity 51 
by typing 


eigenvalues(A); 


I MÀ 
The result of this command is a list containing two elements, each of 
which is itself a list. 


The first of these two lists gives the (distinct) eigenvalues of the 
matrix, and the second list gives their respective multiplicities. 


In this case, each eigenvalue has multiplicity 1. 


The eigenvalues obtained should be the same as those you found in 
Activity 9 of Unit 11. 


(b) Find the eigenvalues of the matrix B defined in Computer activity 51. 


Internally, the eigenvalues command uses the solve command to find 
the solutions of the characteristic equation formed using charpoly. As you 
may recall from earlier (for example, from Subsection 3.3 of this Guide) 
the solve command cannot always find the solutions of an equation. In 
particular, the solve command cannot always find the solutions of a 
polynomial equation of degree 5 or more. This is because, in general, such 
equations cannot be solved using algebraic manipulations alone. 
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Since the characteristic equation of a n x n matrix is a polynomial 
equation of degree n, this means that the eigenvalues command cannot 
always find the eigenvalues of a matrix larger than 4 x 4. If the 
eigenvalues command cannot find the eigenvalues of a matrix, then an 
error message such as the following is displayed. 


solve is unable to find the roots of 


the characteristic polynomial. 


In such circumstances, it may be possible to find numerical approximations 
to the eigenvalues of a matrix, but this is beyond the scope of this module. 


'To find the eigenvectors of a matrix, as well as the eigenvalues, you can use 
the eigenvectors command as demonstrated in the following activity. 


Computer activity 53 Finding the eigenvectors of a matrix 


(a) Find the eigenvectors of the matrix A entered in Computer activity 51 
by typing 


eigenvectors(A); 


The result of this command is a complicated hierarchy of lists! 


The first element of the main list is 
[[1,2,3]1, [1, 1, 11] 


This is a list of two lists. The first of these lists gives the (distinct) 
eigenvalues of the matrix, and the second list gives their respective 
multiplicities. This is the same as the result of the eigenvalues 
command in Computer activity 52(a). 


'The second element of the main list is 
11 
[ [12,21]. [(1,1,11], [[1,1,21] ] 
(Here, additional spaces have been included to help identify the 
different elements of the lists.) This is a list containing three lists, one 
for each of the three eigenvalues. Each of these three lists contains an 
eigenvector corresponding to the appropriate eigenvalue. 


The eigenvectors obtained should be the same as, or multiples of, the 
ones you verified in Activity 9 of Unit 11. Remember that any 
multiple of an eigenvector corresponding to an eigenvalue is also an 
eigenvector corresponding to that eigenvalue. 
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(b) Find the eigenvectors of the matrix B entered in Computer activity 51. 


'The output is similar to that for the matrix A, except that many of 
the lists contain only two elements, since B has only two (distinct) 
eigenvalues. However, there is something else to notice. Consider the 
second element of the main list, which is 


[ PPOs T, 11]. [[1,0,1]1, [0, 1, -11] ] 


This is a list of two lists, one for each of the two eigenvalues, in order. 
The first of the two lists contains an eigenvector corresponding to the 
eigenvalue —2 (which has multiplicity 1), but the second list contains 
two eigenvectors corresponding to the eigenvalue 4, which has 
multiplicity 2. These two eigenvectors are not scalar multiples of each 
other. 


In general, an eigenvalue with multiplicity 2 can have up to 2 
corresponding eigenvectors that are not scalar multiples of each other. 


More generally, an eigenvalue with multiplicity n can have up to n 
corresponding eigenvectors that together have a property called linear 
independence. You may learn about this concept in higher level 
modules. 


The Maxima commands for working with characteristic polynomials, 
eigenvalues and eigenvectors are summarised below. 


Characteristic polynomials, eigenvalues and eigenvectors 


Operation Command Example 


Find the characteristic polynomial | charpoly (matrix , variable) | charpoly(A,m); 
of a matrix 

Find the eigenvalues of a matrix eigenvalues (matrix) eigenvalues (A); 

Find the eigenvectors of a matrix eigenvectors(matrix) eigenvectors(A) ; 


If you began studying this subsection from Activity 10 of Unit 11, now 
return to the unit to continue your study. 


9 Solving linear recurrence systems 


In this section you will learn how to use Maxima to solve 
constant-coefficient linear recurrence systems. First, let’s revise how to 
input and use recurrence systems for sequences in Maxima. You first met 
the commands for this in MST124 Unit 10. 
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Computer activity 54 Working with recurrence systems and 
sequences 


(a) Consider the first-order recurrence system 
$55, du 914-9 me2,3:,..) 
which you met in Example 12 of Unit 12. 


(i) Define this recurrence system in Maxima by entering the two 
commands below. 


x[1]:5; 
x[n] :=2x*x[n-1]+3; 


Notice that you assign the value of the first term, x [1], using the : 
operator that is used to assign values to variables, and you specify the 
recurrence relation using the := operator that is used to define 
functions. 


When you enter a recurrence system in Maxima, you do not need to 
enter a range of values for the index variable, as Maxima assumes a 
suitable range. Maxima can calculate terms in a properly specified 
sequence for integer values of the index variable greater than or equal 
to the values of the index variables of the initial terms. 


If you make a mistake when defining x, you need to remove x from 
Maxima's memory using kill(x); before entering the correct 
definition. There are more details on this in the box following this 
activity. 


(ii) Find the second term of the sequence defined by the recurrence 
system above by entering 


31215 
(iii) Find the tenth term of the sequence. 
(b) (i) Define the second-order recurrence system 
ug =4 uie). un = hai bün- (m2,3:4,...) 


in Maxima using commands similar to those used in part (a). 
(You met this recurrence system in Example 16 of Unit 12.) 


(ii) Find the fourth and tenth terms of the sequence defined by this 
recurrence system. 
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Warning 


When you are using recurrence sequences in Maxima, the system 
remembers all the values of the sequence that it has previously 
calculated, and bases subsequent calculations on them. This improves 
the efficiency of the calculations. However, it also means that 
recurrence systems cannot be changed easily. 


'To change the definition of a recurrence system, you first need to 
remove the system from Maxima's memory using the kill command, 


and then redefine it. 


The commands for defining and working with recurrence sequences are 
summarised as follows. 


Recurrence sequences 


Operation 


Define a recurrence 
sequence 


Calculate a term in a 
recurrence sequence 


Command 


sequence [initial term number]: initial value 
sequence [n] :- expression containing 

previous terms 
sequence [term number] 


Now let's consider how to solve a recurrence system in Maxima, that is, 
how to find a closed form for the sequence defined by the recurrence 
system. The Maxima command for solving a recurrence system is 
solve.rec. This command is not part of the basic Maxima collection of 
commands, so before using it you have to load an additional Maxima 
package, the solve.rec package. You can do this using the 1oad command. 
The package needs to be loaded only once in each Maxima session. 


'The use of the command is demonstrated in the following activity. 
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Example 


x[1]:1; 
x[n] :=2*x [n-1] +x [n-2] ; 


x[50] ; 


9 Solving linear recurrence systems 


Computer activity 55 Solving recurrence systems 


(a) First load the solve.rec package by entering 


load(solve.rec); 


Note the underscore character (_) in the name of this package. 


'The output of this command is a line similar to 
(401) C:/PROGRA~1/MAXIMA~1.0/share/maxima/5.30.0/share/ 
solve_rec/solve_rec.mac 


which indicates the location of the package on your computer. 


(b) Solve the recurrence system 
=b; Ii=2%n1+F3 (2 = 2,340) 
from Computer activity 54(a) as follows. 
(i) First, remove any definition of x from the system by entering 


kill (x) ; 


If you did Computer activity 54 in the same Maxima session as this 
activity, then Maxima will already have values stored for various 
terms in the sequence (xn). This will cause an error in the next part 
of this activity when you try to find a closed form for £n. To avoid 
this, you must first remove x from Maxima's memory using the kill 
command. 


(ii) Solve the recurrence system by entering 


solve_rec(x[n]=2*x[n-1]+3, x[n], x[1]=5); 


The first argument of the solve_rec command is the recurrence 
relation for the sequence, entered using the syntax that you saw in 
Computer activity 54, with one difference, as follows. 


(=, which is used to specify an equation) in the recurrence relation, 


system in Computer activity 54. 


of the sequence that you want to solve for. 


The third argument is the value of the initial term. If you omit this, 
then the closed form is given in terms of an arbitrary constant 4k. 


The closed form obtained should be the same as the one found in 
Example 12 of Unit 12. 


Since you want to solve the recurrence system, you use an equals sign 


not the colon-equals sign (:=) that you used to define the recurrence 


The second argument of the solve_rec command is the general term 
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(c) Solve the recurrence system 


ug =4u =9; Un = Sty =6Un a (2 = 2,3,4,...) 

from Computer activity 54(b) by entering the following two commands. 
kill(u); 
solve rec(u[n]-5*u[n-1]-6*u[n-2], u[n], u[0]24, u[1129); 


Since this is a second-order recurrence system, two initial terms are 
needed. You include both of these as the final arguments of the 
solve_rec command. 


The closed form obtained should be the same as the one found in 
Example 16 of Unit 12. 


'The solve.rec command for solving a recurrence system is summarised 
below. You can practise using it in Computer activity 56. 


Solving recurrence systems 


Solve a recurrence solve rec(recurrence relation, solve rec(a[n]-2*a[n-1], 
system general term, initial term(s)) aln] ,8[1]21) ; 


Note: the solve_rec package needs to be loaded before this command 
can be used. 


Computer activity 56 Solving more recurrence systems 


Solve the following recurrence systems. 


(a) 


gii, du 2 =) WH 25,446.) 


(b) y =12, yn =—3¥n-1t+1 (n=2,3,4,...) 


to = 244 =T Un = Tuüsi-—1249. (n—2,9,4,...) 


(d) Q2 0T Un = Cigar — Vin 2 (Nn = 2,3,4,...) 


If you began studying this subsection from Activity 31 of Unit 12, now 


ret 


94 


urn to the unit to continue your study. 


Maxima accessibility guide 


This section outlines issues associated with the use of Maxima if you 
interact with a computer using only a keyboard, use a screenreader or 
require different colour settings. The information should help you decide 
whether the wxMaxima interface is suitable for you, or whether you should 
consider using an alternative interface, such as the command-line interface. 
It also explains how to configure Maxima to make it easier for you to use 
and contains alternative versions of Subsection 1.2 and Section 2 
describing how to use the command-line interface. 


Using wxMaxima with a keyboard alone 


'The majority of the wxMaxima interface is accessible using a keyboard 
alone. On a Windows computer you can access the wxMaxima menus by 
pressing and releasing the key, navigating through the menus using the 


keyboard arrow keys, and then pressing to make your selection. 
Alternatively, you can press and hold down the [Ait] key, and then press the 


keyboard key that corresponds to the underlined letter in a menu entry to 
select it. In both cases, you can leave menu navigation without making a 
selection by pressing the key. 


On an Apple Mac you can access the wxMaxima menus by pressing and 


releasing the key while holding down the and keys (or just 


on some keyboards), navigating through the menus using the keyboard 


arrow keys, and then pressing to make your selection. You can 
leave menu navigation without making a selection by pressing the key. 


'The wxMaxima toolbar cannot be accessed using the keyboard alone, but 
this is not used extensively in this Guide and alternative typed commands 
for accessing the functions are given as required. 


Changing wxMaxima colours and fonts 


You can change the appearance of the text used by wxMaxima as follows. 
First, click on the following toolbar icon. 


"Uv 


P 


Alternatively, select Configure from the Edit menu (or, on an Apple Mac, 
select Preferences from the wxMaxima menu). 


In the configuration window that opens, select the following ‘Style’ icon. 
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You can then change the colour, style, and sometimes the font and size, of 
individual elements of the wxMaxima interface by first choosing the 
appropriate element from the ‘Styles’ list, and then selecting your 
preferred settings. 


Where a selected element has a greyed out font button this is usually 
because it shares its font with other elements of the interface. For 
example, those elements that form part of the mathematical output of the 
interface share a font that you can change by first clicking on the Math 
font ‘choose font’ button and then selecting your preferred font and size. 
Most of the remaining elements with a greyed out font button share a font 
that you can change in a similar way by clicking on the Default font 
‘choose font’ button. 


Unfortunately the Configuration window is not keyboard accessible, and so 
if you only use a keyboard you may need the assistance of a non-medical 
helper to initially configure the system. If you choose a large font size, 
some horizontal scrolling may be necessary to read all the text displayed in 
a worksheet. 


The Configuration window does not change the colours and fonts used in 
graphs produced by Maxima. The ‘Changing graph properties’ subsection 
below explains how to change the colours, fonts and line thickness used in 
graphs. You should read this subsection in conjunction with Subsection 3.2 
of this Guide, where plotting graphs is revised. 


Using a screen magnifier or screenreader 


The text in wxMaxima can become pixelated when you use a screen 
magnifier. An alternative is to increase the size of the text font as 
described above and/or use the zooming facilities on the Edit menu, which 
you can also access using Alt-I (zoom in) and A1t-0 (zoom out). 


'The wxMaxima interface is not accessible to a screenreader. If you use a 
screenreader you might wish to consider using the command-line interface 
to Maxima instead of the wxMaxima interface. You should ensure your 
screenreader is set to read all the punctuation. Details of using Maxima 
with the command-line interface are given in the section beginning on 
page 98. The graphs produced by Maxima are not accessible to a 
screenreader. You might need assistance when using these. 
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Changing graph properties 


You can change the thickness of curves plotted in graphs by including an 
additional argument in the wxplot2d command, of the form 


[style, [lines, 3]] 


where the number indicates how thick the line should be. By default, lines 
have thickness 1. 


So, to plot a graph of the equation y = 2x with a thicker line, enter, for 
example, 


wxplot2d(2*x, [x, 0, 1], [style, [lines, 311); 


Typing this additional argument each time you want to plot a graph can 
be cumbersome, so you might like to assign this argument to a variable 
that you can reuse. For example, to assign your preferred style to the 
variable thick enter 


thick: [style, [lines, 31]; 
'Then you can plot graphs using commands such as 
wxplot2d(2*x, [x, 0, 1], thick); 


To change the colour of a curve, you can include an argument of the form 
[color, red] as described in Subsection 3.2. Alternatively, you can add a 
further element to the styles argument described above. For example, 
including 


[style, [lines, 3, 5]] 


produces curves with thickness 3 and colour number 5, which is black. The 
colour number codes used by the style argument are given in the table 
below. 


style argument colour numbers 


Colour 
blue 

red 
green 
magenta 
black 
cyan 


Number 


Note that the value of any variable you define to help set curve properties 
(such as thick above) will be lost when you close Maxima at the end of a 
session unless you save your session. The style argument also only 
changes properties of the curve plotted, not other aspects of the graph 
such as the axes and labels. 


Changing the axes and labels, and ensuring such changes are automatically 
used each time Maxima is started is more complicated. A computer file is 
provided in the Accessibility section of the OU Maxima website 
learn1.open.ac.uk/site/maxima which you can download and configure 
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to set the graph properties to your needs. These settings will then be 
automatically used for each graph in every Maxima session. Details of how 
to download and use this are given in the Accessibility section of the OU 
Maxima website. 


Starting, testing and configuring Maxima using the 
command-line interface 
This section describes how to start, test and configure Maxima when using 


the command-line interface. If you are using this interface you should read 
this section instead of Subsection 1.2 of this Gwide, which this replaces. 


Activities in this section are numbered to match the corresponding activity 
in Subsection 1.2. The format of the solutions given at the end of this 
Guide may not exactly match the output you obtain using the 
command-line interface. 


Computer activity 2 Starting Maxima 


Start the command-line interface to Maxima on your computer and keep it 
open to work with as you study this section. 


If you cannot remember how to do this on your computer, see the 
OU Maxima website: learni.open.ac.uk/site/maxima. 


The command-line interface to Maxima is illustrated in Figure 12. 


E Command line Maxima 


Maxima 5.30.0 http://maxima.sourceforge.net 

using Lisp GNU Common Lisp (GCL) GCL 2.6.8 (a.k.a. GCL) 
Distributed under the GNU Public License. See the file COPYING. 
Dedicated to the memory of William Schelter. 

The function bug report() provides bug reporting information. 
(%i1) 


Figure 12 The initial command-line interface window 

At the top of the window is some introductory text, followed by 

(hil) 

This is how Maxima labels lines. The i stands for ‘input’: this is input line 


number 1 and is where you can type your first command. These commands 
either instruct Maxima to perform calculations or change system settings. 
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The percentage symbol (4) is used by Maxima to indicate objects built 
into the system. You will learn more about this later. 


The following activity shows you how to use Maxima to perform a simple 
calculation and checks that your system is working properly. Follow the 
steps using Maxima as you read the activity. 


If you encounter unexpected problems when working through this activity, 
check the Frequently asked questions (FAQ) section on the OU Maxima 
website: learni.open.ac.uk/site/maxima. 


Computer activity 3 Your first Maxima calculation 


Follow these steps to use Maxima to calculate 2 4- 3 x $. 


1. Click anywhere on the Maxima window, to ensure that the text you 
type next is directed to the correct place. 


2. Type the following exactly as it appears here: 
2+3*4/5; 


including the semicolon (;) at the end. 


The text you type will appear next to the input line number (411), as 


indicated by a flashing horizontal line known as the editing cursor. 


3. Press , 


If you forget to type the semicolon (;) at the end of the command, 
the editing cursor simply moves to the next line, waiting for you to 
type more. Maxima will continue to do this until you indicate you 
have reached the end of the command by typing a semicolon. So if 


you forgot the semicolon, then type it now and press again. 


Pressing instructs Maxima to calculate the expression. This 
is known as evaluating the expression. 


The result is displayed under the input line, preceded by the label 
(401). The o stands for ‘output’, so this label identifies the line as 
output line 1. It gives the output corresponding to input line 1. 
Throughout each Maxima session the line numbering starts at 1 and 
then increases by one for each new input (and output). 
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Maxima then displays another input line number, ready for you to 
enter your next command. 


E Command line Maxima 


Maxima 5.30.0 http://maxima.sourceforge.net 
using Lisp GNU Common Lisp (GCL) GCL 2.6.8 (a.k.a. GCL) 
Distributed under the GNU Public License. See the file COPYING. 
Dedicated to the memory of William Schelter. 
The function bug report() provides bug reporting information. 
($11) 243*4/5; 

22 
($01) == 

5 
($12) 


If, in Computer activity 3, Maxima behaved as described, then you have a 
working system to use throughout MST125. 


The result of the calculation in Computer activity 3 was displayed as 
formatted mathematics over several lines, replicating how mathematics is 
usually written. This can, however, cause problems if you are using a 
screenreader. The display behaviour can be changed for every Maxima 
session as described in the following activity. 


Computer activity 4 Configuring the Maxima display of 
mathematics 


You should only do this activity if you wish to change how Maxima 
displays mathematics. 
(a) First exit Maxima by entering the command 
quit(); 
(on some installations you may also have to close the window). 


(b) Turn off the two-dimensional display of mathematics by adding the 
Maxima command display2d:false; to your maxima-init.mac 
configuration file, as follows. 


(i) Follow the instructions in the Accessibility section of the OU 
Maxima website to locate your maxima-init.mac configuration 
file, or create one in the correct location as appropriate. 
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(ii) Open your maxima-init.mac file in a plain text editor (such as 
Notepad on a Windows computer, or TextEdit on an Apple Mac), 
and add the text 


display2d:false; 


on a new line. There are further details on how to do this in the 
Accessibility section of the OU Maxima website. 


(iii) Save your modified maxima-init.mac, and close your text editor. 


(c) Restart Maxima and enter the calculation from Computer activity 3 
again, that is, enter 


243*4/5; 


Notice that Maxima now displays the result on a single line: 


22/5 


Closing wxMaxima 


When you have finished using Maxima, close it by using one of the 
following methods. 


e = Enter quit; 


e Ona Windows computer, click the usual small cross button at the top 
right-hand corner of the Maxima window. 


Getting started with Maxima using the 
command-line interface 


This section revises how to use Maxima with the command-line interface. 
If you are using this interface you should read this section instead of 
Section 2 of this Guide, which this replaces. Differences between the 
wxMaxima interface and the command-line interface in later sections are 
minor and are mentioned within those sections when necessary. 


Activities in this section are numbered to match the corresponding activity 
in Section 2. The format of the solutions given at the end of this Guide 
may not exactly match the output you obtain using the command-line 
interface. 


If you encounter unexpected problems when working through this section, 
remember to check the Frequently asked questions (FAQ) section on the 
OU Maxima website: learni.open.ac.uk/site/maxima. 


Maxima accessibility guide 


101 


Maxima accessibility guide 


Calculating with Maxima 


The way in which a calculation or a command has to be entered in 
Maxima is known as its syntax. The Maxima syntax for the basic 
mathematical operations and other functions is summarised below. You 
can practice using some of these commands in Computer activity 5. 


Mathematical operations and functions 


Operation or function Syntax Example 
Addition t 443; 
Subtraction - 4-3; 
Multiplication * 4*3; 
Division / 4/3; 
Brackets ( and ) 2* (3*4) ; 
Powers, for example 2? ^ 2^3; 

OT ** 2**3; 
Square root, for example V5 sqrt( ) sqrt(5); 
Exponential, for example e? | %e” 4e^3; 

or exp( ) | exp(3); 
Natural logarithm, In log( ) log (8) ; 
Magnitude (absolute value) | abs (F) abs(=3) ; 
sin sin( ) sin(1); 
COS cos( ) cos (3*%pi/2) ; 
tan tan( ) tan(Api/4); 
cosec csc( ) csc(2) ; 
sec sec( ) sec(%pi/3) ; 
cot cot( ) cot (/pi/4); 
sin”? asin( ) asin(sqrt(3)/2); 
cos”! acos( ) acos(1/sqrt(2)); 
tan atan( ) atan(1/2); 


In this Guide, we will use ^ for powers. Many Maxima commands take the 
form of a command name, such as sqrt or abs, followed by a pair of round 
brackets containing one or more objects, such as numbers, that the 
command operates on. Such an object is called an argument of the 
command. In the summary tables in this Guide, arguments of commands 
are often represented by | , as illustrated in the previous table. 


Like many computer systems, Maxima assumes that all angles are 
measured in radians. The names of the inverse trigonometric functions 
sin^!, cos”! and tan”! in Maxima are asin, acos and atan respectively. 
This notation is used by many computer systems; it is short for arcsin, 
arccos and arctan, the alternative names often used for sin^!, cos”? and 
tan-l. 
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The symbols 4e and %pi in the table above are the Maxima 
representations of the mathematical constants e and 7. The % symbol 
indicates that this is the name of a special quantity built into Maxima. 


The constants e and 7 


€ ^e 4e^2; 
T Api 2X pi; 


You should end all Maxima commands with either a semicolon (;) or a 
dollar sign ($). Ending a command with a dollar sign tells Maxima to 
perform the instruction, but not to display the answer. 


Maxima generally ignores all spaces that you type while entering a 
calculation or command, so you can use spaces to make a command easier 
to read. You can also enter line breaks, using or Shift-Enter, 
which are similarly ignored by Maxima. 


Warning: Do not omit multiplication signs! 


A common error when first using Maxima is to omit multiplication 
signs. For example, you might type 


e  2sqrt(2) rather than 2*sqrt (2), or 
e 2(3+4) rather than 2* (344). 


Doing this will result in an error. The displayed error message might 
include one of the following statements: 


e  parser: incomplete number; missing exponent? 


e incorrect syntax: Syntax error 
incorrect syntax: Too many )'s 
incorrect syntax: Premature termination of input at ;. 


e  ... is not an infix operator 


If you receive one of these error messages, first check for missing 
multiplication signs! 


Warning: Take care with question marks (?) 


Do not type question marks when asking Maxima to evaluate an 
expression. The symbol ? is a special symbol in Maxima which 
provides access to underlying systems. 


The one exception to this is when you are using the commands for 
obtaining help from the system, which you will see later. 
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You can practice using Maxima for mathematical calculations in the 
following activity. 


Computer activity 5 Calculating with Maxima 


Calculate the following using Maxima. 


(a) V32 


Since /32 = 323 this can be found by entering 
32^(175)5 
(The brackets in this expression are required since Maxima calculates 


expressions using the BIDMAS rules for the correct order of 
evaluation.) 


Remember to press { Enter )to evaluate the command. 
P 


(b) sin (5) 


Remember that to enter the constant 7 in Maxima, you type Api. 


Maxima returns an exact answer, in the form of a surd. 


(c) tan 302 


Remember to convert the angle to radians, by multiplying 30? by em 


Remember, the Maxima symbol for e is ^e. 


Since there is no simpler exact way of expressing e^, Maxima leaves 
this expression unevaluated. 


(e) 1n3.4 


Remember that the Maxima command for a natural logarithm is 
log( ). Maxima has no command 1n( ). 


Here, since the argument is a decimal number, Maxima returns a 
decimal approximation to the answer. 
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In Computer activity 5 you saw that Maxima usually tries to calculate 
quantities exactly rather than approximating them by decimal numbers. It 
does this by manipulating the symbols in the expression according to the 
rules of mathematics. This is sometimes known as symbolic 
computation. If an exact answer cannot be obtained, the expression will 
be returned unevaluated, as in Activity 5(d). The exception to this is if 
the expression to calculate contains a decimal number, such as 

Activity 5(e) when Maxima will work in terms of decimal numbers. 


If you wish to have a result displayed as a decimal number, you can 
instruct Maxima to do so with the float command. 


Converting to decimal numbers 


Convert to a decimal number | f1oat( ) | float(sqrt(2)); 


The name of the command float arises from the fact that the method 
used by computers to store decimal numbers internally using the binary 
digits 0 and 1 is called a floating-point representation. The float 
command instructs Maxima to convert a number stored symbolically as a 
mathematical expression to one stored using a floating-point 
representation. 


Maxima usually displays decimal numbers to 16 significant figures (though 
it suppresses trailing zeros at the end of the decimal part of a number; for 
example 0.125 is displayed simply as 0.125 rather than as 
0.1250000000000000). This is also the precision to which Maxima 
performs decimal calculations. You will see later how to change the 
number of significant figures displayed. 


Computer activity 6 Simplifying and approximating numerical 
expressions 


For each of the following expressions, use Maxima to simplify it and then 
find a decimal approximation for it. 


(a) 325 (b) 5200 


The expression in part (b) doesn't involve decimal numbers and so 
the result is evaluated exactly. If the result is too long to fit on one 
line, Maxima ends the first line with \ (to indicate that the result 
continues) before continuing to display the result on the following line. 


The decimal form of the result 6 .2230152778611415E+139, is given 
in the command-line interface's way of displaying scientific notation. 
It means 6.2230152778611415 x 101%, 
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Troubleshooting Maxima 


If you are using Maxima and it does not seem to be responding, try the 
following suggestions. 


Troubleshooting Maxima 


If Maxima does not seem to be responding, the following might be 
helpful. 


e Check you have ended your previous command with either a 
semicolon (;) or a dollar ($). Maxima will not do anything 
without them! 


e Check that Maxima is not waiting for you to type something. Is 
a question being displayed that you have not yet answered? 


e If Maxima is taking too long to evaluate a command, you can 
interrupt it by typing Ctrl-C. 


If you do this, you may see the following output, or something 
similar, which you can safely ignore. 


Maxima encountered a Lisp error: 
Console interrupt. 
Automatically continuing. 
To enable the Lisp debugger set *debugger-hook* to nil. 


Note that pressing Ctr1-C at other times may close your Maxima 
session. 


If you encounter other problems, check the Frequently asked questions on 
the OU Maxima website: learni.open.ac.uk/site/maxima. 


Reusing and editing commands 


In Computer activity 6, you may have entered each expression twice, once 
to evaluate it exactly and once for a decimal answer. This is, however, not 
necessary. When you have entered and evaluated an expression, you can 
edit it and then evaluate it again, as shown in the following activity. 
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Computer activity 7 Editing a command 


(a) Calculate 2/3 + 54/27 by entering the following line into Maxima. 
2*sqrt(3)*5*sqrt (27) ; 


(b) Edit your entry in part (a) to calculate 2/3 — 5/27 as follows. 


1. Press the keyboard up arrow key once, so that the previously 
entered command appears against the current input line number 
(on an Apple Mac you may need to hold down the key while 
pressing the keyboard up arrow key once). 


You can access all the commands you have previously entered in your 


Maxima session using the up and down keyboard arrow keys. 


2. Edit the expression by moving the editing cursor using the left 
and right keyboard arrow keys, typing new text and using the 


(backspace) or key to make deletions. 
3. Press : 


Maxima also permits the result of one calculation to be used in another 
calculation. The symbol % represents the result of the last evaluated 
command. So, for example, if you enter sqrt (4) ; Maxima calculates the 
square root of the last evaluated output. You can also use the output (or 
input) of other cells, by typing the cell line label in a calculation. For 
example, if you enter %0573; Maxima calculates the cube of the expression 
in output line 5. 


Computer activity 8 Using the results of previous calculations 


(a) Enter sin(5*%pi/6) ; into Maxima and evaluate it. 


(b) On the following line type 2*%; and enter it. 


Maxima returns 1, which is twice the previous answer. 


(c) Enter 4o| ^2; with F replaced by the output line number 
corresponding to the result of sin(5*%pi/6); entered earlier. 


'This command finds the square of the result of the line referenced. 
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Assigning values to variables 


Maxima allows you to assign a value to a variable, and then use the 
variable in further calculations 


For example, you can assign the value 23 to the variable a, and then 
evaluate an expression in a, such as a? — 3a +2. You can also assign 
expressions, both numerical and (as you will see later) algebraic, and 
indeed any other Maxima object, to a variable, and use that variable in 
subsequent commands. 


To assign a value to a variable you use a colon (:). For example, the 
command 


a:23; 


assigns the value 23 to the variable a. The colon can be read as ‘is 
assigned the value’. 


Warning 


You cannot use = to assign a value to a variable. This symbol is used 
for equations, as you will see later. 


Variable names can be any combination of letters and numbers that begins 
with a letter. For example, Maxima will accept any of the following 
variable names. 


a A solution solution2 x2b 


Variable names are case-sensitive; for example, the variable x is different 
from the variable X. 


Once a variable has been assigned a value, Maxima will remember this 
value for the rest of your session. This may sometimes surprise you. It is 
easy to forget that earlier in the day you assigned a value to, say, x; you 
may then be confused later on when a calculation involving x gives an 
unexpected result. You can, however, tell Maxima to ‘forget’ about a 
variable to which you have previously assigned a value. You can also ask it 
to list all the variables that you have assigned values to. These commands 
are demonstrated in the following activity. 
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Computer activity 9 Working with variables 


(a) Assign the value of v8 to the variable a by entering a: sqrt(8); 


The value of a is displayed as output, in an equivalent form. 


(b) Enter a; 


This displays the current value of a. 


(c) Assign the value of a2 to b. 


Don't forget to include a multiplication sign when inputting aV2. 


The value of b is simplified when displayed. 


(d) Edit the line where you assigned \/8 to a, to assign the value V7 to a. 
(e) Enter b; to display the value of b. 


Notice that the value of b has not changed. The variable b was 


defined when a had its original value. 


(f) Enter values; to list the names of all the variables that are currently 
assigned values. 


The variable names are given within square brackets, separated by 
commas. 'This is how Maxima displays lists. 


(g) Enter kill(a); to remove the variable a from the system. 


You should obtain the output done. 


(h) Enter a; to display the value of a. 


The variable name is output, as a no longer has a value. 


(i) Enter b; to display the value of b. 


b is still defined. 
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Notice that in Computer activity 9(d) and (e), when the value of a was 
changed, the value of b was not affected, since it was defined when a had 
its original value. If you want to update the value of b using a different 
value of a, then you should enter a new value of a then use the up and 
down keyboard arrow keys to obtain the line defining b and then press 


Enter ) to re-evaluate it. 


Working with variables 


Operation Command Example 
Assign a value to a variable : a:23; 
Display the value of a variable | variable a; 

List all user assigned variables | values values; 
Remove an assigned variable kill( variable) | kill(a); 
Remove all assigned variables | kill (all) kill(all); 


Note: here the placeholder variable represents any variable name. 


You can practise using variables in the following activity. 


Computer activity 10 Using variables 


(a) Define the variable a to have the value 42. 


(b) Define the variable b to be equal to sin? a + cos? a. 


2 


To enter, for example, sin^ a in Maxima, you have to type sin(a)^2 


(c) Find the decimal approximation of the value of b. 


System variables 


'There are some variables built into Maxima whose values affect the 
behaviour of the system. These are called system variables. You met 
one system variable in Computer activity 9; the variable values holds a 
list of the names of all the variables you have defined. 


Another system variable, fpprintprec, specifies the number of significant 
figures of decimal numbers that are displayed. The name fpprintprec is 
an abbreviation of ‘floating point print precision’. 


To change the system behaviour so that, for example, only 4 significant 
figures are displayed, use the command fpprintprec:4; to assign the 
value 4 to the variable fpprintprec. 


You can set the variable fpprintprec to any value between 2 and 16. 
Also, setting it to 0 restores the default behaviour of displaying 16 
significant figures. The value of fpprintprec sets the requested number of 
significant figures, but Maxima does not always exactly meet the request. 
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Computer activity 11 Displaying a specified number of significant 
figures 


Display the decimal approximation of 7 to 8 significant figures by doing 
the following. 

(a) Set the system variable fpprintprec to be 8. 

(b) Display the decimal approximation of 7. 


You can reset the values of fpprintprec and all other system variables to 
their original values by using the reset () command. Try entering 
reset(); followed by float (Api); and check that the value of 7 is 
displayed to 16 significant figures. Note that the reset command also 
resets your Maxima line numbers to start at 1 again. 


Resetting system variables 


Reset all system variables 


Another useful system variable is leftjust. Setting this to be equal to 
true using 


leftjust:true; 


makes Maxima display the result of a command left-justified against the 
output line number. The default value of leftjust is false, which 
displays output centred in the line. 


Note that the list of variables displayed when typing values; does not 
include any system variables. If one of your variables does not appear in 
the list displayed by values;, then you have probably used a variable 
name that is also the name of a system variable and hence changed the 
value of that system variable. In extreme circumstances this may change 
the behaviour of Maxima, which can be restored by resetting all system 
variables as described above. 
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Saving and printing your work 


After working with Maxima you will probably want to save your 
calculations so that you can reuse or read them at a later date. 


When you are using the command-line interface you cannot save the 
session directly, as you do when saving other files. Instead there are two 
stages in saving your work: 


e If you want to keep a record of the contents of the command-line 
window, then you need to save a transcript of the session. You can 
then print this transcript if needed. 


e If you want to save the values of all the variables, so that you can 
continue working with them in a new Maxima session, then you will 
need to save the current state of the system. Note that this does not 
give you access to the commands used to define the variables, so you 
cannot edit them in the new session. 


'These two processes are illustrated in Computer activity 13. 


When saving your work, you need to tell Maxima exactly where to save 
the file containing the session state or your transcript by giving the 
complete pathname of the file, which includes not only the name of the file 
but also the hierarchy of folders in which the file is located. Unless you do 
this, Maxima will try to save the file in the location in which the system is 
installed on your computer. Not only is this not recommended, but as a 
computer user, you may not have sufficient permission to do this. 


'The pathname of a file depends on its location, the type of computer you 
are using, and your computer username. In the following activity, we 
assume that a user with the computer username ‘bill’ wishes to save a 
transcript of his Maxima session in a file called transcript.txt within a 
folder called maxima which he created in Computer activity 1. The maxima 
folder is within Bill's ‘Documents’ area on a recent Windows computer 
(that is, Windows vista, 7 or 8). In this case, the full pathname of his file is 


C:/Users/bill/Documents/maxima/transcript.txt 


(Notice that Maxima separates parts of the pathname with a forwards 
slash (/) rather than the usual Windows backwards slash ().) 


If Bill were using an older Windows computer (such as Windows XP), the 
equivalent pathname would be 


C:/Documents and Settings/bill/Documents/maxima/transcript.txt 
or using an Apple Mac the equivalent is 


/Users/bill/Documents/maxima/transcript.txt 
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Computer activity 13 Saving your work 


(a) Create a transcript of your current session by following these steps. 


(i) Enter 
writefile("C:/Users/bill/Documents/maxima/transcript.txt"); 


replacing the pathname with one appropriate to you and your 
system. 


'The path of the file is given in double-quote marks. 


This tells Maxima to save a copy of all the subsequent input and 
outline lines to the file named transcript.txt in your chosen folder. 


'The output 


Starts dribbling to 
C:/Users/bill/Documents/maxima/transcript.txt 


is given, where the pathname is the one you gave in the command, 
followed by the current date and time (in a probably incorrect time 
zone!). Some computer implementations may produce a different 
output. 


(ii) Enter 


playback() ; 


This command redisplays each of the input and output lines of your 
session, and when used in conjunction with the writefile command 
results in a copy of these lines being stored in your chosen file. 


To playback only a selection of the lines from your session, you can 
use a command of the form 


playback([first line number, last line number ]) ; 


which plays back only those lines within the given range. 


(iii) Enter 
closefile(); 


This stops the recording of the session into the file. 


The output Finished dribbling to 
C:/Users/bill/Documents/maxima/transcript.txt is given. 
Again, some computer implementations may produce a different 
output. 
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(iv) Find the file transcript.txt on your computer and open it 
using a text editor, for example, Microsoft Word. 


The file contains a copy of the input and output lines from your 


session. (When opened in some text editors the file may appear all on 
one long line.) 


(b) Save your current session then reload it, by following these steps. 
(i) Enter 


save("C:/Users/bill/Documents/maxima/mysession.lisp", all); 


replacing the pathname with one appropriate to you and your 
system. 


This tells Maxima to save your current session to the file 
mysession.lisp in your chosen folder. (Lisp is the computer 
programming language used to create Maxima, and the save 
command saves your work using commands from this language.) 


'The path of the file is output. 


(ii) Close Maxima by entering 


quit (O; 


(iii) Restart Maxima in the usual way. 
(iv) Load your previously saved session by entering 
load("C:/Users/bill/Documents/maxima/mysession.lisp"); 


using the same pathname you used in part (b)(i). 


‘|| 


This tells Maxima to restore the state of the system to that saved in 
the specified file. 


The error 


assignment: cannot assign to gf_data(characteristic, 
exponent, reduction, primitive, cardinality, order, 
factors_of_order) 


-- an error. To debug this try: debugmode(true) ; 


might be given, which can be ignored. 
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(v) Enter 
values; 


to display the names of all known variables. 


You should see that all the variables from your previous session are 


known. You might like to check their values are as expected too. 


An alternative way of using the writefile and closefile commands is to 
enter the first at the very start of your session and the second at the end, 
so that everything you do is recorded. 


Saving your work 


Operation Command 

Start writing a transcript writefile("file path"); 
Stop writing a transcript closefileO; 

Playback all the input and playback) ; 


output in a session 
Playback a range of input and | playback([/start line number, ‘end line number ]); 


output 
Save the state of a session save(" file path", all); 
Reload a session load("file path"); 


You can also obtain an image of the Maxima window by taking a 
‘screenshot’. To do this on a Microsoft Windows computer, first click on 
the Maxima window, then press A1t-PrintScreen (that is, hold down the 
keyboard key while pressing the PrintScreen key, which may be 
labelled (_PrtScr_) , (Prt Scrn) or something similar). This stores the image of 
the Maxima window, which you can then paste into a suitable application. 
If you are using Windows Vista, Windows 7 or Windows 8, you might like 
to use the ‘Snipping Tool’, available by typing ‘snip’ into the Start menu. 


To take a screenshot on an Apple Mac computer, press Cmd-Ctrl-Shift-4 
(that is, hold down the command) , (ctrl) and (1) keyboard keys while pressing 
the key). Next press the space bar, move the camera pointer over the 
Maxima window and click. This stores the image of the Maxima window, 
which you can then paste into a suitable application. You may also like to 
use the Grab application from the Utilities folder. 


Getting help 


Maxima includes commands that you can use to obtain help. These are 
demonstrated in the following activity and summarised in the table that 
follows. 
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Computer activity 14 Getting help! 


(a) Enter ? float; to show the help for the float command. 


Note the space between the ? and float, which is needed. This is one 
of the few times in Maxima when spaces matter! 


The help information for the command may include more detail than 
you need! After giving the help information, Maxima displays the 


output true. This means that the command you entered was 
successful. 


(b) Suppose that you cannot remember the Maxima command for square 
root, but you do remember that it was something like sqr. 


Type ?? sqr; to list all the help information titles that contain the 
letter sequence ‘sqr’. 


A list of possible titles is displayed. 


Title number 1 is what you were looking for, so enter 1 to display the 
relevant information. 


If more than one title is displayed in response to a ?? command, then 
Maxima will not continue until you enter one of the following: 


e the appropriate number, to display your chosen information 


e several numbers, separated by spaces, to see help information 
on several topics 


e all to display the information on all topics listed, or 


e none to see no information. 


Help commands 


Operation Command Example 


Get help on a command ? command ? float; 
or describe (command ) describe(float); 
Find help information whose | ?? ?? flo; 
title contains the given text | or describe( , inexact) | describe(flo, inexact); 


Note 1: the space after the ?, which is needed. 


Note 2: the second pair of commands displays a list of the titles of all 
help information pages whose title contains the given text. To obtain 
help on a particular topic listed, type the number corresponding to the 
required title in the list. 
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There is further help with Maxima in the Frequently asked questions 
section of the OU Maxima website, at learni.open.ac.uk/site/maxima. 


If you are reading this section as an alternative to Section 2, you should 
nou continue with Section 3. Note, however, the information below 
regarding the plotting of graphs, which you will need when studying 
Subsection 3.2. 


Plotting graphs using the command-line interface 


When using the command-line interface, you should use the plot2d 
command rather than the wxplot2d command described in Subsection 3.2. 
'These commands work similarly, and have the same arguments. The main 
difference is that plot2d opens a new computer window containing the 
graph. 


On some computer systems, the plot2d command opens two windows, one 
containing the graph and a second containing a command-line interface to 
gnuplot, which is the system used by Maxima to plot graphs. To prevent 
this second window opening each time, enter the following command. 


gnuplot.view args:  "-persist ~s"; 


The Accessibility section of the OU Maxima website contains details of 
how to configure Maxima to automatically prevent the display of this 
second window. 


'The graph window can be closed by pressing the usual cross at the top 
right-hand corner, or by pressing Alt-F4. You will not be able to enter 
any commands into the Maxima command-line interface until the graphing 
window is closed. 


The graph can be printed using the print button on the graph window, or 
saved as an enhanced metafile (.emf) file using the save button. These files 
can be opened with many image viewing programs and can be inserted 
into word-processed documents. 
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Computer methods for CAS Activities 
in Books A-D 


Computer solution to Unit 4, Activity 22 


First, load the implicit.plot package to enable the use of the 
wximplicit_plot command for plotting curves given by equations in 
implicit form. 


(4i1) load(implicit.plot); 
(%01) C:/PROGRA-- 1/MAXIMA~1.0/share/maxima/5.30.0/share/contrib/implicit_plot.lisp 


(The output from this command is the location of the package on your 
computer system, and may differ from that shown here.) 


(a) The curve represented by 2x? — 3xy + 5y? — y — 2 = 0 can be plotted 
as follows. 


(4i2) wximplicit plot(2*x^2-3*x*y45*y^2-y-2-0, [x,-2,2], [y,-2,21, 
[gnuplot_preamble, "set size ratio -1"]); 
2 T T T T T T T 


1.57 7 
1- = 
0.5 r J 
0 bs E | 
-0.5 + J 
-1 F E 
-1.5 F 4 


(4t2) 


-2 1 1 1 i 1 L 1 
-2-1.5-1-0.500.51 152 
X 


(%02) 
Here, axis ranges have been chosen to ensure the whole curve is visible 
and the axes have been set to have equal scales. 


If the curve in your plot is not fully visible, or is too small, adjust the 
ranges of your axes until a reasonable figure is obtained. 


The curve appears to be an ellipse. 
The equation 
223? — 3xy + 5? —y—2=0 


is of the form Az? + Bzy + Cy? + Dz + Ey + F = 0 with A = 2, 
B = —3 and C = 5. This gives 


B? — 4AC = (—3)? — 4 x 2 x 5 = 9 — 40 = —31. 
Since B? — 4AC < 0, and A Æ C, the conic is indeed an ellipse. 
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(b) The curve represented by z? + 5xy — y? + 1 = 0 can be plotted as 
follows. 
(13) wximplicit_plot (x*2+5*x*y-y*%2+1=0, [x,-2,2], ly,-2,2], 
[gnuplot_preamble, "set size ratio -1"]); 
2 


1.5 
i 
0.5 
0 
-0.5 
4 
15 


(413) 


-2 
-2-1.5-1-05005 11.5 2 
X 


(403) 
'The curve appears to be a hyperbola. 


'The equation 
xr? +5ry— y +1=0 


is of the form Az? + Bzy + Cy? + Dz + Ey + F = 0 with A = 1, 
B = 5 and C = —1. This gives 


B? — 4AC =5?° — 4x 1x (= = 25 +4 = 29. 
Since B? — 4AC > 0, the conic is indeed a hyperbola. 


Computer solution to Unit 4, Activity 35 


(a) The parametrisation of the line joining (—2, 4) and (3, 1) found in the 
solution to Activity 33(b) of Unit 4 is 


z=—2+5t, y=4-3t (0<t<1). 


This can be plotted as follows. 


(Aii) wxplot2d([parametric, -2+5*t, 4-3*t, [t,0,1]], 
[x,-3,4], [y,0,5], [xlabel,"x"], [ylabel,"y"]); 
5 


- 
ti) 2L 


(101) 
Here, axis ranges have been chosen to clearly show the ends of the line 
segment, and suitable axis labels added. 
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(b) The parametrisation of the circle with centre (—3, 1) and radius 4 
found in the solution to Activity 34(a) of Unit 4 is 


x = —3 + 4cost, y=1+4sint (0<t< 2r). 


This can be plotted as follows. 


(4i2) wxplot2d([parametric, -3+4*cos(t), 1+4*sin(t), [t,0,2*%pi]], 
[x,-8,2], [y,-4,6], [xlabel,"x"], [ylabel,"y"], 
[gnuplot_preamble, "set size ratio -1"]); 

6 


(2 ^" gL | 


(402) 

Here, axis ranges have been chosen to clearly show the entire curve, 
and suitable axis labels added. Since a circle is being plotted, equally 
scaled axes have been used. 


(c) The parametrisation of the semicircle in Activity 34(b) of Unit 4 
found in the solution to this activity is 


zx —2--2cost, y=1+4+2sint (-5s:«z). 


'This can be plotted as follows. 


(%i3) wxplot2d([parametric, 2+2*cos(t), 1+2*sin(t), [t,-/pi/2,%pi/2]], 
[x,0,6], [y,-2,41, [xlabel,"x"], [ylabel,"y"], 
[gnuplot_preamble, "set size ratio -1"]); 


T T T T T 
3 X "m 
2 L al 
(63) | | 
-1 L B 
29 1 L 1 L 1 
O 1 2 3 4 5 6 
Xx 
(3103) 


Here, axis ranges have been chosen to clearly show the entire curve, 
and suitable axis labels added. Equally scaled axes have been used. 
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Computer solution to Unit 7, Activity 37 


(a) From the ‘Maxima reference guide’, the commands for sinh z and 
cosh z are sinh(x) and cosh(x) respectively. So the graphs of 
sinh x, cosh x and ie" can be plotted as follows. 


(Aii) wxplot2d([sinh(x), cosh(x), exp(x)/2], [x,-2,2]); 


-3 F si 
-4 i ' i ' 
2 15 -1 -0.5 0 0.5 1 1.5 2 
x 
(%01) 
Alternatively, since sinh z = 3 (e* — e^?) and cosh x = 3 (e” + e^), 


the graphs of sinh z, cosh x and ie" can be plotted, with an 
appropriate legend, using 


wxplot2d([(exp(x)-exp(-x))/2, (expGO *exp(-x))/2, exp(x)/21,[x,-2,21, 
[legend,"sinh x","cosh x","exp(x)/2"]); 


(b) When z is positive and of large magnitude, the expression e~” is close 
to 0, and e? is positive with large magnitude. So sinh x and cosh are 
both positive with large magnitude, and both close in value to ie". 

(c) When z is negative and of large magnitude, the expression e* is close 
to 0, and e^* is positive with large magnitude. So sinh x is negative 
with large magnitude, and is close in value to jet , whereas cosh x 


is positive with large magnitude, and is close in value to se". 
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Solutions to Computer activities 


Solution to Computer activity 3 
(hit) 2+3*4/5; 
(A401) > oa 
So 2+3x 555 

Solution to Computer activity 5 


(a) (411) 32^(1/5); 
(401) 2 
So, 1/32 = 2. 

(b) (412) sin(%pi/3) ; 


(402) E 


om v3 
So, sin (5) => 
(c) (413) tan(30*%pi/180) ; 
(408) + 
"V3 
So, tan 30° = e 


(d) (414) her4; 
(%04) het 


(e) (4i5) log(3.4); 
(%05) 1.223775431622116 


So, In 3.4 = 1.223775431622116. 


Do not be surprised if your line numbering differs from ours. Each of 
our solutions is prepared as a separate Maxima session with the line 
numbering starting at 1 (unless it is a continuation of a previous 
activity); your numbering will continue to increase for the duration of 
your session. 
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Solution to Computer activity 6 


(a) (111) sqrt(325); 
(hot) 5/13 
(412) float(sqrt(325)); 
(402) 18.02775637731995 
So V325 = 5/13 ~ 18.02775637731995. 
(b) (13) 5^200; 
(403) 622301527786114170714406405378[80 digits]551224821247160434722900390625 


(414) float (57200); 
(%04) 6.2230152778611415 10139 
So 5200 ~ 6.2230152778611415 x 101%, 


Solution to Computer activity 7 
(a) (Aii) 2*sqrt(3)*5*sqrt (27) ; 
(Y01) 1743 
So 2/3 + 5427 = 1743. 
(b) (412) 2*sqrt(3)-5*sqrt (27) ; 
(%02) -1343 
So 2/3 — 5/27 = —134/3. 
Solution to Computer activity 8 


At the end of this activity, your worksheet should look like the following, 
although line numbers may differ. 


(i4) sin(5*/4pi/6); 
E T 

(404) 2 

(hib) 2X: 

(405) 1 

(416) pum 

(406) 4 


Solution to Computer activity 9 


(a) (411) m: sqrt (6): 
(3,01) 23/2 


(b) (4112) a; 
(3502) 23/2 


(c) (413) b: a*sqrt(2); 
(403) 4 


(d) (414) a: sqrt(7); 
(04) V7 
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(e) (415) b; 
(405) 4 


(f) (416) values; 
(406) [a,b] 


(g) (417) kill(a); 
(WoT) done 


(h) (418) a; 
(408) a 


(i) (419) b; 
(409) 4 


Solution to Computer activity 10 


(a) (hi1) a:42; 
(Y 01) 42 


(b) 02) btain(a)^24eog(a)^2; 
(402) sin(42)?*cos (42)? 


(c) (X13) float(b); 
(403) 1.0 


Solution to Computer activity 11 

(Aii) fpprintprec:8; 

(%01) 8 

(This sets the number of significant figures to be shown in the display of 
decimal numbers.) 


(412) £loat (pi); 

(%02) 3.1415927 

So the value of 7 is 3.1415927 (to 8 s.f.). 
Solution to Computer activity 12 


The circumference of a circle of radius 5 is 31.42 (to 2 d.p.) and the area is 
78.54 (to 2 d.p.). 
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Solution to Computer activity 14 
(a) (11) ? float; 


wna 


-- Function: float (<expr>) 
Converts integers, rational numbers and bigfloats in <expr> to 
floating point numbers. It is also an ‘evflag’, ‘float’ causes 
non-integral rational numbers and bigfloat numbers to be converted 
to floating point. 

There are also some inexact matches for ‘float’. 

Try ‘?? float’ to see them. 


(%01) true 


(h12). TT sor: 
0: isqrt (Functions and Variables for Number Theory) 
1: sqrt (Root Exponential and Logarithmic Functions) 
2: sqrtdenest (Package sqdnst) 
3: sqrtdispflag (Functions and Variables for Display) 
Enter space-separated numbers, ‘all’ or ‘none’: i; 
-- Function: sqrt (<x>) 
The square root of <x>. It is represented internally by 
‘<x>^(1/2)?. See also ‘rootscontract’. 
‘radexpand’ if ‘true’ will cause nth roots of factors of a product 
which are powers of n to be pulled outside of the radical, e.g. 
fsqrt(16*x^2)? will become ‘4*x’ only if ‘radexpand’ is ‘true’. 


(402) true 


Solution to Computer activity 15 
(a) (hi1) p:2*x*sqrt (143*x^6) ; 


(fol) vV3x8+1+2x 


(b) (412) aubat(3,*.p); 


(402) 2 / 547+6 
So, when x = 3, 2x + v1 + 32% = 24547 + 6. 


(c) (413) f(x) :=(x%4+x+1) / (x^3-13*x*12) ; 


x^-x1 


(503) f(x) {= 312 


(d) hi4) £02; 


125 


Solutions to Computer activities 


(e) (4i5) £(3); 
expt: undefined: 0 to a negative exponent. 
#0:  f(x-3) 
-- an error. To debug this try: debugmode(true) ; 


(f) 085 gedd); 
(406) g(1) 


(g) (hi7) functions; 
(%07) [£G0] 


Solution to Computer activity 16 


(a) (411) expand( (1+sin(x))*(x+4)%2 ); 

(%01) x?sin(x)48 x sin(x)+16 sin x) *x?48 x416 

So, 

(1+ sin z)(x + 4)? = x? sin(x) + 8x sin(x) + 16 sin(x) + 27+ 8x +16. 
(b) (412) factor( 3*x^3*11*x^2*8*x-4 ); 

(%02) (x*2)? (3x-1) 

So 319 + 113? + 8x — 4 = (x + 2)? (3x — 1). 
(c) (413) fullratsimp( (2*x^3-34*x^241)/(x-1) J; 

(403) 2x?-x-1 

21? — 33? 1 
So, diet ab UN TS 
g-—1l 

(d) (414) radcan( log(x^3-3*x-2)-2*1og(x*1) ); 

(404) log(x-2) 

So, In(z? — 3a — 2) — 21n(z + 1) = ln(z — 2). 

You might like to think about how you could obtain the same answer 

by hand. 


(e) (415) logcontract( log(2*x)-2*log(x)+log(y) ); 
2y 
(%05) 108(2) 


2 
So, In(2z) — 21n z + In y = ln (2). 
a 


(f) (16) trigexpand( tan(3*x) ); 

(106) 3tan(x)-tan(x)? 
i 1-3 tan(x)? 
3tanz — tan? z 
1 — 3tan?z 

(g) (Liv) trigreduce( cosGO^4 ); 

CLoT) cos (4x) *4cos (2x) *3 
8 


So, cos! x = A 


So, tan(3x) = 
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(h) (418) trigsimp( cos(x)%2*(1+tan(x)%2) ); 


(408) 1 
So, cos? z (1+ tan?(z)) — 1. 


(i) (419) trigrat( tan(x)/sin(2*x) ); 


VA 9 LLL 
(409) cos(2x)+1 
tan x 1 


^ Sin (22) ~ cos(2x) + 1° 


Solution to Computer activity 17 
(a) (1411) wxplot2d(x*sin(x), [x,-5,51); 


2 T T 
1 - 
0r 
(itl) £ 2- 
-3 + 
4b 
-5 
x 
(%01) 


(b) (%i2) f(t) :=(t^2+2*t+3)/(t^2+2*t-3); 


t?+2t+3 
162) Eme 
(402) f(t) apes 


(113) wxplot2d(f(£), [t,0,3]); 
3 T 


000 ; r 
& 2000 } 
& 1000 + 
ai 0 4 
= -1000 + 
(4t3)  -2000 f 
& -3000- 
& -4000 + 
= -5000 + 
-6000 
0 05 1 15 2 
t 
(403) 


Solutions to Computer activities 
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Solutions to Computer activities 


(c) ORE weplot2d(f(00, [6,0,31, Ey,-10,10]D4 
plot2d: some values were clipped. 


10 
e» 
5 
T 
2 
x 0 
(yta) pi 
T = 
` 
-10 
t 
(404) 


Solution to Computer activity 18 


(a) (4i1) p(x) :2x*cos (x); 
(%01) p(x) :=x cos (x) 


(412) q(x) :=x-x%3/2; 


x3 


(402) q(x) x7 


(413) wxplot2d([p(x),q(x)], [x,-3,31); 
15 T T 


x*cos(x) 
10 x-x43/2. — —— 4 


(43) 


(403) 
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Solution to Computer activity 19 


(hi4) wxplot2d([p(x),q(x)], [x,-3,3], [color, black, green], 
[legend, "p^, Mgt] " [ylabel,"y"]); 


15 i T T T T 
10 
5 
(44) 
-10 
-15 1 1 | 1 1 
-3 -2 -1 0 1 2 3 
X 
(404) 


Solution to Computer activity 20 
(a) (Aii) load(implicit_plot) ; 


Solutions to Computer activities 


(ol) C:/PROGRA-- 1/MAXIMA--1.0/share/maxima/5.30.0/share/contrib/implicit plot.lisp 


(The output from this command is the location of the package on your 


computer system, and may differ from that shown here.) 
(x^ 24y at, [x,-2,21, ly, 2,2105 
T T T T 


(b) (412) wximplicit.plot 
2 T 


T 
1.5 H E 4 
1 lE 


05r 


Uim > + 


(%02) 
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Solution to Computer activity 21 
(It is assumed the implicit. plot package has already been loaded.) 


(11) wximplicit.plot(x^2*y^2-1,[x,-2,2], [y,-2,2], [gnuplot.preamble,"set size ratio -1"]); 
2 


eL ee 


Qux + OT 


-2 1 1 1 i 1 1 L 
2 15 -1 -05 0 05 1 15 2 


x 


(hol) 


Solution to Computer activity 22 


(It is assumed the implicit_plot package has already been loaded.) 


(a) 


(11) set_plot_option([gnuplot_preamble, "set size ratio -1"])$ 


(b) (12) wximplicit plot((x-10^249^2-8, [x;-1,3]. Ly;-2,2]1); 
2 


H 


nee c 


Gey + "D 


-2 
1-05 005 1 15 2 25 3 
x 


(402) 
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Solution to Computer activity 23 


(It is assumed the implicit plot package has already been loaded, and 

the plot settings changed to use equal scales for the z- and y-axes.) 

(Lii) wximplicit plot([y-x^2-6*x*9, (x-3)^2*(y-2)^2-3], [x,0,6], [y,0,6]); 
6 


y = x^2-6"x+9 : 
5 L (y12)/2+(x-3)2 = 3 |] 
4 iz 4 
3 E al 
(Ati) 
2 = 4 
1 H E 
0 L Í 1 L 
0 1 2 3 4 5 6 
x 
(401) 


Solution to Computer activity 24 
(a) (hi1) solve( x^3-4*x^242*x*4-0 ); 

CLo1) [x=1-/3,x=/3+1,x=2] 

So, the three solutions are x = 1 — V3, z = 1 + V3 and z = 2. 
(b) (412) eqn:y-log(u*x*v); 

(402) y=108 (u x+v) 


(13) solve(eqn, x); 
hel -v 


(403) [x=———] 
M 


e?" —v 


So, the rearranged equation is y — 
u 


(c) (414) solve(tan(3*x-1)=sqrt(3)); 
solve: using arc-trig functions to get a solution. 
Some solutions will be lost. 
7*3 
(ho4) [x=] 


T+3 
So one solution of the equation is x = — 


(d) (415) solve(x^2-4*x*20-0); 

(05) [x=2-4%i,x=4/i+2] 

The two complex solutions are x = 2 — 4i and x = 2 + 4i. 
(e) (4i6) solve(x^9*2*x-1-0); 

(406) [O-x?*2x-1] 

Maxima was unable to solve this equation. 
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Solution to Computer activity 25 
(a) (Aii) solns:solve(x^2-x-1-0); 


1 
(o1) [x p iX 
So the solutions are z = : E and x — Í 
(b) (412 solns [1]; 
(402) — 


(c) (413) rhs(solns[1]); 


(d) (414) rhs(solns[2]); 


(io) A 


Solution to Computer activity 26 
(a) (i) (il) eq1:3*x+2*y=1; 
(%01) 2y+3x=1 
(/12) eq2:5*x+3*y=3; 
(%02) 3y+5x=3 


(ii) (413) solve( [eqi,eq2], [x,y] ); 
(%03) [[x=3, y=-4]] 
The solution of the equations is z = 3, y = —4. 
(b) (414) solve([y=x"2-6x*x+9, (x-3)^2* (y-2)*2=3] , [x,y]); 
V5-5  v/5+*3 _V5+5 V5-3 


(404) [ i, y SS l], [x= a r M. 


/5-7 = = 3 _V5+7 V5+3 
Bh a Bo 


[x=- —], [x= 


So the solutions (the seta of intersection) are 
5-v5 — 3-vb — 5 E v5 — 3- v5 
= P) Y= 2 T= Y= 2 > 
|7T-vV5 — 3- V5 7 t VE 3-4 V5 


g 47 —3 and z = 2E 
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Solution to Computer activity 27 


(a) 


(i) Chii) f(x) :=x*9+2*x-1; 
(%01) f(x) :=x%+2x-1 


üt) (2412) wxplot2d(f (x), [x,=1,11); 
2 


& 
(112) à 
X 
(402) 


Since f (0) « 0 and £(1)> 0, there must be a solution in the 
interval (0, 1). 
Gii) (13) tind oroot(f(D, x, 0, 10; 
(403) 0.49904018035697 
So a solution of z? + 2x — 1 = 0 is 0.499 (to 3 s.f.). 
(114) g(x) :=cos(x)-sqrt (x) ; 
(h04) g(x) :=cos(x)-./x 
(Aib) wxplot2d(g(x), [x,0,1]); 
1 


0.8 F 
3 0.6 |. 
S 0.4 + 
- 02r 
(t5) = 
[77 0r 
fe) 
9 .02- 
0.4 + 
-0.6 1 į 1 fi 
0 0.2 0.4 0.6 0.8 1 
Xx 
(%05) 


Since g(0)> 0 and g(1)< 0, there must be a solution in the 

interval (0, 1). 

(If you chose an interval other than (0,1) then your graph will look 

different from the one given above. If, in particular, your graph does 
not cross the x-axis, then you will need to plot the graph of g again 

using a larger interval.) 


(18) find root(g(x), x, 0, 10; 
(06) 0.64171437087288 
So a solution of cos z = yz is x = 0.642 (to 3 s.f.). 


Solutions to Computer activities 
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Solution to Computer activity 28 


(a) 


134 


(Hil) diff (cos(x)/(x224+3), x); 


sin (x) 2xcos(x) 
Chol) -— = ur a 
x*+3 (x2+3) 
i 2 
So the derivative of f(x) = 23 is f'(x) = "ES — ast 


(4i2) diff (u*2*sin(u), u, 2); 

(%02) -u? sin(u)+2 sin(u) -4u cos (u) 

So, the second derivative of g(u) = u? sin u is 
g" (u) = —u? sinu + 2sinu + 4u cos u. 


(13) difft(lopg(t)/t^2, t, 395 
26 241og(t) 


(103) —- > 
T t 
Int 26  24]nt 
So, the third derivative of h(t) — > is ns m = . 
(%i4) integrate(x*sin(x)^2, x); 
(oa) PMA teos (00-2 
So, n — pu 4 c. 


(4i5) integrate(t^2*/4e^t, t, O, 1); 
(%05) %e-2 


1 
so, f tel dt =e — 2. 
0 
(416) integrate(1/x^2, x, 1, 2); 
1 
y 6 = 
(406) 5 


2 

1 1 

S — dr = +. 
o f 2341-75 


(417) integrate(tan(Ae^-x), x, 0, 1); 


1 
(LoT) f tan(%e Y) dx 
0 


1 
Maxima cannot evaluate / tan(e”) dx exactly. 
0 


Solutions to Computer activities 


Solution to Computer activity 29 
(a) (hi1) fGO:-x*sqrt(1*4*x^2); 


(fot) f(x) :2xvV 1+4 x? 


(b) (1412) d£ (D :ediffCF GO x): 
(402) df (x) :=diff (f (x) ,x) 


(c) (13) dE (3); 
diff: second argument must be a variable; found 3 
#0: df (x=3) 
-- an error. To debug this try: debugmode(true) ; 


(d) (aa) dfGc) :='' (diff (f (x),x)); 


4x? 
(%04) df (x) :- V 1+4 x2+ 
V1+4 x? 
Az? 
So fs) = Vi F + == 
o f'(x) = v1 + 42? VEEE 


(e) (4i5) df (3); 


73 
%05 Ma 
LE 


So f’(3) = —=. 
(f) (416) F(x) :='' (integrate(f(x),x)); 


: ere 
(406) F(x) o 


(g) (417) F(1)-F(0) ; 
: 53/2 1 
(LoT) 49 12 


(418) integrate(f(x),x,0,1); 
3/2 


; "S 
We) a ae 
: 53/2. ] 
so, | x) dz = F(1) — F(0) = — —- —. 
fle) ae = FQ) - FO) =~ 5 


Solution to Computer activity 30 


(a) (i) (iD kill (a); 
(401) done 


(ii) (412) integrate(x^n, x); 
Is nti zero or nonzero? nj 
n+1 


7 x 
(502) ud 
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Solutions to Computer activities 


(b) (i) (413) assume (notequal(n,-1)); 
(403) [notequal(n,-1)] 


(ii) (4i4) integrate(x^n, x); 
nti 
(%04) 


n+1 


(ii) (115) factsO; 
(405) [notequal (n,-1)] 
So Maxima knows n 4 —1. 

(iv) (4i6) forget (notequal (n,-1)); 
(506) [notequal (n,-1)] 


(v) AXX7) facts0; 
(%07) O 


No facts are now known. 


Solution to Computer activity 31 
(a) (i1) integrate(log(sin(x^3)), x, 0, 1); 
1 
(hot) f log (sin(G?))dx 
0 
Maxima cannot evaluate this integral exactly. It returns the integral 
unevaluated. 
(b) (4i2) quad.qags(1og(sin(x^3)), x, 0, 1); 
(%02) [-3.024256571599236,3.3575992781634104 10°*°, 189, 0] 
(Remember 3.3575992781634104 1071? means 3.3575992781634104 x 10715.) 


1 
So T In(sin(z?)) da = —0.302 (to 3 s.f.). 
0 


Solution to Computer activity 32 
(a) (Aii) quotient (59,8) ; 
(01) T 
The quotient on dividing 59 by 8 is 7. 
(b) (412) remainder (59,8) ; 
(402) 3 
'The remainder on dividing 59 by 8 is 3. 
(c) (413) quotient(3^100, 2^100); 
(403) 406561177535215237 
(i4) remainder(3^100, 2^100); 
(o4) 503611859755855824366132007889 
So 31% = 406561177535215237 x 219? + 503611859755855824366132007889. 
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Solutions to Computer activities 


(d) (415) gcd(108,93); 
(%05) 3 
So the highest common factor (greatest common divisor) of 108 
and 93 is 3. 


(e) (416) gcd(3^100, 2^100*11); 
(%06) 9 


So the highest common factor (greatest common divisor) of 312° and 
2100 + 11 is 9. 


Solution to Computer activity 33 

(a) (Aii) gcdex(108,93); 
(%01) [-6, 7,3] 
So the highest common factor of 108 and 93 is 3, and 
108 x (-6) 4-93 x 7 =3. 


(b) (412) gcdex(3^100, 2^100*11); 
(402) [27354783488210640033683046227 , -11121392986187780374611343715203446360889730414 , 9] 


So the highest common factor of 3!°° and 21% + 11 is 9, and 
3100 x 27354783488210640033683046227 + (2109 + 11) x 
(—11121392986187780374611343715203446360889730414) — 9. 


Solution to Computer activity 34 
(a) (411) wxplot2d( [parametric, t+4, 2*t+3, [t,-1,2]] ); 


2*t-3 


Chel) 


N wo A OC! O N 
T 


(hol) 
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Solutions to Computer activities 


(b) (412) wxplot2d( [parametric, t*4, 2*t+3, [t,-1,2]], 
[ylabel,"y"] ); 


(162) 


(%02) 


O-—nmco co o 


[x,0,8], 


[y,0,8], 


[xlabel,"x"], 


(c) (413) wxplot2d( [parametric, 3*t^2, 6*t, [t,-2,2]], 


(443) 


(%03) 


10 F 


[55595 435]; 


rabei; "z" 


[ylabel,"y"] ); 


[x,-5,15], 
15 


5 L 


o 


Note that axis ranges have been chosen to clearly show the curve 
represented by the given equations and range of values of t. 
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Solution to Computer activity 35 

(a) (i) (hii) x(t) :=12-12*sin(t)"2; 
(%01) x(t) :=12-12 sin(t)? 
(112) y(t) :712*cos (t) ; 
(402) y(t) :=12 cos (t) 


(ii) (A4i3) wxplot2d( [parametric, x(t), y(t), [t,0,%pi]], 
[x,-5,15], [y,-15,15], [xlabel,"x"], [ylabel,"y"] ); 
15 T T 


10 + J 
5 F al 
(t3) | 
-10 + 4 
15 
-5 0 5 10 15 
X 
(403) 


(b) (1414) x(t):-t*cos(t) ; 
(404) x(t) :=t cos(t) 


(%i5) y(t) :=t*sin(t); 

(405) y(t) :=t sin(t) 

(%36) wxplot2d( [parametric, x(t), y(t), [t,0,%piJ], [xlabel,"x"], [ylabel,"y"] J; 
2 T T T 


(4t6) > 0.8 


w 
a 
[^ 
; 
m 
a 
m 
E 


(406) 
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Solution to Computer activity 36 


(a) (i) | (i01 x(t) :=t%3; 
(201) x(t):st* 
(112) y(t) :=sin(10*t) ; 
(402) y (t) :=sin(10 t) 


(413) wxplot2d( [parametric, x(t), y(t), [t,0,5]], [xlabel,"x"], [ylabel,"y"] J; 
1 T T T T T 

0.8 7 

0.6 7 

0.4 4 

0.2 4 

> 0 7 

(43) -0.2 4 

-0.4 4 

-0.6 4 

-0.8 4 


0 20 40 60 80 100 120 140 
(403) 
(111) 


(414) wxplot2d( [parametric, x(t), y(t), [t,0,5]], [nticks,200], 
[xlabel, "x"] g [ylabel, Meth] ) : 


qa ^" go i 


0 20 40 60 80 100 120 140 


(%04) 
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(b) 
(4i5) wxplot2d( [parametric, t*sin(2*t), t*cos(2*t), [t,0,10]], [nticks,100], 


[xlabel, get] : [ylabel, ip] ) A 
10 T T T T T T T T 


(/t5) 


(4,05) 


(16) wxplot2d( [parametric, sin(2*t), cos(3*t), [t,0,20]], [nticks,400], 
[xlabel, "x"] , [ylabel, wpa] ) ; 
, 


(t6) 


-1 -0.8-0.6-0.4-0.2 0 02040608 1 
x 


(406) 


Solution to Computer activity 37 
(a) (1) Qui Ge) Sees; 

(401) x(t) :=3 t? 

(%i2) y(t) :=8*t02; 

(102) y(t) :=8 t? 

(1413) £(x) :=2*x02; 

(%03) f(x) :=2 x? 


(y (jie) pi[parametric, a(t), y(t); Dt,-1,1]]; 
(%04) [parametric, 3t?, 8t?, [t,-1,1]]; 
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Gu) (1415) wxplot2dt [p,f(3)]1, [x,-2,2], [mbicks,50] 2j 


3'1^3, 8'^2 — — 
2x2 /7 


(t5) 


O-—nmo HK ODN o 


(405) 


(b) (416) p: [parametric, sin(t)^3, cos(t)^3, [t,0,2*/^pill; 
(406) [parametric, sin(t)?, cos(t)?, [t,0,27]]; 
(417) q: [parametric, cos(t)^2*sin(t), sin(t)^2*cos(t), [t,0,2*/4pill; 
(%07) [parametric, cos(t)? sin(t), cos(t) sin(t)?, [t,0,27]]; 


(18) wxplot2d([p, ql, [nticks,50]); 
1 


1 1 1 
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 08 1 


(408) 


Solution to Computer activity 38 
(a) 


(Aii) wxplot2d( [parametric, 3+2*sec(t), -5+tan(t), [t,-4pi/2,3*^pi/211l, 
[xlabel A sett] : [ylabel A wyn] ) A 
2e+016 T T T T T T T T T 


1.5e+016 
1e+016 + 
5e+015 - 
0 os 
-5e+015 p 
-1e«016 + 
-1.5e+016 r 


-2e+01 6 L L 1 L ji L L L L 
-1.5e-401666015356e4 063036 29180936 861836016 


X 


(ti) 


(hol) 
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(412) wxplot2d( [parametric, 3+2*sec(t), -5+tan(t), [t,-4pi/2,3*^pi/211l, 
[x,-2,8], [y,-10,0], [nticks,200], [xlabel,"x"], [ylabel,"y"] ); 
0 T 


-2 - 4 
-4 Lr 4 
(t2 C gl | 


(402) 


(%i3) p:[parametric, 3*2*sec(t), -5+tan(t), [t,-4pi/2,3*74pi/211; 
(403) [parametric, 2sec(t)+3, tan(t)-5, is i 


(4i4) wxplot2d( [x/2-13/2, -x/2-7/2, pl, [x,-2,8], [y,-10,0], [nticks,200], 
[xlabel, ingen] : [ylabel, 10] ) n 
0 
x/2-13/2 ——— 


-x/2-7/2 — — | 
2*sec(t)+3, tan(t)-b — — 


(t4) T e 


(%04) 
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(d) (4i5) wxplot2d( [parametric, -1-4*t^2, 4+8*t, [t,-2,2]], 


[xlabel, ty] : [ylabel, m1] ) x 
20 


15 
10 


(t5) 7 o 


(%05) 


(e) 


(16) wxplot2d( [parametric, 1+2*cos(t), 2+3*sin(t), [t,0,2*%pil], [nticks,100], 
[xlabel, usc] , [ylabel, Het] ) ‘ 
T T 


(446) 


(406) 
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Solution to Computer activity 39 
(a) (hi1) partfrac( 7*x/((x*5)*(x-2)), x 5; 
Co) EA 
: xtb x-2 
Ta 5 " 2 
(z-5)r—c2) z45 x-2 
(b) (412) £actor(0D; 
(102) — > 


(x-2)(x*5) 


So, 


Solution to Computer activity 40 
(a) (Aii) partfrac( 9/((2*x*1)*(x-1)), x); 


3 
yi ï po 
et) x-1 2x+1 
9 


08 6 
'(2z--1)(r—-1) z-1 2z-4 
(b) (112) partfrac( 3*x/(x^2+x-2), x); 

(iun cac 

(02) —~+—— 

x+2 x-1 

3x 2 1 
AE ee i 7 
we+a-2 +2 a-1l 
(Note that the denominator, z? + x — 2, of the given expression can be 
factorised as (x + 2)(x — 1).) 


(c) (413) partfrac( (4*x+8)/(4*x-3)%2, x); 
11 


1 
(103) 8 ara? 
da 4-8 1 11 


4-39] ^ 42-3 mF 
(d) (414) partfrac( (8*x-4)/((5b*x-3)*(x-1)^2), x); 


So 


So, 


So 


55 11 6 
y, 4 — + 
(ho4) 5x-3 x-1 (x-1)? 
8r +4 55 11 6 
So, = | 


(5x —3)\(£x— 1)  5r—-3 x-1 (1-12 
(e) (16) partfrac( 28/((x*5)*(x^2431)) , x); 


a 5- x 
TUR xhi xb 

26 5-2 1 
So, 


(£ +5)(£? +1) 2? +1 PO 

In each case, the denominators of the partial fractions are the linear 
and quadratic factors of the denominator of the original expression, 
except that where the denominator of the original expression has a 
repeated factor, there is an additional partial fraction whose 
denominator is the square of this factor. 


You will learn more details about this in Unit 7. 
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Solution to Computer activity 41 


(a) (Aii) quotient (x*5+x+1, x%3+2); 
(%o1) x? 
(/12) remainder (x*5+x+1, x%3+2); 
(%02) -2x?*x41 
z?--r--1l —277+2+1 
So, OSEE =A o 
(b) (i) (413) quotient (3*x*4-7*x, x-1); 
(%03) 3x?*3x?43x-4 
(i4) remainder (3*x*4-7*x, x-1); 
(%04) -4 
So, an cda 
x—1l r—1l 


ii) (15) quotient(x^6*x^5b, 2*x*2+1); 
q 


4x^-Ax?-2x?-2x41 
(505). ——————— 
(i6) remainder (x^6+x^5, 2xx22+1); 
(06) E 
LO 
8 
qb gr Ag? + 4r? — 2x? — 2x +1 
So, = eee 
2z? + 1 8 
^ Ë zz? m 1 


Solution to Computer activity 42 


(2z — 1)/8 
22? +1 
2r—1 


8(2z? + 1) 


(a) (i) i1) eqn:'diff(y,x) - y = exp(x)*sin(x) ; 


d 
(%01) —y-y=he* sin(x) 
dx 


(ii) (%i2) ode2(eqn, y, x); 
(402) y=he* (hc-cos (x)) 


So the general solution of the differential equation is 


y = e? (c — cos(z)) , where c is an arbitrary 


constant. 


(b) (413) ode2(3*'diff(y,x) + y = x*y22, y, x); 


yes 


Yc= (-3x-9) 4e 3 


(414) fullratsimp(/); 


(403) y= 


(504) YT yes x8 or x8 


(Remember, % refers to the result of the most recently evaluated 


command.) 


So, the general solution of the differential equation is 


1 
Y ce yrr3 
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where c is an arbitrary constant. 


Solutions to Computer activities 


(c) (4i5) ode2('diff(y,x) = exp(cos(x))-1, y, x); 
(05) y= [o dx -x- c 


So, the general solution of the equation is 


y= | *àz- se 


where c is an arbitrary constant. 
Here, Maxima has calculated the solution in terms of an integral that 
cannot be evaluated analytically. 


(d) (416) ode2('diff(y,t) = try?3-1, y, t); 
(406) false 


'The ode2 command failed to find an analytic solution to this equation. 


Solution to Computer activity 43 
(a) (hi1) eqn:t*'diff(y,t)42*y-t^2; 
d 
%01) t| —y )+2y=t? 
(401) (E) y 


(b) (412) sol:ode2(eqn, y, t); 
4 


t 
a 
(402) y= zz 
t 
So the general solution is y = 4 mo where c is an arbitrary constant. 
. . i? [e 
'This can also be written as y — T + z 


(c) (413) ici(sol, y=1, t=1); 


t*43 
vA 3 = 
(403) y 1d 
"e i443 
So the required particular solution is y = E 
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Solutions to Computer activities 


Solution to Computer activity 44 

(a) First find the general solution of the equation. 
(111) eqn:'diff(ly,x)-y-(x-2)^2; 
(hoi) Z yy- (x-2)? 


(12) sol:ode2(eqn, y, x); 
(402) y=(-(-x?-2x-2) khe * +4 (-x-1) ke * +4 %e"* +%c) he* 


'Then find the particular solution required. 


(13) p:ici(sol, y=1.9, x=0); 

rat: replaced -0.1 by -1/10 = -0.1 
, 4e*-10x?*20x-20 

(403) yt — d. 


This solution can be simplified using the expand command. 
(414) p2:expand(p); 
(404) q +x2-2x+2 
So the required particular solution is 
e 


+r? Dn 


Y=- io 


(b) The expression to plot is the right-hand side of the equation 
represented by the variable p2. 


(415) wxplot2d(rhs(p2), [x,0,5], [y;0;5]]; 
5 T 


T T T 

a 

E xdi 

y 

NX 3r 

x 
C%t5) = 2L 

x 

$ tt 

x 

(%05) 


Here, a suitable range for y has been chosen. 


(c) To find the value of the solution when x = 4, substitute 4 for x in the 
right-hand side of p2 


(16) subst(4,x,rhs(p2)); 
het 


(%06) 10- i 


4 


So, when z =4,y=10—<. 
10 
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, 


false) 


5 


/ PROGRA-—1/MAXIMA-1.0/share/maxima/5.30.0/share/diffequations/drawdf.mac 


(Some output may appear here.) 


0101) € 
(b) (412) wxdrawdf(y-(x-2)^2, [x,y], field arrows 


Solution to Computer activity 45 


(a) (Aii) load(drawdf) 


5 


false) 


10 


[y,-5,10], field.arrows 


[x,-3,6] > 


1 
" 


VN N= 7101 


ME 


(413) wxdrawdf(y-(x-2)^2, 


(402) 


\ 
y 
\ 
\ 
y 
\ 
\ 
y 
\ 
1 
Hd 
\ 
\ 
y 
\ 
\ 
H 
i 


porron 
LAI“ 
Liceo 
NAL LLL ooo 
ANNA LL 1122 Aree ere a 
SASNAN LTDA AAP Re 0" - 
SSNS NC Mob tt 
SN CE ul cy te ee 
SARNA AAA gru 
SSAA SEN il LA RI 
SSSR NG oue et mom 
SSA SS NON CNUN NU p ume 
SS SSN RA LG Pe euo 
SS ANA YT ans me 
SSA RAN D RP est att AA 
AUSUNCRONCNC NS ou rac tmn 
CONNU oa ute tmt mm 
SNA E bao gata en 
NANA IAA ee eee ee He 
VV LIS Loco 
LISA se eesre ern nH HH HK 
a Aa a el | 
A eee eee) 
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(%03) 


IVities 


Solutions to Computer act 


From the direction field, the particular solution passing through 


the point (0,0) seems to be as shown below. 


(d) (i) 


Peeples eae ae Sed Se lees 
PP uf ah gh OEE A eae BE Br ae pe 
PPP PPP ese enn nH HHH HK 
Nill P eererer7r7-7----- H] 
NANA LIS AS 4 4 e227 ee em mem Re 
SANNA UN DAIS FS ORS Mo mK oe eee 
ESS ASSN N DL LDL Se eee ee eee 
"eom Dandi si at, NA qo gr urat a 
"wo. UON NON ACA Jur utres enm 
SS eS S URSUS RUN. poo a ur univ stet 
[+= e NN AA LEA eee 
—NNNNNMM LI ee 
Mee tots Ai ae a RP ae 
Pe SR RY d E Pee ee v on] 
SAA SSA SRN NE et 
SSAA LSANSARANNA LTD SA eee 
la ~ NOMCNONIN OM A qf Lot At 
NC SERNA t vu 
== ED A e 
SS S XE qi dew 
Es E geee 
SS 1111414 
NN VIP E 
Fg PO EE 
iF T E EEE 
r PAPES Ree AAA 
Eat [asisto a | elec eeu pene pee eoe epe 
un oq Y 


10 
8 
6 
4 


3 


the point (0,4) seems to be as shown below. 


(ii) From the direction field, the particular solution passing through 


Vea Sales de 2 Ru ewes 
Lorca 
Litro 
Rll P erererr--------- E 
NNNAV LDA AS ee ee eee eee em nme 
NO A p off abut! annm mtm m 
SS IL IDSA ee eee eee c 
SSUES PLA Ae aS ee ee 
O A gx HAE eo ee 
Momo OU RN LE oec ue ur nue tt 
[AN NNNNMLP gg 742777 
SRSA SASSANID DA SAAS Oe 
Mme ey iy e UN ANA o ue ven v 
[=== => ASNSNSANSN ADD PX eT 
SSSA SSSR A eee tu 
mee CK NN A gh a atv A 
SANNNNAA TD Dae eeerr ae 
SONS NV 11444090422" 
SNA VILLE 
= NVWL LAA 44 42447272 
^ WV M Patentes 
S Ll 10ccccccc==- 
\ lA ee 
H PPP 
/ Vn HESS GIU tero ve iaa as 
^ AAA AA A A A e A e 
AAN EA ESA ESAS AAA 
O © O FTN o N Y 
= 


false, 


[x,-3,6], I[y;-5,10], field.arrows 


(14) wxdrawdf(y-(x-2)^2, 


--A--Ad--4--243--1--1--I. 
ERRADA AAA 
Litro 
PULL ccoo 5 
NAAA LEI ccoo 
*NRCRONON T Ee) mmm mm m 
SS SENNA UL I SAA Oe ee ee ee el 
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false, 


[x,-3,6], [y,-5,10], field arrows 


(15) wxdrawdf(y-(x-2)^2, 


3 


7 


NN--d77 717 


[trajectory.at, 0,0], [trajectory.at, 0,4]) 
10 


NN et er 7 


Ls sr—=~NrNV 


1 


H 
FPP ROP ae aa XX, 


/ 
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/ 
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See Pe eee 
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(%05) 


false, 


[x,-3,6], [y,-5,10], field arrows 


(416) wxdrawdf(y-(x-2)^2, 


, 


LLANO 


REALI 


[trajectory.at, 0,0], [trajectory.at, 0,4], [trajectory.at, 0,1.9]) 
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Solutions to Computer activities 


Solution to Computer activity 46 
(It is assumed the drawdf package has already been loaded.) 


(a) The direction field, together with the approximate particular solution 
passing through (0,0), can be plotted as follows. 


(Ail) wxdrawdf(exp(cos(x))-1, field .arrows-false, [trajectory_at, 0,01) 
10 


(t1) 0 


(%01) 

For the particular solution displayed, y generally increases as x 

increases, but there is some oscillation about this trend. 

Every particular solution of this equation will display the same 

behaviour, since the direction field slopes do not vary as y changes. 
(b) The direction field, together with the approximate particular solutions 

passing through (0,0), (1, 1) and (—1, —1), can be plotted as follows. 


(412) wxdrawdf(t*y^3-1, [t,-2,2], [ly,-2,21, field_arrows=false, 
[trajectory.at [trajectory.at, 1,1], [trajectory.at, -1,-1]); 


PAPA 
EVA 


Wir Arr do E 
yA at A A A tersrs-r--- 


(162) 0 


o 
al 
ry 


^44 g44 


eee eer yy sss 


SS SR A AR SR SOR RR 


SSS sk Sa RA G 


p I 


coo» 


PIE 
Soe ee yy sss ss 
eee eee sss | 4 n4 4 
wee ere tI ttt 
== td 
a 1 Yr//y 


DAA 
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- 
a 
N 


(402) 

For the particular solution that passes through (0, 0), y decreases as t 
increases, with the gradient of the curve being steeper the larger the 
magnitude of t. 


The particular solution that passes through (1,1) has a local minimum 
near t = 1. The particular solution that passes through (—1, —1) has a 
local maximum near t = —1. 


All the particular solutions of this differential equation will have 
behaviours similar to one of the above. 
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Solution to Computer activity 47 


(a) 


(b) 


(511) pbsirk(y-(x-2)^29, y 1.9, lx 05.5, 1125 
(%01) E[[0.0,1.9],[1.0,0.70833333333333] , [2.0,1.18923611111111] , 
[3.0,2.783347800925924] , [4.0,3.97573362750771] , [5.0,0.66344524116672]] 


(412) asol:-exp(x)/10*x^2-2*x42; 
0j WX 


(02) JE 4x2-2x42 


(413) wxplot2d([asol, [discrete, pts]], [x,0,5], [legend, "Analytic", "Numerical"]); 
5 T 


4.5 
4 
3:5 


(%t3) 25 


(%03) 


(414) pts2:rk(y-(x-2)^2, y, 1.9, [x, 0, 5, 0:11)$ 

(Note that this line has been ended with $ to prevent the display of 

the calculated points.) 

(Aib) wxplot2d([asol, [discrete, pts2]], [x,0,5], [legend, "Analytic", "Numerical"]); 
5 T 


4.5 
4 
3.5 


(%t5) 25 


(%05) 
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Solution to Computer activity 48 


(a) 


o 
= 


154 


(Ail) pts:rk(exp(cos(x))-1, y, 0, [x,0,10,0.1])$ 
(412) wxplot2d([discrete, pts]); 
3.5 


3 L 4 
25+ J 
2 L 4 
> 
(462) 15+ | 


(402) 
This solution is consistent with the direction field that you found in 
Computer activity 46(a). 


(413) pte:rk(6*y^3-1, y, 1.5, [t,0,1,0.01])$ 


Note that here the independent variable is t. A step size of 0.01 has 
been used to ensure that a smooth curve is obtained. 


'The approximate solution can be plotted, with suitable axis labels and 
vertical range, as follows. 


(14) wxplot2d([discrete, pts],[xlabel,"t"], [y,0,2.5]); 
2.5 


2 L x 


aa ^ al | 


(04) 
This solution is consistent with the direction field that you found in 
Computer activity 46(b). 


Solutions to Computer activities 


Solution to Computer activity 49 
(a) Du Praeteix( 12,383,315 21,21, [3,91,4] 25 


2: 3. 3 
(%01) |-2 1 2 
1 -1 4 


(b) (412) Q: matrix( 
[5,-3] , 


(o) (018) Rmatrix DH 2,1), 14, 7,0]1, [55, 1,315; 


i 3. 
(403) |4 7 0 
"b, d 3 


Solution to Computer activity 50 
(a) (4i4) P+R; 


3 =5. 2 
(404) | 2 2 
-4 7 


8 
0 
3 —5 2 
So P+R = 2 8 2|. 
(b) (%i5) 3*P; 
6 -9 9 
(%05) |-6 3 6 
3 -3 12 
6 —9 9 
So 3P = | —6 3 6]. 
3 —3 12 


(c) (6) RaQ; 
7 -3 
(%06) | 13 -26 
-26 25 
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(di QUT Poem 
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NS 


— 


== 


13 -12 12 
(407) |-4 5 4 
8 -8 17 
13 -12 12 
SoP?=| -4 5 4 
8 —8 17 
(/18) determinant (P); 
(408) -15 
So det P — —15. 
(%i9) invert(R); 
5 7 
2 6 6 
1 2 
(409) |-2 -2 -- 
3 3 
13 3 5 
2 2 2 
7 5 7 
2 6 6 
SoR-!- 2 z 2 
3 3 
13 3 5 
2 2 2 


Gato) z:matriz( [0]; H], [i]: 
[°] 
(%010) |1 


Cali) Pixs 


0 
(aii) 13 
3| 
0 0 
Since Px = | 3 | = 3x, where x = | 1 |, it follows that x is an 
3 1 


eigenvector of P and the corresponding eigenvalue is 3. 


Solutions to Computer activities 


Solution to Computer activity 51 
(à) G) ON Axmaters(CIO,1, 11; E-1,2,1], -1,-1,405; 


(401) ^ 2 1 


(ii) (412) charpoly(A,m) ; 
(402) -((2-m) (4-m) *1)m-2m*6 


(iii) (413) expand CD ; 
(%03) -m?46m?-11m*6 
So the characteristic polynomial of A is =m? + 6m? — 11m +6 
and the characteristic equation is =m? + 6m? — 11m +6 = 0. 
'This could be expressed in terms of A as 


—M + 6A? — 111 4- 6 2 0. 


(b) (ae) B:matrix( [4,0,0]; [2,1,-3],5 [2,-3, 102 ; 


4 0 0 
(404) |3 1 -3 
3 -3 1 


(Aib) expand(charpoly(B,m)); 

(%05) -m?*6m?-32 

So the characteristic polynomial of B is —m? + 6m? — 32 and its 
characteristic equation is =m? + 6m? — 32 = 0. This could be 
expressed in terms of A as 


=) 26 ag eq 


Solution to Computer activity 52 
(a) (416) eigenvalues(A); 
(406) [[1,2,3] , [1,1,11] 
So the eigenvalues of A are 1, 2 and 3, each of which occurs once. 
(b) (417) eigenvalues (B); 
COT) [[-2,4] , [1 ,21] 
So the eigenvalues of B are —2 (which occurs once) and 4 (which has 
multiplicity 2). 


157 


Solutions to Computer activities 


Solution to Computer activity 53 
(a) (418) eigenvectors(A); 
(%08) [[[1,2,3],[1,1,1]], [E1545] s Ct 51,017; 001,1,2111] 


So the eigenvalues of A are: 


1 
e 1, with corresponding eigenvector i : 
1 
2 
1 
e 2, with corresponding eigenvector | 1 |, 
il 
1 
e 3, with corresponding eigenvector | 1 
2 


These are the same as the eigenvectors verified in Activity 9 of 
Unit 11, except that the eigenvector corresponding to the eigenvalue 1 
is half of that given in Activity 9. 
(b) (419) eigenvectors (B); 
(109) LLE-2,4] ,[1,2]],[[[0,1,1]],[[1,0,1],[0,1,-11]]] 


So the eigenvalues of B are: 


0 
e  —2, with corresponding eigenvector | 1 |, 
1 
e 4 (which has multiplicity 2), with corresponding eigenvectors 
1 0 
QO] and | 1 
1 —1 


Solution to Computer activity 54 
(a) G) ui t8; 
(%01) 5 
(412) x[n] :=2*x [n-1] *3; 
(402) Xp :=2Xp-1+3 


Gy 4359) x[2]; 
(403) 13 
So x» = 13. 
(111) (414) x[10]; 
(404) 4093 
So x19 = 4093. 
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(b) G) i5) u[01:4; 
(405) 4 
(416) uli]: 
(406) 9 
(17) ulm] :=5*u [n-1] -6*u [n-2] ; 
(LoT) u,:75u41-6u4-2 


Ko) 


d 


(ii) (%i8) u[4]; 
(408) 129 
(1419) u[10]; 
(409) 62121 
So U4 = 129 and uio = 62121. 


Solution to Computer activity 55 


(a) (Aii) load(solve.rec); 
(401) C:/PROGRA~1/MAXIMA~1.0/share/maxima/5.30.0/share/solve_rec/solve_rec.mac 


(b) Q) (%i2) killGO; 
(402) done 


(ii) (413) solve .rec(x[n]-2*x[n-1]*3, x[n], x[1]=5); 
(403) x,=2”*?-3 
So the closed form for the sequence given by the recurrence 
system is £n = 2^*? — 3. 
(c) (1414) kill(u); 
(%04) done 
(Aib) solve.rec(u[n]-5*u[n-1]-6*u[n-2], u[n], u[0]24, u[1129); 
(405) u,73**3 2” 
(Note the space between 3 and 2”, which indicates multiplication. 
wxMaxima does not display a multiplication sign.) So the closed form 
for the sequence given by the recurrence system is Un = 3" +3 x 2”. 
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Solution to Computer activity 56 
(It is assumed the solve. rec package has previously been loaded.) 
(a) (Aii) solve_rec(x[n]=2*x[n-1]-1, x[n], x[1]=3); 
(fol) x,=2"+1 
So the closed form for the sequence given by the recurrence system is 


Ln = 2” +1. This is the same as the solution obtained in Unit 12, 
Activity 20(a). 


(b) (%i2) solve_rec(y[n]=-1/3*y[n-1]+1, yin], y[1]=12); 
2-n = n 
(402) y Ca 43 * C1) 87 


This result can be simplified as follows. 
(413) fullratsimp(%) ; 
4152135 (245 
YA 3 = 
(103) y. fag 


So the closed form for the sequence given by the recurrence system is 
| 3ntt 135(-1)” 
Hus 4 x 3n 
g^ 135 , (-1)” 
x 3n 4 3” 


~ 4 
29 4. < EU" 
4 4x3 gn=1 
3 de q A 
zc 


This is the same as the solution obtained in Unit 12, Activity 20(b). 


(hi4) solve.rec(u[n]27*u[n-1]-12*u[n-2], ulm], u[0]22, u[1127); 
(404) u,24* 43^ 

So the closed form for the sequence given by the recurrence system is 

Un = 4" +3”. This is the same as the solution obtained in Unit 12, 
Activity 30(a). 


(/15) solve.rec(u[n]26*u[n-1]-9*u[n-2], u[n], u[0]22, u[1]=7); 
(%o5) u,=(5+2) 3" 

So the closed form for the sequence given by the recurrence system is 

Un = 4 +2) x 3”. This is the same as the solution obtained in 

Unit 12, Activity 30(b). 
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Mathematical operations and functions 


Operation or function 


Addition 

Subtraction 
Multiplication 
Division 

Brackets 

Powers, for example 2? 


Square root, for example v5 


Exponential, for example e? 


Natural logarithm, In 


Magnitude (absolute value) 


sinh 
cosh 
tanh 
cosech 
sech 
coth 
sinh“! 
cosh~! 
tanh“! 


Equality (within an equation) 


Syntax 


ct 
em, 
wert ONES ONES ONES we 


Example 


4+3; 
4-3; 
4*3; 
4/3; 

2* (3+4) ; 
2^3; 
2**3; 
sqrt(5); 


%e^3; 
exp(3); 
log(8); 


abs (-3); 


sin(1); 


cos (3*4pi/2); 


tan(/pi/4); 
csc(2); 

sec (/pi/3); 
cot (/pi/4); 


asin(sqrt(3)/2); 
acos(1/sqrt(2)); 


atan(1/2); 


sinh(0); 
cosh(0) ; 
tanh(0); 
csch(1) ; 
sech(0); 
coth(1); 
asinh(0); 
acosh(1); 
atanh(0); 


5-2*x*3; 


Maxima reference guide 
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Mathematical operations and functions (continued) 


Operation or function 


Matrix addition 

Matrix subtraction 

Scalar multiplication 

Matrix multiplication 

Matrix powers 

Determinant (of a square matrix) 
Matrix inverse (of a square matrix) 


Factorial 


TOR 


vp. 


Mathematical constants 


€ 


he?2; 
T Qx*hpi; 
i i 2+3*%i; 
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Syntax 


+ 


* 
^^ 
determinant (matrix) 


invert (matrix) 
or matrix^^(-1) 


factorial( ) 
or o! 


binomial (in, xr) 
or combination(n,r) 


permutation(n,r) 


Note: the combination and 
permutation functions require 


the functs package to be loaded. 


A.B; 

A^^3; 
determinant (A); 
invert (A); 
A^^(=1); 


factorial(5); 
5!; 


binomial(5,2); 
combination(5,2); 


permutation(5,2); 


Maxima commands 
Operation 


Get help on a command 


Find help information whose 
title contains the given text 


Prevent the evaluation of a 
command 
Force the evaluation of a command 


Load a package 


Assign a value to a variable 

Display the value of a variable 

Define a function 

Evaluate a function at a value 

Remove an assigned variable 
or function 

Remove all assigned variables 
and functions 

Reset all system variables 


Convert to a decimal number 


Substitute 


Expand brackets 

Factorise 

Simplify 

Simplify something involving 
exponentials and logarithms 

Combine logarithms 

Expand trigonometric functions 
of sums, differences and 
multiples of angles 

Express powers of sines and 
cosines in terms of sines and 
cosines of multiple angles 

Simplify trigonometric expressions 

Simplify algebraic fractions 
containing trigonometric 
functions 

Find the partial fraction expansion 
of a rational expression 


Command 


? command 
or describe (command ) 
27 


or describe( , inexact) 


BEd D) 


load (package name) 


variable 


function ( ) 
kill (variable ) 
kill (function ) 
kill (all) 


reset () 


float (D 


subst (value , variable , 
expression ) 

expand( ) 

factor( ) 

fullratsimp( ) 

radcan( ) 


logcontract( ) 
trigexpand( ) 


trigreduce( ) 


trigsimp( ) 
trigrat( ) 


partfrac (expression , 
variable ) 


Maxima reference guide 


Example 


? float; 
describe(float); 

7T? flo; 

describe(flo, inexact); 


'diff(y,x) 


gGO:-2''(diff(x^4,x)); 
load(implicit. plot); 


a:23; 

a; 

f (x) :=2*x+3; 
£(1); 
kill(a); 
kill(f); 
kill(all); 


reset(); 


float(sqrt(2)); 


subst(4, x, x^2*1); 

which substitutes 4 for x in x^2*1. 
expand( (x*1)^2 ); 
factor( 2*xt+4*x%2 ); 
fullratsimp( (2*x+4*x%2)/x ); 
radcan( log(x^2) ); 


logcontract (log(a)+log(b)) ; 
trigexpand(sin(A+B)) ; 


trigreduce(sin(x)^2); 


trigsimp(sin(x)^2*cos(x)^2); 
trigrat((sin(2x))/cos(x)); 


partfrac(l/(x^2-1),3); 
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Maxima commands (continued) 


Operation 


Quotient on dividing one number or 
polynomial expression by another 
Remainder on dividing one number or 
polynomial expression by another 
Highest common factor (greatest 
common divisor) of two numbers 
Bézout's identity: 
given integers a and 6 find the 
highest common factor d 
and integers v and w such that 
av + bw — d 


Left-hand side of an equation 
Right-hand side of an equation 


Solve an equation, exactly 
Solve an equation for a variable 
Solve simultaneous equations 


Solve an equation numerically 


Extract an element of a list 


Specify a matrix 

Find the characteristic polynomial of 
a matrix 

Find the eigenvalues of a matrix 

Find the eigenvectors of a matrix 
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Command 


quotient( , m) 
remainder( , ) 
gcdC , ) 


gcdex( , ) 


lhs ( equation ) 
rhs ( equation ) 


solve (equation ) 
solve(equation, variable) 
solve (list of equations, 
list of variables) 
find root (expression , 
variable, 
interval start value, 
interval end value) 


list [index] 


matrix(row,row,...) 


charpoly (matrix , variable ) 


eigenvalues (matrix) 
eigenvectors (matrix) 


Example 


quotient (6,4) ; 
quotient (x*2,x-1); 
remainder (6,4); 
remainder (x*2,x-1); 
gcd(6,4); 


gcdex (6,4) ; 


lhs (4*x+1=2*x-2) ; 
rhs (4*x+1=2*x-2) ; 


solve (2*x*%2-1=0) ; 

solve (2*a*b-3*b=0, a); 

solve([2*x+y=4,x-2*y=1], 
Ei 

find. root(2*xx^2-1,x,0,1); 


A[2] ; 


A:matrix([1,2],[3,4]) ; 
charpoly(A,m) ; 


eigenvalues (A); 
eigenvectors (A); 


Maxima commands (continued) 


Operation 


Differentiate 
Differentiate multiple times 


Integrate 
(find an antiderivative) 
Evaluate a definite integral 


Find an approximate value 
of a definite integral 


Make an assumption about a 
variable 

State two things are equal 

State two things are not 
equal 

Forget a property 

List all known facts 

List all known facts about 
a particular variable 


Solve a differential equation 
analytically 


Apply initial conditions to 
the analytic solution of a 
differential equation 


Plot the direction field of a 
differential equation 


Solve a differential equation 
numerically 


Command 


diff( 
diff( 


, variable ) 

, variable , positive integer ) 

integrate( , variable) 

integrate ( 
lower limit , upper limit ) 


quad_qags( , 


lower limit , upper limit ) 


, variable , 


variable S, 


assume( ) 


equal( , ) 
notequal (=,=) 
forget (property) 
facts() 


facts (variable ) 


ode2 ('equation , dependent variable, 
independent variable ) 


ici (solution , variable =value , 
variable = value ) 


wxdrawdf (expression, 
[variable , variable], 
field arrows-false, ...) 

Note: this command requires the 
drawdf package to be loaded. 


rk( expression , dependent variable, 
initial value, 
[independent variable , 
initial value, final value, 


step size ]) 
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Example 


diff(sin(x^2),x); 
diff(logGO,x,3); 


integrate(x^2,x); 
integrate(sin(x),x,0,1); 


quad.qags (exp(-x^2) ,x,0,1); 


assume (a»0); 


assume (equal(n,3)); 
assume (notequal(n,-1)); 


forget (a>0) ; 
facts(); 
facts(a); 


ode2('diff(y,x)=x, y, x); 


ici(sol, y=1, x=0); 


wxdrawdf (x+y, [x,y], 
field_arrows=false) ; 


rk(x^2,y,1,[x,0,10,0.1]) ; 
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Maxima commands (continued) 


Operation 


Plot a graph 


Plot several graphs 


Plot a curve from 
parametric equations 


Plot a graph represented 
by a set of points 


Plot a curve represented 
by an equation in 
implicit form 


Plot several curves 
represented by 
equations in implicit 
form 


Set the default plotting 
options 
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wxplot2d (expression , 


horizontal range,... 


wxplot2d(list of expressions , 


horizontal range,.. 


wxplot2d([parametric, 
expression for r, 
expression for y, 


.) 


[parameter , min, [max]], 


wxplot2d([discrete, 
list of points]) 


(use plot2d similarly if not 


using wxMaxima) 


wximplicit_plot (equation , 


horizontal range , 
vertical range,...) 


wximplicit_plot ( 
list of equations , 
horizontal range, 
vertical range ,...) 


(use implicit plot similarly if not 
using wxMaxima) 


Note: the wximplicit plot and 
implicit.plot commands require 
implicit.plot the package to be 


loaded. 


Set plot option( option) 


Example 


wxplot2d(x^2, [x,0,1]); 
which plots the graph of x? 
for x between 0 and 1. 
wxplot2d(x*2, [x,0,1],ly,0,2]); 
which plots the graph of x? 
for z between 0 and 1 and 
vertical axis values between 0 and 2. 


wxplot2d([x^2, 2*x], [x,0,1]); 
which plots graphs of z? and 2x 
for x between 0 and 1. 


wxplot2d([parametric, 
2*t, 3*t*1, [t,0,1]]); 


wxplot2d([discrete, 
[[0,0], [1,1]] 1); 


wximplicit.plot(x^2-*y^2-1, 
[x,-1,1], l[y,-1,1]); 


wximplicit_plot ( 
[x^2*2*y^2-4, y-x^3], 
[x,-1,1], [y,-1,1]); 


set_plot_option([gnuplot_preamble, 
"set size ratio -1"]); 


Maxima commands (continued) 


Operation 


Define a recurrence sequence 


Calculate a term in a 
recurrence sequence 


Solve a recurrence system 


Command 


sequence [initial term number]: 
initial value 

sequence [n] :=expression 
containing previous terms 


sequence [term number ] 


solve_rec(recurrence relation, 


general term, initial term(s) ) 


Note: this command requires the 
solve_rec package to be loaded. 


Maxima system variables 


Variable 


fpprintprec 


functions 
values 


Description 


The number of significant figures of 
a decimal number to display 
The list of all user defined functions 


Example 


The list of all user assigned variables | values; 


Maxima graph plotting options 


fpprintprec:3; 


functions; 
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Example 


x[1]:1; 


x[n] : 2*x [n- 1] +x [n-2] ; 


x [50] ; 


solve rec(a[n]-2*a[n-1], 
a[n] ,a[1]721) ; 


'The following options can be added as additional arguments to plotting 


commands. 


Option 


Legend text 


scales.) 


Vertical range 
Curve colour 


Turn legend off 
Horizontal axis label 
Vertical axis label 


Argument 


(y, B. li 
[color, | ,...] 
[legend," ",...] 
[legend, false 
[xlabel, " "] 
[ylabel, " "] 


Change the relative lengths/scales | [gnuplot preamble, 
of graph axes. (A positive 
number specifies the ratio of 
axis lengths, a negative number 
specifies the ratio of axis 


Change the number of points [nticks,  ] 
plotted on a curve from 
parametric equations 


Example 


[y, 0, 5] 

[color, red] 
[legend, "Car A"] 
[legend, false] 
[xlabel, "t (h)"] 
[ylabel, "s (km)"] 
[gnuplot.preamble, 


"Set size ratio number"] "set size ratio -1"] 


which sets the axes to 
have equal scales. 


[nticks, 100] 
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Maxima direction field plotting options 


Operation Argument Example 
Set the range of a variable neag, NN. NI [x,0,5] 
Plot the graph of a particular solution | [trajectory_at, F, l| [trajectory.at, 0, 1] 


Set the maximum step size to use when | [tstep,  ] [tstep, 0.001] 
approximating a solution 


Set the number of steps used when [nsteps,  ] [nsteps, 200] 
approximating a solution 


Maxima packages 


Maxima package | Functions added 

drawdf Plotting direction fields of differential equations 

functs Additional mathematical functions including combinations and permutations 
implicit plot | Plotting curves given by equations in implicit form 

solve_rec Solving recurrence systems 


wxMaxima keyboard sequences (Microsoft Windows) 


Key sequence Function 
Zoom in Alt-I 
Zoom out A1t-0 
Copy Ctrl-C 
Interrupt a calculation Ctrl-G 
Re-evaluate the worksheet | Ctr1-R 
Close wxMaxima Ctrl-Q 
Paste Ctrl-V 
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wxMaxima menu commands (Microsoft Windows) 


Operation Menu name | Menu option 
Close wxMaxima File Exit 

Save your work File Save or Save As 
Print the worksheet File Print 

Open a worksheet File Open 

Configure wxMaxima | Edit Configure 

Copy Edit Copy 

Paste Edit Paste 

Interrupt a calculation | Maxima Interrupt 


Restart Maxima Maxima 
Maxima help window | Help 
Re-evaluate worksheet | Cell 


Insert text cell Cell 
Insert title cell Cell 
Insert section cell Cell 


Insert subsection cell | Cell 
Enter a matrix Algebra 


Restart Maxima 

Maxima Help 

Evaluate All Visible Cells 
Insert Text cell 

Insert Title cell 

Insert Section cell 

Insert Subsection cell 
Enter Matrix... 


Maxima command-line interface commands 


Operation 


Close Maxima 

Interrupt a calculation 

'Turn off 2D formatting of maths 

Display output left justified 

Start writing a transcript 

Stop writing a transcript 

Playback all the input and 
output in a session 

Playback a range of input and 
output 

Save the state of a session 

Reload a session 


Command 


quit O; 

Ctrl-C 
display2d:false; 
leftjust:true; 
writefile("file path"); 
closefile(); 
playback() ; 


playback( [start line number, end line number] ); 


save("file path", all); 
load(" file path"); 
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Index 


abort calculation 
using command-line interface 106 
using wxMaxima 20 
abs 15 
absolute value 15 
accessibility 95 
changing wxMaxima colours and fonts 95 
using a screen magnifier 96 
using a screenreader 96 
using wxMaxima with keyboard alone 95 
acos 15 
acosh 161 
adding text to a wxMaxima worksheet 27 
addition 15 
of matrices 85 
annotating a wxMaxima worksheet 27 
approximation, numerical 
of a definite integral 57 
of a solution to a differential equation 79 
of a solution to an equation 49 
argument of a command 15 


asin 15 
asinh 161 
assign 


a derivative or integral to a function 54 

function 32 

variable 23 
assume 56 
assumption 

about a variable 56 

equals 56 

forget 56 

list all 56 

not equals 56 
atan 15 
atanh 161 


Bézout’s identity 59 
binomial 162 
binomial coefficients 162 
brackets 15 

expanding 35 


CAS see computer algebra system 
cell, in wxMaxima 18 

marker 11, 18 

section 27 

subsection 27 

text 27 

title 27 


characteristic equation of a matrix 86 
characteristic polynomial of a matrix 86 
charpoly 87 
closefile 115 
closing 
command-line interface 101 
wxMaxima 13 
colour of graph 39 
combination 162 
command 
force evaluation of 55 
prevent evaluation of 70 
command-line interface 8, 98 
abort calculation 106 
closing 101 
creating transcript 112 
cursor, editing 99 
editing commands 106 
editing cursor 99 
graphs 117 
help 115 
input line number 98 
interrupt calculation 106 
output line number 99 
plotting graphs 117 
printing 112 
reusing commands 106 
saving work 112 
scientific notation 105 
transcript 112 
computer algebra system 5 
configure wxMaxima 13, 18 
copy, using wxMaxima 21 
cos 15 
cos 15 
cosec 15 
cosech 161 
cosh 161 
cosh ! 161 
cot 15 
coth 161 
csc 15 
csch 161 
cursor 
editing 
in command-line interface 99 
in wxMaxima 21 
horizontal, in wxMaxima 21 


Index 


171 


Index 


decimal approximation of a number 19 equation 45 
define changing the subject of 47 
differential equations 70 differential 70 
functions 32 direction field 74 
matrices 82 initial condition 72 
recurrence systems 91 inputting 70 
variables 23 particular solution of 72 
definite integral 51 solving analytically 70 
approximate value of 57 solving numerically 79 
delete left-hand side 48 
function 34 parametric, plotting 60 
recurrence system 92 plotting 37 
variable 25 rearranging 47 
derivative 51 right-hand side 48 
assign to a function 54 simultaneous, solving 48 
describe 32 solving 45 
determinant 85 solving numerically 49 
determinant of a matrix 85 erf 53 
diff 51 evaluate 
differential equation 70 an expression 12 
direction field 74 force command to 55 
initial condition 72 prevent command from evaluating 70 
inputting 70 re-evaluate a cell using wxMaxima 21 
particular solution of 72 re-evaluate worksheet using wxMaxima 22 
solving exp 15 
analytically 70 expand 35 
numerically 79 expanding brackets 35 
differentiate 51 exponentials 15 
direction field of a differential equation, plotting 74 expressions 
division 15 evaluating 12 
of polynomial expressions 69 logarithmic and exponential, simplifying 35 
drawdf simplifying 35 
command 74 trigonometric, simplifying 35 
package 74 
fact see assumption 
editing commands factor 35 
using command-line interface 106 factorial 162 
using wxMaxima 20 factorising 35 
editing cursor facts 56 
using command-line interface 99 file types, wxMaxima 28 
using wxMaxima 21 find.root 49 
eigenvalues 88 float 19 
eigenvalues of a matrix 88 force evaluation of a command 55 
multiplicity 88 forget 56 
eigenvectors 89 forget assumption 56 
eigenvectors of a matrix 89 fpprintprec 26 
element of a list 48 front-ends 7 
equal 56 fullratsimp 35 
equals functions 34 


assumption 56 
within an equation 45 
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functions 32 


assign 32 
define 32 
delete 34 


hyperbolic 161 

inverse hyperbolic 161 
inverse trigonometric 15 
plotting 37 
trigonometric 15 


gamma_incomplete 53 

gcd 58 

gcdex 59 

gnuplot 43 

graphs 37 
changing default behaviour 43 
changing line thickness 97 
changing properties 39, 97 
colours 39 
equal scales 42 
from a set of points 81 
from equations in implicit form 41 
from parametric equations 60 
legend 39 
plotting 37 
using command-line interface 117 
vertical range 38 

greatest common divisor 58 


help 
using command-line interface 115 
using wxMaxima 30 
highest common factor 58 
horizontal cursor, in wxMaxima 21 
hyperbolic functions 161 
inverse 161 


ici 72 
implicit form of an equation, plotting 41 
implicit plot 
command 41 
package 41 
indefinite integral 51 
initial condition, for a differential equation 72 
initial value problem 72 
input line number 12 
input prompt, wxMaxima 11 
inputting 
differential equations 70 
matrices 82 
recurrence systems 91 


Index 


installing Maxima 9 


integral 
assign to a function 54 
definite 51 


approximate value of 57 
indefinite 51 
integrate 51 
interfaces 7 
command-line 8, 98 
wxMaxima 8,9 
interrupt calculation 
using command-line interface 106 
using wxMaxima 20 
inverse hyperbolic functions 161 
inverse of a matrix 85 
inverse trigonometric functions 15 
invert 85 


kill 
function 34 
recurrence system 92 
variable 25 

Kutta, Martin Wilhelm 79 


left-hand side, of an equation 48 
legend, of a graph 39 


lhs 48 
line number 
input 12 
output 12 
Lissajous curve 63 
list 48 


element of 48 
of assumptions 56 
load 
package 41 
previous session, using command-line interface 115 
loading packages 41 
log 15 
logarithmic and exponential expressions 
manipulating 35 
simplifying 35 
logarithms 15 
logcontract 35 


magnitude 15 
manipulating 
expressions 35 
logarithmic and exponential expressions 35 
trigonometric expressions 35 
mathematical operations 15 
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Index 


matrix 82 partfrac 67 
addition of 85 partial fraction expansion 67 
characteristic equation of 86 particular solution, of a differential equation 72 
characteristic polynomial of 86 paste, using wxMaxima 21 
determinant of 85 permutation 162 
eigenvalues of 88 playback 115 
eigenvectors of 89 plot2d 37 
input 82 using command-line interface 117 
inverse of 85 plotting 37 
multiplicaton of 85 changing default behaviour 43 
multiplicity of eigenvalue 88 changing line thickness 97 
powers of 85 changing properties 39, 97 
scalar multiplication of 85 colours 39 
subtraction of 85 direction field 74 
matrix 82 equations in implicit form 41 
Maxima 5 graphs 37 
instaling 9 legend 39 
menus, wxMaxima 10 parametric equations 60 
multiplication 15 number of points used 63 
matrix 85 set of points 81 
scalar and matrix 85 using command-line interface 117 
multiplicity of an eigenvalue of a matrix 88 using equal scales 42 
vertical range 38 
natural logarithm 15 powers 15 
not equals, assumption 56 of a matrix 85 
notequal 56 prevent evaluation of a command 70 
nsteps 78 printing 
nticks 63 using command-line interface 112 
numerical approximation using wxMaxima 29 
of a definite integral 57 prompt, wxMaxima 11 
of a solution of a differential equation 79 
of a solution of an equation 49 quad.qags 57 
numerical solution quadrature 57 
of a differential equation 79 quotient 
of an equation 49 on dividing integers 58 
on dividing polynomial expressions 69 
ode2 70 quotient 58, 69 
ordinary differential equation see differential equation 
oudrawdf 74 radcan 35 
output line number 12 re-evaluate 
ouwxdrawdf 74 a cell using wxMaxima 21 
a worksheet using wxMaxima 22 
packages 41 rearranging an equation 47 
drawdf 74 recurrence systems 90 
functs 168 changing 92 
implicit.plot 41 defining 91 
loading 41 deleting 92 
solve.rec 92 solving 92 
parametric equations remainder 
plotting 60 on dividing integers 58 
number of points used 63 on dividing polynomial expressions 69 
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remainder 58, 69 

reset 26 

resetting system variables 26 

reusing command 
using command-line interface 106 
using wxMaxima 20 


rhs 48 
right-hand side, of an equation 48 
rk 79 
roots 
higher 17 
square 15 


Runge, Carl 79 
Runge-Kutta 79 


save 115 
saving your work 
using command-line interface 112 
using wxMaxima 28 
scalar multiplication of a matrix 85 
scales, on graph 42 
Schelter, William 6 
scientific notation 
using command-line interface 105 
using wxMaxima 19 
screen magnifier, use with wxMaxima 96 
screenreader 96 
screenshot 30 


sec 15 
sech 161 
section cell, wxMaxima 27 
sequences 
given by recurrence systems 90 
term of 92 


set_plot_option 43 
significant figures 26 
simplifying 
expressions 35 
logarithmic and exponential expressions 35 
trigonometric expressions 35 
simultaneous equations 48 
sin 15 
sin ! 15 
sinh 161 
sinh"! 161 
solve 45 
Solve.rec 
command 92 
package 92 


solving 
differential equation 
analytically 70 
numerically 79 
equations 45 
numerically 49 
simultaneous 48 
recurrence systems 92 

sqrt 15 

square root 15 

subject of equation, changing 47 

subsection cell, wxMaxima 27 

subst 34 

substitute 34 

subtraction 15 
of matrices 85 

syntax 14 

system variables 26 
resetting 26 


tan 15 

tan? 15 

tanh 161 

tanh ! 161 

term ofa sequence 92 
text cell, wxMaxima 27 
title cell, wxMaxima 27 
toolbar, wxMaxima 10 
trajectory_at 78 


Index 


transcript, creating using the command-line interface 


112 
trigexpand 35 
trigonometric expressions 

manipulating 35 
simplifying 35 
trigonometric functions 15 

inverse 15 
trigrat 35 
trigreduce 35 
trigsimp 35 
tstep 78 


values 25 

variable 23 
assign 23 
assumption 56 
define 23 
delete 25 
system 26 
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Index 


worksheet, wxMaxima 9 
annotate 27 
re-evaluate 22 

writefile 115 

wxdrawdf 74 

wximplicit_plot 41 

wxMaxima 8, 9 
abort calculation 20 
adding text 27 
annotating 27 


cell 18 

marker 11, 18 
section 27 
subsection 27 
text 27 

title 27 
changing colours and fonts 95 
closing 13 
configure 13, 18 
copy 21 
cursor 

editing 21 


horizontal 21 
editing commands 20 
editing cursor 21 
file types 28 
help 30 
horizontal cursor 21 
input prompt 11 
interrupt calculation 20 
menus 10 
paste 21 
printing 29 
prompt, input 11 
re-evaluate 

cell 21 

worksheet 22 
reusing commands 20 
saving work 28 
scientific notation 19 
toolbar 10 
use with keyboard alone 95 
worksheet 9 

wxplot2d 37 
for plotting parametric equations 60 
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